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Wenn man für wirbelfreie Unterschallströmungen einer kompressiblen Flüssigkeit die Hodographen- 
methode verwendet, muß man die Theorie der analytischen Funktionen einer komplexen Veränderlichen auf den 
Fall allgemeinerer linearer partieller Differentialgleichungen vom elliptischen Typus erweitern. Für diesen 
Zweck werden in der vorliegenden Arbeit Integraloperatoren eingeführt. Sie transformieren die Klasse der 
analytischen Funktionen einer komplexen Veränderlichen in die Klasse der Stromfunktionen, welche die 
Unterschallströmungen der idealen Gase um beliebige Profile beschreiben. Die Operatorenmethode wird 
auf den gemischten Fall ausgedehnt. 


By using the hodograph method for irrotational compressible subsonic flow problems, the theory of ana- 

E Iytie functions of a complex variable has to be generalized to the case of more general linear partial differential 

equations of elliptic type. For this purpose, this note deals with integral operators. T’'hese operators transform 

j the class of the analytic funktions of a complex variable into the class of the stream functions representing 

subsonic flows of ideal gases around profils of arbitrary shape. The operator-method is extended to the 
mized case. 


Si l’on emploie la methode de l’hodographe pour les Ecoulements irrotationnels des fluides compressibles 
@ vilesse subsonique, il faut generaliser la theorie des fonctions analytiques d’une variable complexe au cas 
d’equations aux derivees partielles lintaires plus generales du type elliptique. Pour ce but, il est introduit, 
dans cette note, des op£erateurs integrals. Ces operateurs iransforment la classe des fonctions analyliques 
d’une variable complexe dans la classe des fonctions de courant qui decrivent les &coulements des gaz parfaits 
autour de profils arbitraires a vitesse subsonique. 


Ilpu ynorpeönennuu MmeTona ronorpada Aa ÖesBuxXpeBbIX 'I03BYKOBEX TeyeHHÄ C3KHMaeMoH 

3KHIKOCTH HeOOXOAHMO PACHPOCTPAHHTL TeOPHIO AHAJIHTUYeCKUX PYHKIMÄ KOMILIEKCHOTO TTepe- 

MeHHOTO Ha cIy4ah Öollee OÖINAX JIHHeÄHBIX AHbhepeHumassHsIX YPaBHeHHÄ B YaCTHEX IPOu3- 

BOAHEIX SLNHUTHYeCKOTO THUNa. Ü 9ToU MebIo B IIpenaraeMoH CTATbe BBOAATCA HUHTETPAIBHLIE 

Ä 'onepaTops. OHu TpauchopmupymT Kacc AHANHTUYecKuUX BPYHKIMÜ KOMILIeKCHOTO TepeMeH- 

. HOIO B Kllacc PYHKIIHÜ IOTOKA, ONHCBIBAWIIHX N03ByKoBhle TeyeHHA HNEAJIbHBIX  TA80B BOKPYT 
NPOH3BOABHBIX Ipobmseft. Meron omepaTopoB PacnpocTpaHseTca Ha CMeIMaHHBIHa cıyuah. 


sl. Partielle Differentialgleichungen in der Theorie zweidimensionaler Gasströmungen 


Die Tatsache, daß die mathematische Theorie stationärer zweidimensionaler Strömungen 
- — — jnkompressibler Flüssigkeiten im wesentlichen identisch ist mit der Untersuchung gewisser 
Probleme der Funktionentheorie, ermöglicht die Anwendung der hochentwickelten Mittel dieser 
Theorie. Dies ist einer der Hauptgründe für die beträchtlichen Erfolge der Theorie der in- 
kompressiblen Strömungen. Im Zusammenhang mit wichtigen Fragen, die bei äronautischen 
Anwendungen auftreten, hat sich in neuerer Zeit ein starkes Interesse am Studium zweidimensio- 
naler Strömungen kompressibler Medien entwickelt. Die Erfolge der Theorie im imkompres- 
siblen Fall legen es nahe, eine Verallgemeinerung der oben erwähnten funktionentheoretischen 
- Methoden für den kompressiblen Fall zu versuchen. Dabei sind gewisse Modifikationen erforder- 
“ ieh: Im inkompressiblen Fall können wir sowohl in der physikalischen x,y-Ebene (Strom- 
funktion Y als Funktion der Ortskoordinaten x, y) als auch in der Hodographenebene (Strom- 
_ funktion Y als Funktion des Geschwindigkeitsvektors; „Hodographenmethode‘) arbeiten, 
_  denn.sowohl 7(z, y) als auch Y(A4, ©) genügen der Laplace-Gleichung (v, © — Polarkoordi- 
naten des Geschwindigkeitsvektors, A@ = Igv).‘ | | 
® Beim Übergang zur Hodographenebene ändern sich viele Beziehungen’ und gewöhnlich 
_ ist es im inkompressiblen Fall bequemer, in der physikalischen Ebene zu ‚arbeiten. Im kom- 
-  pressiblen Fall dagegen genügt Y lediglich in der Hodographenebene einer linearen Differential- 
Nr gleichung, während sich in der &,y-Ebene eine komplizierte nichtlineare Differentialgleichung 
ergibt. Da die Theorie der nichtlinearen Differentialgleichungen nicht genügend entwickelt ist, 
erweist sich die Hodographenmethode bei kompressiblen Medien in vielen Fällen als bequemer, 
- obwohl beim Übergang zur Hodographenebene viele Beziehungen erheblich verwickelter werden 
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und das Arbeiten in der Hodographenebene bzw. der „logarithmischen Ebene“ verschiedene 
zusätzliche Schwierigkeiten mit sich bringt. 

Die Stromfunktionen sind in der Hodographenebene im allgemeinen auf Riemannschen 
Flächen definiert und besitzen verschiedene Singularitäten. So sind beispielsweise im inkom- 
pressiblen Fall Stromfunktionen Y(A*%, ©), welche in der physikalischen Ebene Strömungen um 
geschlossene Kurven liefern, im allgemeinen mehrwertig und besitzen Verzweigungspunkte, Pole 
und logarithmische. Singularitäten. Z.B. stellt 


vl re eA—U | Za—j)k —_-i9.. . (LI) 


ez* e* _ RU 


in der physikalischen Ebene eine Strömung mit der Grundgeschwindigkeit U 'um eine Ellipse 
mit dem Dickenverhältnis R? dar [11]. Vgl. Bild 1, 2 und 3. Die w, v-Ebene und die A4,0-Ebene 
werden zweimal überdeckt. Der Punkt P als Bildpunkt des Punktes 2= » ist Verzweigungs- 
punkt. Die einander entsprechenden Stromlinien sind in den Abbildungen entsprechend mar- 
kiert. Die Funktion 


(—n)e® | (hm)e® _sT, (ton 
| tn en rel) a 


wobei t=t(Z&) eine Lösung von 
Be ee Ul+ne* ıI' 1 | (1.3) en 
e—1l1+n (t—n)® 2raflt — nr en | 


ist und U, a, n, T' Konstanten sind, stellt das komplexe Potential ® + <Y einer Strömung um ein 
Joukowski-Profil dar [6, S.55—57]. Um daher für den kompressiblen Fall Stromfunktionen 
Y(}*, ©) mit ähnlichen Eigenschaften wie im 
inkompressiblen Fall zu erhalten, müssen wir 
mehrwertige Lösungen der Stromfunktions- 
differentialgleichung (2.3) bzw. (3.2) einführen 
mit ähnlichen Singularitäten wie im inkom- 
pressiblen Fall. 


erh — 


Infolgedessen entsteht die Aufgabe, für . 
die Lösungen {u} linearer partieller Differen- 
tialgleichungen eine der Theorie der harmoni- 
schen Funktionen analoge Theorie zu ent- 
wickeln. Die von dem Verfasser in diesem 

. Sinne unternommenen Versuche bezogen sich 


Bild2. In der Hodographenebene Bild3. Bild in der logarithmischen Ebene 


E 


hauptsächlich auf Verallgemeinerungen von Methoden der Funktionentheorie nach folgenden - 
beiden Richtungen: 
1. Methoden, die sich auf den Riemannschen Abbildungssatz gründen oder in enger Be- 
ziehung zu ihm stehen, ey 
2. Methoden, die sich auf die Tatsache gründen, daß analytische Funktionen eine Algebra 
bilden und viele nützliche Eigenschaften besitzen, welche uns gestatten viele Beziehungen zu 
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entwickeln, sowie auf die Tatsache, daß die Operation Re (Nehmen des Realteils) viele Eigen- 
schaften der analytischen Funktionen einer komplexen Veränderlichen erhält. 
Bei 1. ergeben sich die gewünschten Verallgemeinerungen durch Betrachtung der Klasse Ber 
= L?(B) von Funktionen, die ineinem Gebiet B analytisch sind, durch Einführung des Begriffs der zz 
E\ Kernfunktion K5(2,8) und mit dem Nachweis, daß man auf Grund der Eigenschaften der Kern- 4 
e funktion die Resultate erhalten kann, die in der klassischen Theorie aus dem Riemannschen 
 Abbildungssatz hergeleitet werden. In vielen Fällen können wir die funktionentheoretischen 
Sätze in dieser neuen Form unmittelbar verallgemeinern auf den Fall linearer partieller Diffe- 
rentialgleichungen vom elliptischen Typus, vorausgesetzt, daß die Koeffizienten ausreichende 
a Suesiohahten besitzen?). Diese Methoden sollen jedoch hier nicht weiter diskutiert 
werden. 
== 7 Wir wenden uns nun zu den unter 2. genannten Methoden. In der Funktionentheorie “ 
3 werden tiefliegende und interessante Beziehungen hergeleitet, die sich auf solche Fragen be- En 
ziehen wie z. B. das Studium und die Klassifikation von Singularitäten, Beziehungen zwischen Ba 
- den Koeffizienten der Funktionselemente und dem Verhalten der Funktionen im großen, das 
Problem der Darstellung analytischer Funktionen, die auf Riemannschen Flächen definiert 
sind. Mit Hilfe des Operators ‚Re‘ können diese Ergebnisse unmittelbar auf harmonische 7 A 
Funktionen übertragen werden. Das Ziel der unter 2. genannten Methoden ist das Studium Be 
2 solcher Operatoren, welche analytische Funktionen in Lösungen gegebener linearer Differentia- “u 
 —_, gleichungen transformieren und die erwähnten oder andere Eigenschaften erhalten. Dabeihaben — 
J e; wir zwei Fälle zu unterscheiden: 1. Differentialgleichungen, deren Koeffizienten in dem Bereich, z zu 
in denen die Lösungen betrachtet werden, analytisch sind; 2. Differentialgleichungen, bei denen en: 
5 ‚die Koeffizienten singulär werden. Bei den Differentialgleichungen der kompressiblen Ströo- 
- mungen handelt es sich um den zweiten Fall (vgl. $3). x a. 
a - Oft ist der eine oder andere Operator P, besonders geeignet um zu zeigen, daßeine gewisse 
5 Eigenschaft einer analytischen Funktion, auf die ein Operator P, angewandt wird, erhalten 
bleibt. Aus diesem Grunde ist es von Interesse, verschiedene Typen von Operatoren aufzustellen 
und zu untersuchen. 5 > Be; 
E. ‘ Falls die Koeffizienten a, b, c der Differentialgleichung  . 


Be x L(u) = Au+ au, +bw,+ cu=0 EEE a) 


‚ganze Funktionen von x, y sind, liefert die Operatorenmethode Funktionsklassen, die nur von 

‚gewissen Eigenschaften dieser, Koeffizienten abhängen oder überhaupt von ihnen unabhängig 
sind (vgl. z.B. 83). B | ee a i I 
Während bei regulären Koeffizienten a, b, ce der Integraloperator erster Art für die meisten 
Zwecke sich als der geeignetste erweist, scheinen bei singulären Koeffizienten andere (als 
_  Integraloperatoren zweiter Art bezeichnete) Operatoren von großem Nutzen zu sein, da sie 
im Gegensatz zu den Integraloperatoren der ersten Art in Punkten, in denen die Koeffizienten 
der Differentialgleichung singulär sind, definiert werden. er . >. 
Mit Hilfe der Integraloperatoren lassen sich verschiedene Kapitel der Funktionentheorie 
mit einer komplexen Veränderlichen nicht nur als Sätze über harmonische Funktionen sondern 
allgemeiner als Sätze über Lösungen der Differentialgleichungen (1.4) deuten, wobei ab, c 
entweder ganze Funktionen sind oder Sigularitäten der in dieser Arbeit betrachteten Art besitzen. 
Insbesondere sei hervorgehoben, daß verschiedene Ergebnisse, wie unsere Untersuchung 
zeigen wird, für eine ganze Klasse von Differentialgleichungen gelten. Wie das Verhalten har- 
monischer Funktionen typisch ist für das Verhalten der Lösungen der Gl. (1.4) mit ganzen m 
zienten so existiert im Falle singulärer Koeffizienten, wie sie in dieser Abhandlung be- 
t werden, ein Sonderfall, dessen Theorie besonders einfach ist. Wir werden diesen Fall 
ten Fall‘‘ bezeichnen. Das Verhalten der Lösungen des vereinfachten Falles ist 
ang typisch für das Verhalten der Lösungen von Differentialgleichungen einer 

lasse, nämlich mit Singularitäten, wie wir sie in $ 6 beschreiben werden. So gestattet 
n, den allgemeinen Charakter gewisser Gesetze für die Lösungen der Differential- 
deren Koeffizienten entweder ganze Funktionen sind oder’gewisse Singu- 
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analytischen Funktion und umgekehrt. Der Hauptzweck besteht jedoch nicht nur darin, ein ° 
Mittel zur Erzeugung von Lösungen zu haben, sondern es sollen Operatoren gefunden werden, 
die verschiedene Eigenschaften des Operators „Re‘“ behalten. 

Anmerkung: Es sei hervorgehoben, daß die Operatoren definiert werden für einen ge- 
wissen Bezugspunkt, den wir gewöhnlich in den Nullpunkt (Z=0) legen. Bei Änderung des 
gewählten Bezugspunkts Z ändert sich auch die einer festen analytischen Funktion zugeordnete 
Lösung, was beim Operator ‚Ze‘ nicht der Fall ist. 


8 2. Der Operator R, 


Wir betrachten eine isentropische Strömung mit der Druck-Dichte-Beziehung p=o0* 
(p = Druck. o= Dichte; o = const und k=const). Es empfiehlt sich, die Stromfunktion 9 
in der H, ©-Ebene zu betrachten, wobei 


v 


n=[|ı _k1 Dis a re: 2.1) 


Baht, v 
und v die allgemeine, v, die „‚kritischef‘, d.h. der Mach-Zahl M=1 entsprechende Strömungs- 
geschwindigkeit ist. a, ist die Schallgeschwindigkeit an einem Staupunkt. Dann nehmen 
die Gleichungen für das Potential ®(H,©) und die Stromfunktion Y(H,©) die Form an: 


ds = I Or = — (AH) Y, each ar er (2.1), 
B) « 
wobei ® = = usw. gesetzt ist. Durch Elimination von ® aus (2.2) erhält man für Y die 
Differentialgleichung: 
a 1—M? 
S(P) = U(H) a m 0, I(H) = N (23); 


(vgl. [4, (42)]). M ist die örtliche Mach-Zahl. 

In [2, S. 23—27] und [3, $2] hat der Verfasser den durch (2.6) definierten Operator R,(f) 
eingeführt, um Potentialfunktionen ®@(M,®) und Stromfunktionen Y(H,®) für die Gl. (2.2) 
zu erzeugen. Wir führen zunächst die Ausdrücke ein: 


®,1+i,1= (0 + iH)m 


- jo —2| HE | L(H,) u | an rt... | 


3 ! (2.4), 

Ei Fe H H Hs Hy 

2 +ul)erm [am 31) ES an, | U(H,)aH, | dH, + e 

” ß 0 

K On, 2 Pu? =ıi0 (@ 7 cd)" ’ 

H H Hi H, ri Ma 

RB ul) uenam.-s()e— [raman| an, | UHyaH, +... € 

f) 0 0 (2.5). 


Hl (er: EZ fi (HyaHR, 


Im Spezialfall H)=1 sind ® und Y durch die Cauchy-Ri i 

i y-Riemannschen Gleichun 
verknüpft und können als Potential und Stromfunktion einer inkompressiblen Strömung ne 
deutet um = a I(H) sind ® und Y durch die Gl. (2.2) verknüpft und können 
wenn wir =(l— 0? setzen, als Potential und Stromfunkti i iblen 
Strömung gedeutet werden (vgl. [3]). ee ee 

Wir wollen wie oben bereits geschehen, mit (O-+sH)'" und io(®-+iH)®} die A 

; | sch s { usdrücke 

bezeichnen, die aus (9+ vH )" durch die in den G1.(2.4) und (2.5) angegebenen anne entstehen 
Allgemein definieren wir die Operation R, durch die Beziehung: 


r h 
R, z+ iB,) O+iHA)| = m, (+:HA)" +ioß, (O+:H)W] ... . (2.6), . 


wobei %, P„ reelle Konstante sind. R, kann offenbar auf endliche S J 
Ä ummen an 
in manchen Fällen aber auch auf unendliche Summen, und erzeugt Potential uns ge 
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funktion einer kompressiblen Strömung. Der Operator R, hat ein für die späteren Betrachtungen 
wichtige-Eigenschaft: Die Funktionen y, mit k=1,2 mögen Entwicklungen der Form: 


= = Ion ern. 21020, 0, >. nn: (2.7), 

besitzen. Dann konvergiert die Funktion 
Y(H,0) -3 209, 9,1 (H,0) + 3 Pi (Ma) wen (2.8) 
in in Gebiet |0|+0|H]|<©,, wobei o hinreichend groß ist und Y(H,6) die Eigenschaft hat 
200. ,0, FA ee... 29). 


Bers und Gelbart [17] haben unabhängig vom Verfasser denselben Operator gefunden und 


in weiteren Veröffentlichungen [19, 20] untersucht. Sie gebrauchen für die erhaltenen Funktionen 
die Bezeichnung ‚‚3-monogen‘. 


$ 3. Die Differentialgleichungen der Gasdynamik in der pseudologarithmischen Ebene 

Bei Überschallgeschwindigkeit wird die Differentialgleichung (2.3) hyperbolisch; das Ver- 
halten der Lösungen dieser Gleichung ist völlig verschieden vom Verhalten der Lösungen von 
Differentialgleichungen des elliptischen Typus. Wenn wir daher eine weitgehende Analogie zum 
inkompressiblen, von der Laplaceschen Gleichung beherrschten Fall zu erhalten wünschen 
müssen wir uns auf den Unterschallbereich beschränken und haben daher zunächst die Strom- 
funktionen von Unterschallströmungen in der Hodographenebene zu betrachten. Im Unter- 
schallbereich können wir Gl]. (2.3) auf die kanonische Form der elliptischen Differentialglei- 
chungen durch Einführung einer neuen Variablen an Stelle von H bringen. Mit 


1 | 1 R nn at 
——] = = Pe Ai Fe I nenn —- a 5 
eG 2h°|1—hyı- m: Yrri EN 
(vgl. [4, $7]) geht Gl. (2.3) über in 
rear BE Re a Ar 
N aaa 1 Se 


2 
Anmerkung: Wir werden mehrfach Y*— — an Stelle von Y benützen, wobei?) 


1 
k—1 # SED 
— (1 — M23\14 a) N 

H=(1—M) (1 aM 
—= 8, (— 24-711 +8, (— 223 +8, (— 2A) +...) 
gesetzt ist [9, (4.3)]. Y* genügt einer gegenüber (3.2) etwas einfacheren Gleichung, nämlich 
A ER de 
HR 99 


(3.3) 


N EEE . (8.4) 


wobei E 
(k +1) M[—- (8k —1) Mt—4(3— 2%) M? + 16] 
42 16 213 ve 
(1— M) 
BE 5 (3.5), 
ns 78 BA AN ee]; Bor rene M=0... 
ist [4, (70), (71), 9]. 
Für die Potentfalfunktion ® erhält man die Gleichung: 
28 9® od 
ne ne a ee ll 3.6). 
- IE ATi Go 


Anmerkung: Es sei hinzugefügt, daß esin manchen Beispielen vorzuziehen ist auch im 


„gemischten Fall“ (d.h. wenn es sich teilweise um Unter-und teilweise um Überschallgeschwindig- 


23) Durch ein Versehen würde eine Koordinate, vgl. (2.1) und eine Funktion, vgl. (3.3) (die nichts 


miteinander zu tun haben) mit demselben Buchstaben 4 bezeichnet. Dies soll zu keiner Verwechslung 
führen. A 


\ 
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keiten handelt) die Veränderlichen A, © zu verwenden. Dafür M>1 die in Gl. (3.1) eingeführte 
Funktion A in 


—2M)=—i|aretg VRRZT)— + arcetg (HZ) = iA) EURE) 


übergeht, ist es nützlich, die in Frage stehende Funktion zu komplexen Werten des Arguments A 
zu verallgemeinern. Wir setzen A,+?A an Stelle von A und lassen dann den unteren Index 1 
in A, wieder weg. Betrachtet man dann Y(A-HiA,0), so ist A=0 im Unterschallbereich und 
A=0 im Überschallbereich. Im ersten Fall genügt Y(A,@) der Gl. (3.2), im zweiten Fall aber 
der Gleichung 
le A u oY k+1 M“ 
ea 4N, Ey a 0, N Be (3.8). 


Die Halbebene (A=0, A<0, — o<@<oeo) entspricht dem Unterschallbereich, die Halb- 
ebene (A=0, A>0, — © <:9 < oo) dem Überschallbereich. Beide Halbebenen hängen längs 


der Achse A=0, A=0 miteinander zusammen. = 


$ 4. Integraloperatoren zur Erzeugung von Stromfunktionen für Unterschallströmungen 


Die neuerdings in die Theorie der partiellen Differentialgleichungen eingeführten Integral- 
operatoren [1] stellen ein nützliches Werkzeug zur Erzeugung und Untersuchung von Strom- 
funktionen für Unterschallströmungen dar. Ein Integraloperator ist ein Operator, der eine 
Klasse von Funktionen, etwa W, in die Lösungen einer gegebenen partieHen Differentialgleichung 
transformiert. In unserem Falle ist dies die Gleichung (3.2) und die Funktionenklasse X ist bei 
Unterschallströmungen die Klasse der analytischen Funktionen einer komplexen Veränderlichen 


Z=A-+i0 oder Z=A—iO. Der Integraloperator kann entweder in der Form 


t=1 
Fr zZ’ di 
P(f; Z,) = E(Z,Z,ı)-f e (1 u) a NE EL in (4.}) 
mit Z=2Z2—-Z, Z=Z-Z, 
geschrieben werden oder in der Form 
_. ZZ I 
PZ)=H 112) ION fen | gmaz..-. Z 
/ 5 (4.2). 


zZ 
= N nz’, Z : 
{) 


Dabei sind f und g willkürliche analytische Funktionen einer komplexen Veränderlichen und 
E(Z,Z, t), H und die Q® sind gewisse feste Funktionen. Die Funktionen / und g werden als 
„P-zugeordnete“ bzw. „p-zugeordnete“ bezüglich der durch die Integraloperatoren erzeugten 
Lösungen der Differentialgleichung (3.2) bezeichnet. f und g sind dabei durch die Beziehung 


n/2 
(azyır aurg(az) 1 2 [,... „dg(2Zsin®# 
SO—TUR san zT, J nn Eanrar: 


miteinander verknüpft. Für manche Zwecke erweisen sich Operatoren der Form (4.1), für 
manche Zwecke solche der Form (4.2) als besser angepaßt. Der Ausdruck 


R B=A|ı ar Samzeer] A BE NET ER ya (4.3) 


n=1 
£ 


wird ‚als erzeugende Funktion des Operators P bezeichnet. Z, ist der Bez spunkt“ 
. .. * ; n: des 
Wir setzen gewöhnlich Z,=0'und schreiben P(f) und 7(9) an Stelle es KK f: 0) und 


- 


Rz: ‚angew. alt th.M M 


Bd. ‚32 Nr. 2/31 Febr, a 1952 eB ergman, Operatorenmethoden in der Gasdynamik 39 


' Zwei Re le raferen wurden in der Theorie der kompressiblen Strömungen we 


E- er ee nämlich die sog. Integraloperatoren erster Art und zweiter Art?) P, bzw. Psivgb 
Se alt: 2, = und [3—7, I). Die Ber Funktion g(Z) mit dem Realteil 
. PFZD)=-ıKay)=z 2) +3). Z=atiy Z=a-iy.. (457 


E ‚kann erhalten werden, wenn man die Veränderlichen «, y zu komplexen Werten fortsetzt und 
© ‘sie dann in der charakteristischen Ebene z=iy betrachtet. 


Man kann zeigen, daß für jede Differentialgleichung LP = a ze ir Pr— (0, work — | 


‚eine Bene Funktion ist, ein Integraloperator existiert (d. h. die Q(m) re werden können) 


‘derart, daß die der Lösung Y* (Z, Z) entsprechende assoziierte Funktion 9(Z) durch die näm- 
au Beziehung 

gZ)=2W*(Z, 0) + const ee Ar er . (4.6) 
Sachen wird wie im Fall der harmonischen Gleichungen. Dieser Operator ist ein Tier 
operator erster Art. (Der Integraloperator zweiter Art wird nur in dem Fall, indem nur von 
abhängt, definiert.) 


Anmerkung: Der se erster Art kann auch in der Form. 


By a \ 5 er 22 2,23 h = 
eg ZH et [roaz,az, | Hs [Fgaz, 12.32.42, — ER | an BR; 
x j. iR Z Zo Ze Z Zo Ei MT & 
_ geschrieben werden. Re. \ Be 


Mit Hilfe des En ealoneaters erster Art können wir verschiedene Ergebnisse der un 
tionentheorie als Sätze über Lösungen partieller Differentialgleichungen deuten. Diese Resultate 
"haben zwar keine ünmittelbaren Anwendungen in der Gasdynamik, bieten aber eine interessante 
_ Illustration von Ergebnissen, die man mit Hilfe der Operatorenmethode erhalten kann. 

: 1. Die reelle Lösung Y= Re[pı(g)] ist regulär in jedem einfach zusammenhängenden 
Gebiet, in dem g regulär ist, und umgekehrt [23, S. 135, 136]. en 
2. Wenn gin Z=« einen Verzweigungspunkt endlicher Ordnung hat, darn hat rı(g)in a 

Ks einen Verzweigungspunkt derselben Ordnung. Wenn g in & einen Pol hat, dann wird p,(g) in & 

ı derselben Ordnung unendlich, hat-aber einen Verzweigungspunkt unendlich hoher Drdmnnez 

Vir bezeichnen « als „‚polähnliche Singularität‘‘ [23, S. 318, 319] der reellen Lösung. 

83. Für jede Differentialgleichung ‚L(d) = 0, vgl. (1.4), gibt es eine Folge komplexer Br 

Lösungen‘ d,(2,y) = U, (2, Y)+VV, My), v = 02% 2,2.. „von denen jedeeine a Funktion 
von zund y yist, die als Analogon von (+ iy)? betrachtet werden kann. i En 

a) Die ®, verhalten sich ähnlich wie die Potenzen für große Werte von v. Genauer eilt, ; 
-|8,(2,y) —H (, 3) (+ iyP] s(C,(®,Y) Gas Wat wobei H und CO, ganze, von vun 

ge Funktionen sind. — 

Jede ee w, die in einem Kreis + P<o: regulär ist, kann dort in eine > gleich- 


ig ne] Zus) 


= Y von Lı ‚) zuwı oY=0, die in einem einfach zusammenhängenden Beh. 
re h regulär ist, ‚kann dort dürch eine Linearkombination. Re 3 I: a) o,| 5 2% 


v-1 "AO 
on des. Run [3 eschen Satzes). Hier bedeutet C' eine ganze 


ngen (die wieder Beziehungen in der Fin ktlonentheois 
n einer Lösung 7 mit der Reihenentwicklung 
im punkt: und dem nalen der Folge 


25 


Be 


’ 


ei eher k- ter Art(k=1 2). In den Arbeiten Br Bi - 
die dort mit P und p bezeichnet werden) Integraloperatoren Bi. 


onen einer komplexen Veränderlichen in komplexe : 
or (d.h. die erzeugende Funktion) andern 


Igende. wichtige Bigenschaft: ‚Wenn wir die 
mente A = @ en d. h. wenn wir zZ 2 
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a) N A„„2”2" ist regulär in jedem’ Gebiet, in dem 5) Amo 2” regulär ist. Man beachte, 
un na 
daß infolgedessen die Lage von Singularitäten von Y durch die Ano, m=0, 1,:2, . .. bestimmt 
ist, unabhängig vom Koeffizienten Oder Gleichung L,(?)=0 [23, 8.140]. 


b) Wenn 3 |Amo|? < oo, dann ist Y= 3 A„„2”z" im Einheitskreis #+y?< 1 regulär 

m=0 . . . . 
und besitzt ‚‚fast überall‘ Grenzwerte für radiale Annäherung an den Einheitskreis. Wieder 
sind die Punkte, wo die Grenzwerte existieren, unabhängig vom Koeffizienten € [23, S. 141]. 
5. Man kann unmittelbar verschiedene Sätze der Funktionentheorie einer komplexen Ver- 
änderlichen übertragen auf den Fall der Lösungen Y der Ditferentialgleichung LAP)=B-88 


liefern beispielsweise die Hadamardschen Bedingungen, daß eine Funktion  Ano 2" in einem 
Gebiet meromorph ist, die notwendigen und hinreichenden Bedingungen in den Gliedern der 
Unterfolge {A„0} dafür, daß Y= 53 3 Ann?”z" in einem Gebiet B nur polähnliche Singu- 
laritäten hat. en 


6. Mit Hilfe verschiedener Sätze über die Riemannsche Fläche von Funktionen einer 
komplexen Veränderlichen kann man die Theorie der Riemannschen Flächen für die Lösungen ’ 


von Z,(?)=0entwickeln. Wenn Y(z, z) in der charakteristischen Ebene z=0 eine algebraische . 


Funktion wird, ist die entsprechende Riemannsche Fläche von Y eng verwandt mit der 

Riemannschen Fläche der Funktion g, die mit 7 durch die Beziehung (4.6) verknüpft ist. 
Insbesondere kann man mit Hilfe der Theorie der Integraloperatoren der ersten Art die 

Theorie der mehrwertigen Lösungen linearer partieller Differentialgleichungen entwickeln. 


$ 5. Der Integraloperator zweiter Art 


Bei regulären Koeffizienten der Differentialgleichung ist der Integraloperator erster Art 
für die meisten Zwecke am geeignetsten. Wenn aber der Koeffizient N in Gl. (3.2) singulär ist 
und wir das Verhalten der Lösung für A= — oo oder auf der singulären Linie A=0 untersuchen 
wollen, ist der Integraloperator zweiter Art von beträchtlichem Nutzen. Im Gegensatz zum Fall 
der Integraloperatoren erster Art, wo die Koeffizienten Q® Funktionen von zwei Variablen A 
und © sind, hängt in diesem Fall Q(® nur von A ab, so daß die Tabulierung einfacher und der 
Operator für numerische Rechnungen besonders geeignet ist. Die Operatoren zweiter Art sind 


definiert für den Fall, in dem der Koeffizient N nur von/= 4 (Z-+Z) abhängt; sie konyergieren 


8 nur in Unterbereichen des Regularitätsbereichs 
von N. Sie haben jedoch den Vorteil, daß einige 
von ihnen auch noch auf der singulären Linie 
A=0 von N und für A= —oo definiert sind. 
Wenn man A ins Komplexe fortsetzt, d.h. A 

A durch A+vA ersetzt (vgl. Schluß des $3), dann 

istim Unterschallgebiet A= 0 und die Zugeord- 
nete f(A+:0©) ist eine Funktion einer kom- 
plexen Veränderlichen, während im Überschall- 
gebiet A=0 und f(i(A+0))=s(A-+ 0) 
eine Funktion einer reellen Veränderlichen- ist 
(vgl. Bild 4). Auf diese Weise wird auch für die 
Untersuchung des gemischten Falles der Opera- 
„ tor zweiter Art herangezogen. (9,11,12). In 
1 
dem Spezialfall, indem N=-—— — __ gilt (wir 
"Bild 4 5 6(Z+Z) a 
werden ihn als den ‚vereinfachten Fall‘ oder 
als Tricomi-Gleichung bezeichnen), ist die erzeugende Funktion Z= E} des Operators zweiter 
5 : 

Art die hypergeometrische Funktion von u= In diesem Fall wird die Funk- 

tion UM) in (2.2) und (2.3) (HA)= — CH, A* und H sind verknüpft durch die Beziehung 

—H=0(— At)? und es ist 4 — m [9,$. 862]. Die Theorie der hypergeometrischen 


Funktionen ist genügend ausgebaut und wir können verschiedene Eigenschaften der hypergeome- 
trischen Funktionen benützen um eine Theorie der Lösungen der Gl. (3.2) im vereinfachten Fall 
zu entwickeln. Wenn man den in diesem Paragraph beschriebenen Operator benützt, lassen sich 


* 


Unterschallgebiet 


EN ADS 
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verschiedene Sätze ähnlichen Charakters, wie sie im Paragraph 4 angegeben sind. für die Lö 
der „vereinfachten“ GI.(2.3) aufstellen [11, 12]. Der in en Pe ehe 
derfall gibt uns ein Muster, wie man im allgemeinen Fall vorzugehen hat; die Theorie der Inte- 
graloperatoren zweiter Art ist eine Verallgemeinerung bekannter Sätze über hypergeome- 
trische Reihen auf den allgemeinen Fall, in dem E nicht auf eine hypergeometrische Funktion 
emer geeignet gewählten Kombination von Veränderlichen zurückgeführt werden kann. 

Insbesondere gestattet uns die Theorie der hypergeometrischen Differentialgleichung 
die erzeugende Funktion in der Form von zwei unendlichen Reihen darzustellen, die für 
Ar +iO)(— 2A] <1 bzw. |(A++iO)/(— 24+)]>1 konvergieren. In $6 werden wir analoge 
Reihen im allgemeinen Fall herleiten. Die erzeugende Funktion des Integraloperators zweiter 
Art ist in diesem Fall: 


E*+=4,F(, : 


5 RABEN: a ©) 
ezetahggg Sr 


wo A, und B, Konstanten sind. Fist die hypergeometrische Funktion. (5.1) kann dargestellt. 
werden. durch: 


k 


[) n— — + — 
ER DI (eZ+}) ? °Qrmm(22+), (k=1,2; Z=4-+10). . (52) 
für |Z+/2 A| <1 und 
1 2k 
} n— — + — 
Bes (k) ai le (2 Ar) — t2Z+) 2 3 (k =, 2) EN N (5.3) 


für |[Z+/2At]> 1, wobei 1 


Ei Br 


een (64) 
Qg+Hom(2rr—=1 
und a ) ( 3 » er 
ee el en a on 
| Ben 
| x 
ie Er Fee) 
A 5 75 5 | 
| ee 


Q+m2 (2%*) brauchte nicht mit angeführt zu werden, da der entsprechende Operator nach 

Integration identisch Null liefert [9, 84, 11]. 

“ In dem vereinfachten Fall befriedigt die erzeugende Funktion E*+ zweiter Art, betrachtet 
als eine Funktion von %, die’'hypergeometrische Gleichung, d.h eine gewöhnliche Differential- 
gleichung von Fuchsschem Typus. Es war nicht möglich in dem allgemeinen Fall die ent- 
sprechende Gleichung für die Erzeugende auf eine gewöhnliche Gleichung zu reduzieren. Andrer- 
seits kann man zeigen, daß im allgemeinen Fall die erzeugende Funktion viele gleiche Eigen- 
schaften besitzt, wie die Erzeugende der vereinfachten Gleichung?) [12]. 24 


$ 6. Der Integraloperator zweiter Art für die allgemeine Gleichung 
Wir wollen nun zeigen, daß wir in gewissem Umfang dieselbe Lage bei den Gleichungen 
des gemischten Typus haben, die in der Theofie der kompressiblen Strömungen auftreten, wo 
der vereinfachte Fall als ein Muster für die allgemeine Theorie betrachtet werden kann. Es wird 
angenommen, daß N in (3.2) eine Entwicklung der Form 


x var San) ee (6.1) 


34) Diese Betrachtungen verallgemeinern gewisser Kapitel der Theorie der Gleichungen von 


-  Fuchsschem Typus auf den Fall von partiellen Differentialgleichungen. 


Dr 
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hat, die für —A,<A<0 gilt. Wenn r, der Konvergenzradius der Reihe 2 ß,x’ ist, dann 
y= 
ergibt sich A, = r}?. 
Wir stellen jetzt folgende Forderungen an die Funktion N(A): Wir sagen, daß die 
Funktion N(}) der reellen Veränderlichen A zur Klasse N gehört, wenn 


N(A) eine analytische Funktion für — o<A<0O und reell für A<0 ist (6.22); 


N(2) in der Umgebung von A=0 eine Entwicklung der Form (6.1) besitzt, \ (6.2b); 
die für —A,<A<0O giltmit 0O<A,so 


A 
das Integral exp - [2N@d ı =H(4) für alle A existiert - - -.... (6.2c). 


Unter diesen Annahmen läßt sich zeigen, daß lim {(—A)U®H (A)} existiert. Wir bezeichnen 
10 


diesen Grenzwert mit S, und fordern außerdem noch 


8,307: 2 AERh Re Re Se + (6.2d). 
Es seien- ö 
[.) RER PC h 
EP=yN(ez) ’ ’Qmm(2}), (k=1,2; A<0, |0%<3j4%) . 22 (6.5) 
n=0 
und > ER „2% 
EP — 5 (—%Z) 3 gm) (22), k=1.2% |9®>3AN).. ..... (6.4), 
n=0 


wo die Q%,% durch die Rekursionsformeln 


an HH" +argmn—0, QM) —=0, | 
1—k 


1 — (6.5) 
Qu — —AF, F%)=— N(A)— N:0), QM — (2A) | 
und die g%# durch die Rekursionsformeln 
an +4F0)g"m—0, 
® (+3 e) ak + a E= 4 F(A)gr tim = 0 (k = ik 2) ee (6 6) 
verknüpft sind und die in (6.8a) bis (6.8d) gestellten Anfangsbedingungen erfüllen. 
Für die q"-®(}) erhalten wir Reihenentwicklungen 
1 2 
o& R=—r nr) 
gm ®) (2) = > AN er Aa See a (6.7) 
mit 
a 
ed = 218, Non) 2 ©5755), (Oele 
er ee ER ET RER AL (6.8b), 
np > 
EN GH DE a (6.8c), 
ONE A Es RER PT Nr (6.8d). 
Dann erhalten wir für die erzeugende Funktion zweiter Art 
EP,=HE® = 1, 2) ee (6.9) 


(6.3) konvergiert für |0®<3]A]2, (6.4) für |O|?> 3]2]2. (6.3) und (6.4) stellen die Verallgemei- 
nerung der im vereinfachten Fall erhaltenen Entwicklungen (5.2) und (5.3) auf den allgemei- 
nen Fall dar [9, 85; 12]. 


$ 7. Der modifizierte Integraloperator 


Bei der Verwendung der Operatoren p, und 9, ,*) zur Erzeugung von Strömungen kompres- 
sibler Medien sind gewisse Modifikationen nötig: 

Die Integraloperatoren erster und zweiter Art erzeugen Stromfunktionen, die in den 
meisten Fällen die nämlichen allgemeinen Eigenschaften besitzen wie die Funktionen, auf die 


Vgl ne (4.2); ?,,ı wird ebenso wie p, gebildet, jedoch sind die Q)(Z; Z) durch QW)(A) zu ersetzen. 
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sie angewandt werden. Wenn beispielsweise die Funktion g einen Verzweigungspunkt n-ter 
Ordnung hät, wird die erzeugte Lösung der Stromfunktionsgleichung einen Verzweigungspunkt 
__von derselben Ordnung und im gleichen Punkt wie g haben. Es gibt jedoch einen wichtigen 
Auspahmefall, der eine Modifikation bei der Anwendung des Verfahrens in der Gasdynamik 
erfordert. Wenn wir nämlich den Operator p auf eine Funktion g anwenden, die einen Pol oder 
eine Singularität besitzt, dann führt der Operator eine wesentliche Modifikation des Verhaltens 
der Lösung in dem in Frage stehenden singulären Punkt herbei derart, daß sowohl der Realteil 
als auch der Imaginärteil von p[lg (Z— Z,)] und p[(Z—Z,)"], wobei n eine positive ganze 
Zahl ist, im Punkt Z=Z, unendlich-vielwertig werden. Da es beim Übergang zur physi- 
kalischen Ebene wichtig ist Singularitäten zu haben, für welche wenigstens eine der Funktionen ® 
oder Y im singulären Punkt einwertig ist, müssen wir unsern Operator in der folgenden Weise 
- modifizieren: > 
Wenn wir eine Funktion g mit einem Pol oder einer logarithmischen Singularität haben, 
zerlegen wir sie in zwei Teile. Der erste Teil g, soll weder Pole noch logarithmische Singularitäten 
haben, der zweite Teil sei 


7 
1 


= a \ SAZZEZFH N KEZ—Z):... 2.02... (2 
S a - i k=1 k=m+1l B 


wobei A, Konstante und v, ganze positive Zahlen sind. 
Mittels des Integraloperators erster Art p, erhalten wir eine Lösung von (3.2), etw. 


erento) en | 
0 Durch Einführung sog. Fundamentallösungen: = 
= EN Z) AZ EZ BZ N 0) 


von (3.2) oder (3.4), wo Aund B gewisse einwertige Funktionen sind, die man leicht bestimmen. 
kann, erhält man Lösungen von (3.2) oder (3.4), die dasselbe Verhalten zeigen wie lg |Z—Z,. 


Siehe [4, S. 45; 6, S.43; 7, 8]. Differenziert man diese Funktion nach dem Realteil und Imagi- 
“närteil des Parameters Z,, so bekommt man einwertige Funktionen 741) und 4.2) mit einer 

Unendlichkeitsstelle erster Ordnung. Wir bezeichnen mit Y12.2 (Z; Z,) die Stromfunktion, die 
der Grundlösung 842) (Z; Z,)= ((Z; Z,)1g|Z2— Z,|+D(Z;Z,) der Potentialgleichung ent- 
- spricht. Vgl. (5.2b), (2.8) und (2.15a) von [8]. - > | 


Anmerkung: Die durch Yun, ya, Y4D) und %4,2) erzeugten Strömungen haben 


ri 
r 
“ 


wie g. Wir werden diesen Prozeß zur Bestimmung von Stromfunktionen eines kompressiblen 


a 2 z Wr g(Z) =; 
[14] der 


nungsebene eine Strömung um einen nahezu elliptischen Zylinder. 


Strömung mit Schallg :schwindigkeitsdurchgang (für den verein- 


ae IN 


der Theorie der inkompressiblen 
weitgehend gleicht. Im folgenden 
orie der inkompressiblen Strömungen 


Es (A 


Operator Pz beispielsweise auf die Funktion (1.1)angewandt worden. 


licht uns die Entwicklung einer 
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I. Bedingungen dafür, daß eine in der 2, 0-Ebene definierte Stromfunktion in der 
Strömungsebene eine Strömung um eine geschlossene Kurve liefert 


d E Se 
Bei der inkompressiblen Strömung ( 77 Strömungsgeschwindigkeit = veie) hat man 


Z = [e-2*an = [est na ee (75) 
mitw=®-+iY, Z2=14—-i0, A* =Igv. Nimmt man dann die reellen und imaginären Teile 


und berücksichtigt die Cauchy-Riemannschen Gleichungen für ® und Y, so erhält man 
für x und y Formeln, die den im kompressiblen Fall geltenden Beziehungen 


u i in® 
2 1 eos 8.2, ° = .w.|a0), 
oll ® 


—_M2 
=[41-: BE Re 


w dv+ lcos 0:9, — 
(7.52) 


) 
27 dv-+ Isin 8:97.42 v.) 10) 


v2 v v 

entsprechen. i 
Wenn nun der Definitionsbereich B einer Stromfunktion Y (in der A, ©-Ebene) einfach 

zusammenhängend ist, dann ist dafür, daß die Randkurve des Bildes von B in der %, y-Ebene 

eine geschlossene Kurve ist, notwendig und hinreichend, daß die Summe der Residuen der 

Singularitäten [e* ES im Bereich B gleich Null ist. 


Dieses Ergebnis kann auf den Fall der kompressiblen Strömung verallgemeinert werden. 
Da die Integranden von (7.5a) vollständige Differentiale sind, hängen die Werte der Integrale 
vom Integrationsweg nicht ab, vorausgesetzt, daß kein singulärer Punkt überschritten wird. 
Wenn wir daher in (7.5a) die Fundamentallösung (7.2) bzw. ihre Ableitungen nach Re Z, oder 
Im Z, einsetzen, können wir jeder solchen Singularität zwei reelle Zahlen 


29 (Zy Fund. Terz): SE (7.6) 


zuordnen. (Im inkompressiblen Fall sind diese Zahlen der Realteil und Imaginärteil des Residu- 
ums in Z,.) Wenn der Bereich B mit der Randkurve b einfach zusammenhängend ist, hat man 


2 
3 4, Te (Z)=0, RE Br EEE (7.7) 
v=0 ? 


als notwendige und hinreichende Bedingung dafür, daß das Bild von 5b im der x, y-Ebene eine 
geschlossene Kurve ist. Dabei sind die A, die in. (7.3) auftretenden Koeffizienten [4, $ 14; 6; 
8, 86]. 


I. Verallgemeinerung der Blasiusschen Formeln 
Eine weitere Anwendung unserer Methode ist eine Verallgemeinerung der Blasiusschen 
Formeln für den Auftrieb F und das Moment M, die auf den umströmten Körper ausgeübt 
werden. Wir nehmen an, ein Profil B sei in eine kompressible Strömung eingetaucht und die 


Stromfunktion in der A, O-Ebene habe die Form (7.3). Dann erhalten wir für die auf den Körper B 
wirkende Kraft: 


iR 
F=F(N)— — j v(A)eiedP + io [ hi u a] VAR? (7.8) 
mit £ ü 
1% 1-M „dv 9 dv .W 
da=2ee | = Port = laa+lw/®+:2elae) FR N re 3 (7.9). 


o ist die in $2 (erster Absatz) eingeführte Konstante. 
C ist eine hinreichend glatte und kleine geschlossene Kurve um den Punkt Zu In Ver- 
allgemeinerung der Blasiusschen Formeln erhält man so im kompressiblen Fall: 


F=AF(Fn) LAFPAD AL AFPUD) . 2.2.0... (7.10). 
Dabei sind F(Y), ... Vektoren, die sich dadurch ergeben, daß man in (7.8) an Stelle 


von Y die Funktionen YD), ... einsetzt. Eine analoge Formel ergibt sich für das Moment M 


[8; 86]. 
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Ein Beitrag zur Statik der Kreisplatten 
Von H. Jung in Stuttgart-N 


Durch die Anwendung der Hankeltransformation ist es gelungen, Differentialgleichungen auf einem 
neuen einfachen Weg zu lösen. Das Verfahren wurde zur Bestimmung der Durchbiegung von Kreisplatten 
bei rotationssymmelrischer Lastverteilung herangezogen. Auf demselben Weg lassen sich Verfahren zur 
Behandlung der Druckverteilung unter Kreisplatten, die auf einem elastischen Medium aufliegen, finden. 
Ebenfalls kann auf diesem Weg die ‚‚Wärmebeulung von Kreisplatten“ behandelt werden. 


Applying the Hank el transformation, it is recently possible to solve dijfferential eguations in an easy 
way. The author makes use of this method to determine the deflexion of circular plates under rotationally- 
symmetrical strain. The same method is applicable to the determination of the pressure distribution under 
elastically bedded eircular plates. Tihe ”buckling by heat” of circular plates can be treated in ihe same way. 


En employani la transformation de Hankelon a r&cemment reussi & rösoudre des Equations diffErentielles 
d’une manitre nouvelle et simple. La methode est employde a la determination du courbement de plaques 
circulaires & une distribution de charge symmeirique & la rotation. Dela m&me maniere on peut trouver des 
methodes pour la manipulation de la distribution de la pression sous des plaques circulaires couch&es Elasti- 
quement. La ‚‚bosselure‘‘ de plaques circulaires causde par la chaleur peut Eire traitde de la m&me manitre. 


IIpnmeuemme Tpauchopmaunn Taukena a0 BO3MO3KHOCTB pemeHnuna AuhhepeHmuasnsHslX 
ypaBHeHHÜ HOBBIM IPOCTBIM IyreM. Ci1oco6 ÖsIı HpHuBIeyeH 1A OMpereneHHA Iporuda KPy- 
TOBBIX INIACTHH IPH BPAImaTeIBHO-CHMMEeTPHYecKOM pacnpererreHHu Harpysku. TeM :xe caMbIM 
IIyTeM MorTyT ÖbITb HaÜNeHBI CIHOCOÖBI TPAKTOBEH PAcIıpere.reHuA MaABIeHHA Meikay KPYTOBBIMH 
ILNIACTHHMH, TIIOKOAIIHMHCH Ha yupyrof cpene. ITHM ;Ke IyTeM MOKeT ÖBITb HccJIeNoBaAHO 
„TEILIOBOe MATbE KPYTOBbIX INIACTUH““. 


1. Einleitung 


Wird eine Kreisplatte durch eine rotationssymmetrische Belastung normal zu ihrer Mittel- 
fläche beansprucht, so erhält man die Durchbiegung f(r) aus der Differentialgleichung 


n 
aa EN 7 a a 
In (1) bedeuten J/() die Durchbiegung, 
»(r) _ die Belastung, 
m?Eh? et RE B 
N= m) die Biegesteifigkeit, 
ed 
A= gar Tan 


Wendet man auf (1) die Hankeltransformation (Ziffer 2) an, so läßt sich wie in Ziffer 3 gezeigt 
wird, ein neues Verfahren zur Berechnung der Durchbiegung aufbauen. Der Vorteil dieses Ver- 
fahrens gegenüber dem Verfahren von Pichler!) ist, daß in die Rechnung nur die Randbedin- 
gungen der Platte eingehen. In. Ziffer 4 werden einige wichtige Unterfunktionen für Belastungen 
hergeleitet und der Zusammenhang dieser Unterfunktionen mit den-in Ziffer 3 auftretenden 
Integrationskonstanten gegeben. Wie die auftretenden Integrale ausgewertet werden können 
wird an Hand von einigen Beispielen in Ziffer 5 gezeigt. In Ziffer 6 wird zuerst die rotations- 
symmetrische Knickung der Kreisplatte konstanter Biegesteifigkeit streng auf einem neuartigen 
Weg behandelt und dann ein Iterationsverfahren für Kreisplatten mit stückweis konstanter 
Biegesteifigkeit mit Hilfe des in Ziffer 3 angegebenen Verfahrens hergeleitet. . Ziffer 7 behandelt 
die dicke Kreisplatte mit rotationssymmetrischer Belastung. % 


2. Hankeltransformation 


Die Hankeltransformation wird auf dem von Weyrich?) gezeigten Weg aus der 
Fouriertransformation hergeleitet und dann erstmalig ihre Anwendung auf zwei einfache 
Differentialgleichungen gezeigt. 

Gegeben sei eine Funktion f(z, y), die im endlichen Gebiet stückweise stetig und zweimal 


‚ stückweise stetig differenzierbar vorausgesetzt wird. Ist die Funktion f(x, y) noch im ganzen 


Gebiet quadratisch integrierbar, so läßt sie sich durch ein Fourie rintegral darstellen. 
f(® Na ef /(s en 
» Am » dsdtdudv ...2.. 2) 


1) O.Pichler: Dissertation. Stuttgart 1927. Berlin 1928 
2) R. Weyrich: Die Zylinderfunktionen und ihre Anwendung. Leipzig 1937, 8.109. 
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In (2) führt man an Stelle der cartesischen Koordinaten x, y und s,t die Polarkoordinaten r,@ 
und 0,0 ein. In komplexer“Schreibweise wird damit 


(2) +ily)= Yr Fo —2rocos (Po) eir, 
wobei 
= =—=/N 
Ist: RS 
Mit 
u+iv=deiftw 
geht dann das Integral (2) über in 


+7” 
Fr, p9)= a [[reran;t a See dß da de da Run (3). 
f() f) nr —ı 
Die Integraldarstellung von Hansen 


Srel mern Bo) ger cos pr +e® —Broos(o=a) dß 
führt (3) über in 
F(r,9) = — ei PR oF(o,a) Ig(A,Yr®+ 0? —2rocos(p—o))dadedA . : . . (A). 
Das aniihkcen a: 
a0) SE | | #e > In (4,7) I1(4,0) f F(e,a) eire=@) dadedi. . . (8). 
Macht man nun für F(r, u). den speziellen Ansatz 23 N 


F(r,0)= © (r) e!mr 
und geht damit in (5) ein, so wird unter Berücksichtigung der Orthogonalität der Kreisfunktionen 


= [Ina 4 [®(0:07.0,0)de EEE Bi. 
0 0 z 
Schreibt man (6) in Form einer Transformation, so wird s 
ga) |): J.@0)de, 
3n-] HA) -A I, (Ar) dA 


Es soll nun gezeigt werden, wie die Transformation (7) zur Lösung von Differentialgleichungen 
herangezogen werden kann. 


Vorgegeben sei die Differentialgleichung 


Macht man für f(r) die Voraussetzungen 


1. f(r) stückweise stetig und. stetig differentierbar, 
2. f(r) beschränkt, 
OO in Osrxe 
und b<r<o 
so läßt sich f(r) lenmiällig darstellen. Mulebact man nun (8) mit r-J,(A,r) durch und inte- 


 griert partiell über r, so wird 


ie 2) [10 len) 09 all) — 20) 3.85) 


+f(a) aA Jı (Aa) — f(b) 6A Jı, 94 ne @ JA, e) de 


d 
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Setzt man in (9) zur Abkürzung 
a 


g9(A)—= ca Jy(A,a) —dıb Ju (A,b) + 6zaAJ; (A,a) —d; bAJ,(A,b) + [ w(e) . oJ, (A,o)de - -. (10). 
0 
so ist durch (7) rein formal die Lösung gegeben 


oo 


fo) = j au Al (A) Re N-T1D. | 


yE 


In (10) sind zwei Konstanten durch die Randbedingungen bestimmt. Um die noch freien Kon- 
stanten zu finden, ist eine zusätzliche Überlegung notwendig. * 


Die Darstellung (11) ist nur dann möglich, wenn die Funktion e L beschränkt bleibt. | 


Damit erhält man die Bedingung 


Aus (12) ergeben sich Gleichungen, um die noch freien Konstanten beziehungsweise die Eigen- 
werte zu bestimmen. Das so erhaltene Ergebnis wird dann noch durch Einsetzen in die Dif- 
ferentialgleichung überprüft. Durch die Wahl des Parameters n der Besselfunktionen lassen 
sich auch andere Arten von Differentialgleichungen mit Hilfe der Hankeltransformation 


& lösen. Multipliziert man z.B. die Differentialgleichung g 
| ER 
| 


mit r Ja(A,r) durch, so ist die Lösung gegeben mit 
- co b 
4 ” 1 ” 
= (= | »\| 10-0 Id 0)de— 5b (a6) + Abdı Fi (Ab) 
r | 


a 
N 


+a@Js(d,a) — AdzaJg(d,a)|AJ, (Ar) dA. 


Im folgenden soll nun die Anwendung dieser neuen Rechenmethode gezeigt werden. 


3. Transformation der Differentialgleichung 


Wird eine Kreisplatte durch rotationssymmetrische Kräfte normal zu ihrer Oberfläche 
beansprucht, so genügt ihre Durchbiegung der Differentialgleichung (1). Das in Ziffer 2 ange- 
gebene Verfahren führt bei konstanter Biegesteifigkeit N auf die Durchbiegung 


so=y|[s@ ER! NÄA ERBE LEE l; 
wobei | ( 


g(A) = [ro "Js, E)do+aa Jlha) + aAJı(,a)+ aM Ild,a)+maAJı(A,a) (14) 


und a der Radius der Platteist. Für eine Kreisringplatte mit dem Innenradius a und dem Außen- 
radius b wird 


b EN 
sa= [Pl FT (&)de+abIlhb) had (da) + Abd (Ab) a. 
draht) HRaha)takbhh|  —— 


Aus den Randbedingungen und der Bedingung (12) lassen sich die Konstanten bestimmen. 
Wird eine Kreisplatte veränderlicher Biegesteifigkeit durch eine „Stufenplatte‘‘ ersetzt, 
so ar sich e En Platte ebenfalls das beschriebene Verfahren anwenden. 
ie schon Marcus?) gezeigt hat, läßt sich die Differentialgleichung (1) aufspalten in di 
beiden Differentialgleichungen 2. Ordnung 5 er an 


AMW)=—pt), 2 (16) 
1-0, er an 


Be 8) H.Marous: Theorie elastischer Gewebe. Berlin 1924. 


2. angew. Math. Mech. S ’ 
Bd.32 Nr. 2/3 Febr./März 1952 Jung, Ein Beitrag zur Statik der Kreisplatten 49 


In (16) ist 
Mm 
Mer m) 
mit 


dr " gmdr 


een 


m dr? rdr 


ee 1 al 


Die Lösungen von (16) und (17) sind gegeben durch 


M(r) = [na NENNE IE a a ae (18), 
fir) = ! BLIEB (19). 

In (18) und (19) ist 
na) = [ro BEN Baader (20), 


%+1 


M 
= 2 | SAN do uchlo)-tanid.l JB . (2). 
Für M(r) wurden in (18) dieselben Voraussetzungen gemacht wie für f(r). 


4. Unterfunktionen von Belastungen und Deutung der Integrationskonstanten 


Wird eine Kreisplatte durch eine konstante Flächenlast von der Größe p, kg/cm? innerhalb ; 
eines Kreises r<a beansprucht, so lautet ihre Unterfunktion mit Bild 1 


a 
uA)=m j Nena)... (2). 
ö 
Aus (22) erhält man durch einen Grenzübergang die Unterfunktion für eine Einzellast im Ursprung 


i 0% P 
q2 (A) = Bu Paz Jı (A, a)—= Di er ORTE ec (23), 


wobei 
P=limp, «r 
a>0 


Bild1 Bild 2 


Für eine entlang des Kreises r=b angreifende Linienlast wird mit Bild 2 
in mr Haar = RR RE (24). 


Die Integrationskonstante cı stellt somit die Auflagerkräfte in-kg/cm dar. - - - 
Für die Kreisringplatte ist die Unterfunktion durch I) gegeben. Ist die Kreisringplatte am 


% -  Außenrand frei, so ist in 15) 
ee ee = 0. 
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Wirkt entlang des Kreises r—=b-—-e eine Linienlast von der Größe a und auf dem Kreis 
r—=b+-seine Linienlast von der Größe +p, (Bild 3), so lautet mit (25) hierfür die Unterfunktion. 


4(4)= — pı(b—e) I, (A, b—E) + Pa (b-+E) Jo(h,b+E) » : . . - - (25a). 
Macht man in-(25) den Grenzübergang 
4 0) = im q,(%) 
und setzt 
lim p 2e=m, 
e—0 
Ipı 2r (b—e)|= |p*2n(b+e)| =p'2n, 
so lautet die Unterfunktion für ein Ringmoment 
GA) m, 5ER HA,d).-. oe DE a ee (25). 


Die Konstante c, stellt somit ein Ringmoment dar, das entlang des Kreises r—b angreift. Ist 
die Platte am Rand r=a Sn so ist c, das Ein- 
b — 
spannmoment. 
-My | Mn 


lim 2em,=M 
€E>o0 


mit 


Bild 3 Bild 4 


Wird eine Kreisplatte durch Momente dritter Ordnung, die entlang des Kreises r—=b an- 
greifen, nr so lautet hierfür die Unterfunktion 


1 MIR ren (26). 


Die Momente dritter Ordnung sind in Bild 4 dargestellt. Sie bewirken einen Sprung in der 
Plattenkrümmung und einen Knick in der Tangente. 

Die Konstante c, stellt damit ein BInISERE dritter Ordnung dar. Bei der eingespannten 
Platte ist 4, = 0: 


5. Formeln und Beispiele 
An Beispielen soll nun gezeigt werden, daß die in Ziffer 4 auftretenden Integrale streng 
ausgewertet werden können. Zur Bestimmung der Funktionen g(A) verwendet man am besten 
eine Tabelle, die die wichtigsten Korrespondenzen enthält. Die Unterfunktion”wird mit g(A) 
und die Originalfunktion mit f(r) bezeichnet. Zur Berechnung der Unterfunktionen wird die 


Rekursionsformel 
[/ 7} 


| j rn+1 J(A,r)dr = 
b 
herangezogen. Damit läßt sich »* auf die beiden Integrale 


a E; md (A r)dr 


ee 


\ f J,+1 (nA) BED Lit, 4) j 
a A 
D 
zurückführen. Aus der Beziehung 
ji ER 0 r>o, 
Ira Jar) Hr d= (ru 1 RE, (27) 
Fat ar) dan 


erhält man die Unterfunktion für 
E ID=(@—ry 


rer ur eh. 
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Die Unterfunktionen werden zur Auswertung der Integrale mit der Beziehung 


2 
eh = Jr) + Jrıldr) 


umgeformt. Wird eine Platte durch Ringlasten bzw. durch eine Einzellast im Koordinaten- 
ursprung belastet, so treten in der Unterfunktion Glieder von der Form 


(a) ex J, RI) 


ER N (28) 
auf. Durch die Differentation 
Bear) (A, EEE 
A dA 4 1%) 
läßt sich (28) gewinnen und auf die Integrale (27) sowie auf das Integral 
7 Jo (A,a)—1 = .0 r>a 
| IT IA. , (29) 


ö In — r<a 
a 


zurückführen. 

In Tabelle 1 sind die in den Beispielen verwendeten Unterfunktionen zusammengestellt. 
Um die Unterfunktionen schnell umformen zu können, sind in Tabelle 2 die hierzu notwendigen 
Formeln zusammengestellt. Integrale zur Rücktransformation sind dann noch in Tabelle 3 
gegeben. (Zilfer 8) 

Für den durch Bild 5 gegebenen einfachen Belastungsfall wird 


(A) yo 2) +2 I,(2A) +24 2A) +2, MR I(QA) + 2 28 J, (2). (30). 


Die Bedingung (12) führt auf 


Pt a=0, 
Pt2% +25 —=0. Bis 
Geht man damit in (30) ein und formt die Unterfunktion um, so wird 


fo=y7 ERDE 5 | Ne 2 Ar) ga 20, | I(22)- 900,7) (31). 


0 


Wertet man die Integrale in (31) aus, so ist 
= 1 Po 2\2 63 2): 
IO=y aan) 4 AP) Fa EB EEE Ten (32). 


Aus den Randbedingungen läßt sich dann c, und c, bestimmen. Der Vorteil dieser Rechen- 
methode besteht darin, daß die Lösung die Randbedingungen der Platte noch frei läßt. Ist die 
Platte auf dem Rand r=2 elastisch eingespannt, so lauten die Randbedingungen hierfür 
d 
r—=2; 0; Sum). 
Die Federkonstante der Einspannung ist x. 
Aus den Randbedingungen folgt sofort c,=0. Um co; zu Baehnen wird die 2. Randbedin- 
gung noch etwas umgeformt. Es ist 


df nz 1 a. 


dr.“ dr?  mrdr BR (33) 
mertadf > ER i 
mr de dr? 
‚ mit | ee 
Setzt man in (33) noch m=3, so wird | 
| ee 
2(3+20) 


4 
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Die Durchbiegung der Platte ist damit 
— 2)? 
Für den durch Bild 6 gegebenen BL. lautet mit (24) die Unterfunktion 
3 36 
ph a ay+ Pen + Freh+ ZB : : 


en 


Die Durchbiegung wird nun mit Hilfe eines RAREnbIaEtEn ermittelt (s. Ziffer 3 Aus den beiden 
ersten Spalten erhält man die Gleichungen 


2pP+34=0, 


4 
en Dr | 
en 


Bild 6 
3 Dann geht man mit den daraus bestimmten c, und c, wieder in das Rechenblatt ein und erhält 
ER für die Durchbiegung 
a N-f)=2Pp0—r)+2pin + a 70-)+4 2<r<3 
2 (36). 
N -f(r) = 279 —r) —2p4—n)+5 - 3 t2rInz ne 0—-)+u4 0<r<2 
Für r=2 ist 
A —=h; 2> Fi 1 —=f; Dir 
Für den durch 3 7 Belastungsfall wird die ann. für das Monat 
4c 
| a PR ah) + Zar) — EAU Er Q) » . . @N. 
Aus (37) erhält man a Rücktransformation 
16— r? 
. Mr)=m 4 +um+a+G a A EEE 5 A: (38). 
Mit (38) erhält man die Unteffunktion für die Durchbiegung 


1 


2 
92 (A) = =, | Ei (16 —r?)-+d, In T +%+ 4) r I (A,r)dr. 
i 


+ 403 Jg(4 A) + dy Jo(A) 
SR eg! 
In (39) setzt man ep und wertet die Integrale aus. Unter Verwendung der Tabelle 1 ist 
8 
ad antrat ad + nl y—B 122) 
men, a en ) R 


p 15 d 
ng 47 


\ 


h rn) 


4 - . 
1 0 BER 
r N, / (4 16) + dm +%+ 2) r Ju, r) Re (39). 


. (40). 
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Aus (40) erhält man die Durchbiegung, wobei O(—=c,-d, ist 


p 16 —r? 
NA) = 4 (16 PP + [16p —12d,+ 80) 3 
j Tr dı dı % 
+ [169 —4dı+8C+4e;] hr a ea In» 2 sr <s4 
N P AE_ ya. 16 —r? r 
Aal) =, Er) + [16p— 120, +80 —5— + 1169 —4 d+80+45]1nZ 
dı SE ae : 7 4—r? 
a N ee a rt a+0] 5 
7 d d 
— gr-1184+0|n5-74u-m+grn;- l<r <2. 


Die Rechnung läßt sich ebenfalls mit Hilfe eines 
Rechenblattes durchführen (s. Ziffer 8). Ist 
die Platte anihrem Innenrand r=1 frei, so ist 

ı—=0 und die Konstanten lassen sich aus den 


Randbedingungen 
Ash =0; m —=0 
v=] m,—0O 
bestimmen. 


Für das durch Bild 8 gegebene Beispiel sind die Abmessungen durch Tabelle 4 gegeben. 


Tabelle 4. 
r | 0 1 | 2 | 3 


Setzt man für das entlang des Kreises r—=2 angreifende Moment 


M=a+Bß, 
so wird die Durchbiegung 
ee). 
I0)= 1, | 03010" “ —0,3564alnr =,0:fürr>1 
0 5 
4—r? DSE t 
+ 3,456 a — + 0,3456 a In, — 0 für r>2 
pa 
pn! —0 für r>2 
36 =" 1 16 Bm 73,092 In | 
len 


Aus der Randbedingung 
m. — V-JSür ‚ro 
erhält man 


An der Stelle r—=2 tritt ein Sprung in der Momentenverteilung auf, aus dem sich « bestimmen 
läßt | 
RE & —= 0,96606 M. 


- Damit ist nun an verschiedenen Beispielen gezeigt, wie die in Ziffer 3 beschriebenen Verfahren 


anzuwenden sind. 
| | 6. Kniekung der Kreisplatten 


Wird eine Kreisplatte durch rotationssymmetrische Kräfte in der Plattenebene belastet, 
so läßt sich der Spannungszustand z.B. mit dem Grammelschen Verfahren berechnen. Über- 


schreitet die Belastung einen bestimmten Betrag, so beult die Platte aus. Für eine Platte kon- 


ke ET er 


54 Jung, Ein Beitrag zur Statik der Kreisplatten Ba: regt 


stanter Biegesteifigkeit wurden die Beullasten von Meißner‘) bestimmt. Federhofer?) und . 
Gran-Olsson®) geben Verfahren, um die Beullasten bei veränderlicher Biegesteifigkeit zu 
berechnen. Das von Federhofer angegebene Verfahren wurde von Egger‘) numerisch aus- 


gewertet. 
Es soll im folgenden nun gezeigt werden, wie mit Hilfe der Hankeltransformation auf 


einfache Weise die Beullasten bestimmt werden können. Im Anschluß daran, wird einer Be- 

merkung Grammels®) folgend, ein Iterationsverfahren aufgebaut, das auch auf Kreisplatten 

mit sprungweis sich ändernder Biegesteifigkeit anwendbar ist. 
Wird eine Kreisplatte auf ihrem Außenrand durch die Belastung pkg/cm beansprucht, 

so wird durch diese Belastung ein homogener Spannungszustand hervorgerufen. Für diesen _ 

Belastungsfall lautet die Differentialgleichung für das Ausbeulen?) 


Au Sau PAR Fe Bere: 


Multipliziert man (41) mit r.J,(Ar) durch und integriert zwischen den Grenzen r—=0 und r=a 
partiell, so geht (41) über in 
af urdlar) dr+caFlAa)+0gaAJı(la)+c; aA? Ju(Aa) + aA Jı(A a) 
0 
up a nr ee (A 
-E 2 [uerdandr+dahaa)+hardılaa) ? 
ö 


Aus (42) erhält man die Durchbiegung 


=] Ga, Aa)+%aAJı (Aa)+czaA2 J,(Aa)+ aA J, (Aa)] X dA . .(43). 


Da keine Kräfte normal zur Plattenebene auftreten können, so muß in (43) 


210 
sein. f - 
Ist eine Platte auf dem Rand r=a eingespannt, so ist mit Ziffer 3 
G3 = 0 und Cq = 0 . 
Damit geht (43) über in , 
 f &aJı(ka) Kr 
um ER Le er EN (44). 


ps az 
ATTEN 


0 
Mit der Bedingung (12) erhält man aus (44) 


ESR —. 


die transzendente Gleichung zur Bestimmung der Beulbelastung. 
Für eine frei aufliegende Platte ist 


Geht man mit (45) in (43) ein, dann wird 


aaa ne] 2 
ng a ie © (46). 


uP 
5 ee 


#) E. Meissner: Schweiz. Bauztg. 101 (1933), S. 87. 


5) K. Federhofer: He i : : 3 : 
ei IH Ei ve S SE aa 11 (1940), 8.224; und Sitzgeber. Akad. Wiss, Wien, Math. Naturwiss, 


6) R. Gran Olsson: Ing.-Archiv 12 (1941), 8.123 
?) H. Egger: Ing.-Archiv 12 a1), 8 Rz 


8) C.B. Biezeno u. R.G f i ; vr 
- #)- Siehe auch Fußnote Re SER Be By 1gaR 


’ 
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Die Bedingung (12) führt auf die transzendente Gleichung 


uP uP Er up 
rel) 


Ist die Platte auf dem Rand elastisch eingespannt, so ist 


m—1 
9=—Xr Te (3 » 
Bei der Kreisringplatte kann, da kein homogener Spannungszustand mehr vorliegt, die Glei- 
churfg (41) nicht herangezogen werden. 
Geht man von der Differentialgleichung!°) 


a du 
ser aer] an BEN ea re, 1 ee (47) 


aus, so läßt sich die rechte Seite als Belastung einer Platte denken, sofern Ir bekannt ist. 
. r 


Nimmt man nun ein w,(r) an und setzt 


1 Br, le n.@l du, = 
N Aus era | ee ee ne ZN, (48), 


so läßt sich mit dem in Ziffer 3 beschriebenen Verfahren eine Durchbiegung u,(r) bestimmen. 


Mit diesem «,(r) ermittelt man q,(r) und daraus u,(r). Das Verfahren wird nun so lange fort- 
gesetzt, bis 
%+1(r) Au, (r) £ 

ist. Bei der Iteration setzt man p=1. Mit dem so gebildeten u, berechnet man die Schwarz- 
schen Konstanten 


u 
Gi+1 —ı U; YWi+i dr ee ee ar. 7£ (49). 
7% 


Damit erhält man einen ersten Näherungswert für die Beulbelastung 
4 A;; 
ae En et) 
Gi+ı 
Für eine Kreisringplatte, die auf dem Außenrand durch p, und auf dem Innenrand durch >; 
belastet ist, wird 


_Bnzon , Bpnrn (51). 


0, = — —ı nn nm 0 — 0.0.0.0... 
f n—n 


- 


In (51) führt man = = ö ein und erhält 
” 


ra —ör; (1-6) rar; 


= — Pa ar, + Pa a ee (52). 
Geht man mit (52) in (48) ein, so ist 
r ie n—öri a) rar; duo " 53 
Io (r) TI gr FR (( ST EREBe ; — ; A) —) Er ER TEE ( ji 


“Zur Durchrechnung eines einfachen Beispiels wird 


2 EN | 


gesetzt. 


Geht man damit in (53) ein und wählt die Näherungsfunktion 
n Up = 16 FF 12, 
dann erhält man 


64 
a) 45 p. 


Die erste iterierte Durchbiegung wird damit. 


64 16—r?)? 3 En an 
Ben u tt ang engt 4 26, M+am;| sn 


10) Siehe Fußnote 8 von 8. 54. 
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Bestimmt man aus den Randbedingungen 
m = 
r—=4, m —=dO 
6, und a, so wird 


16 —r?)? ER FET, 
u(r) = 15V ee + 0,6603 (16—r?) — 4,2606 In Te In r on (55). 


Bildet man mit «, und u, die Schwarzschen Konstanten, so ist 


» = 0,2242N. e: - 
Aus der bei Biezeno- Grammel!!) angegebenen Kurve liest man ab 
p2* = 0,223 N. e“ 


Für den durch Bild 9 gegebenen Belastungsfall, 
der wohl erstmalig PU CHSeIEES wird, erhält 


Fr 12 = ae 
Pe € 8, = 
- x Pi 7—r? 1 
Bild 9 re = a 
Mit der Ausgangsfunktion 
w=16—r | | (56) 
wird F 
1 
|#° 2<sr s4 
AAu = 
RE 1 <r<2 £ 
|® 3 P = i 


ah (56) erhält man 


: SDR 32 1... 
= u zn) +S It zr Int 7 (16—r) 


192 3 4 
+ & In r ’ 2 sr s 4 
BER! (412)? u | 
rer a 0,31906 (1) — 2,6815 m 4 Ir nz 
wenn] 1X 
Mit den Randbedingungen | 
du > 
r—=1 er dr =D ’ 
re m.=0; : 
berechnet man sich c, und a. Der erste Dr für die Beulbelastung wird 
Der Quotient — — ® schwankt zwischen den er 
2a, Br; % 


so daß d 
er. En re Näherungewent hei einem ern Iterationsschritt sich nur unwesentlich 
1 N 


11) Siehe Fußnote 8 von 8. 54. 
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*. Die dieke Platte 
Bisher wurde angenommen, daß die Plattendicke klein gegenüber ihren sonstigen Abmes- 
sungen war. Diese Einschränkung soll nun aufgegeben werden. Wird eine Platte von der Dicke } 
durch rotationssymmetrische Kräfte beansprucht, so lassen sich die Spannungen und Dehnungen 
mit Hilfe der Loveschen Verschiebungsfunktion angeben. Bezeichnet man mit ® die Lovesche 
Verschiebungsfunktion, so genügt diese der Differentialgleichung'2) 


ERSTER ae (57). 
wobei 
ee 
ar 9 
ist. Die Verschiebungen und Spannungen sind 
a m od 
ml Hr °’ 
m— 1 Mm. ED 
Vesper 
At a) 
na 349 nn) 
ERTCH 5 a ea (58) 
nad r or Ä 
2(2m—1) „2? m @®d 
er) FR Are] nn 
2(m—1) , 2 m ®2® 
ee) z or (40 - m—1 a 
Macht man die Voraussetzung über ® 
| 05 nd 
SS rer ee 


so läßt sich auf (57) die Hankeltransformation anwenden. Setzt man 
P—= | G(rz)rJ,(ar) dr, 

0 
so geht (57) über in 


0:9 ep 
4 —_ EM) u ee ee een Ve 59). 
z 2 022 7 02% Y e 
Die Lösung von (57) ist bekannt und wird in der Form angesetzt 
9=(A(@)+ zB(A)) &ind2+(O(A)-+2D(iA))Eofiz. . . : . EN EN CE (60). 


Mit (60) erhält man die Verschiebungen und Spannungen in der Form 


w=| (2BACof A2-+2DA SinAz)AJ, (Ar) da 
! | 


5 | [A SinAz-+ Bez SinAz +0 EofAz+ DzCofA2]4°J,(Ar) dA, 


0 


= 907 
SE Te er ee): 
= —, IC 2)(B SinAz + Dvd) | (61) 
0 
— Am (A C0fA2+ BzCoj A240 GinAz+ Dz Sin A2)]ArJ(Ar) dA, 


Fa —;, fe# Gofis-ED min). 
an 4 


I mA(A SinA2+ Bz&inAz-0 CofAz-+ DzCofAz)]A2J,(Ar) dA 


"Für den durch Bild 10 gegebenen Belastungsfall lauten. die. Randbedingungen 


Fl für 20, o;=f(r) für z—h, 
Re zes 220.40 --für zen. 


. 12) Siehe z.B. C.B. Biezeno u. R. Grammel: Techn. Dynamik.: Berlin 1939, S.35 ff. 
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Geht man mit diesen Randbedingungen in (61) ein und setzt 


oo 


g()= IR) “0 JA 0) do: = Aa ten a ee 62), 


0 
so erhält man 
(m — 2)8 h SinAh a) na 
26m: R(SmAn—am) I» 
_ m—2 SinAh-+AhCojAh i 
= 2Gm 3 (SineAn — en) I’ 
_ m—2 SinAh+AhCojAh a) 
Gm: 3(SinsAn—am) I» R 
“m—2 h SinAh a) 
2 Gm A (Sin2Ah— 2 12)7 
Die auf diesem Wege berechneten Funktionen stimmen mit den schon bekannten überein?®). 


Geht man mit (63) in die erste Gleichung aus (61) ein, so erhält man, wenn der Integrand 
mit Ausnahme von g(A)-A J,(Ar) in eine Potenzreihe nach A entwickelt wird, für das erste Glied 
dieser Reihe 


A=- 


= 


D= 


1.f ga > SR 
u | rnana+-- für :=7: 
(1) 


wobei N wieder die Biegesteifigkeit einer ‚„‚dünnen“ Platte ist. In (63) muß damit für g(A) die 
in Ziffer (3) gefundene Unterfunktion verwendet 
werden. Berechnet man sich noch weitere Glie- 
| 6,:f(r) der der Reihe, so gibt das zweite Glied den Ein- 
EN fluß der Schubkräfte auf die Durchbiegung. 
Damit ist nun ganz allgemein gezeigt, daß 
ee! die technische Plattenbiegungslehre die erste An- 
Tr EEE Tr 3 näherung an den genauen Spannungsverlauf gibt. 
Der Vollständigkeit wegen?°) sollen hier auch 
Bild 10 noch die Spannungen für einen beliebigen Be- 
lastungsfall angegeben werden. 


Mit den dimensionslosen Veränderlichen erhält man 


N 


zZ r 
A=%, a Bram 


Sinz-+xCojx 


Io | 
.= 7 | |tSinnu—nxGoln“ Sims 


0 
; Einx 
+ nr? Sin n% Sinn) g(x) dk, 


Sinx 


er 
| 


0 


z Sin x+xCoj x 
—n Sin rer | 2J, (0x) g(x)dx, 


INNE Sinx+ x Cof x GE, 
= ze | rei n°a+nx% a ne ee Poser gr, (64). 
0 


kei ' Ein x 
[2&0| nx+nx Sinn] a xJ,(o x) g(x) dx 
: | 
a er 


Pr 


») H. Jung: Ing.-Archiv 18 (1950), 8.178. 


Z. angew. Math. Mech. 


Bd. 32 Nr. 2/3 Febr. März 1952 Jung, Ein Beitrag zur Statik der Kreisplatten 59 
ar Sin® 
ın D FEN AR BREISEFTE, 

| |e NT N any (64) 


1 
x x (er) g)dr+ smrun 


Dabei ist zur Abkürzung in den Ausdrücken für o, und o, gesetzt 


oo 


Men | IIom —2) Sinnxz(Sinx+xEoj x) —2(m—1) Co nx Sin] 


0 


+ nl &j nx%(Sinx+xCoj x) — x? Sinnx Sina]! x Se) Bu IR: 


8. Anhang 
Tabelle 1 


fr)=0 r<b und r>a 


fer) = 3 g(A)-AJı(Ar) a) = [nr Jar) dr 
b 


fe) =e 2 E Bam 0 ) 
r? une = Zn J,(Aa) + = J,(A b) 
& er J,(Aa) H+ —z—- = Isla) — Zee had) I 2 u 
2 2094 109 
ad a)— Zn db) — b(a? — b?)-— Jı(Ab) 
SE a © 1.00) SE (oe — u a0) u N) 


2 
z Tu) — ER Iuar) — (a8 — 89) Iulad) — Sy (0 — De (ab) 


— 4 (a — BR AB) 


EN, Dy Et Di nt J,@Ab) 


En AL = EN Tl a) + ee mt J,(Ab) + ae ,(2) 


+ m Ina) 4% En Ha) 5 J0(2 3) 


Tabelle 2 


BE 2 


. er BAU a)+ 8 Sm J,(Aa) + 8 SE 2 (A a) 


Tank) —ı\ 3a Br 
a ae 422 Jo(Aa) ‘ 


> Tung, Ein Beitrag sur Statik der Kreisplaiten na so au IB voran aan 


Tabelle 3 
% 0 für r>a \ 
I, +1(a)) Io (Ar) a Hai = er) x ag | 
J u+l 0 TEEN TREE für rsa 
) a re 
rar a3 
et 1 ara = 5- für en 
A. “ S 
25 i 0 für r>a 
f de == , 
} 4 2 "a pa nr 
7” für 0<r<a 
En [He@n I, 21a) ar _ 3 
SITE ER für 0<r=a 
2r 
0) für 0<a<r 
Bu 
1 f 
Rechenblatt 1 i 
J,(Aa) I, (Aa) 1 d2 Ju(Aa) —1 { 
72 78 aan 7 
4 


ey —ı ; 

z 

. 

; 

| | 

9c 9c h 
232": (32) +53 82) | 
{ 


pa JB) —1 
ı dr 2 


1 „ Jar 1: CM 

— m 2 ———__ 

7a 90) or Jah) — Ge han +E « 2 + ne J. RE = Ja) 
ee En er 


N =270-r) +2pn +2 gra-+ zo-m ta 2<7<3 


Nf,(r) =2p9—-)—2p4—r)+2 Bring Harn + +20 -9)+4 0<r<2 
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ne ckenblakt 2 


= Re 5 Jılal) 
Zn ar 
— Jı(ah) 
Ar ER (4) 
+anan 
10) 
et Bu: (A) 
pe J,(44) + = AD Bla Re I = Su TotA) 3 


DERSA 


16 — r2 1 2 Sr = 
Mr) & Po 4 z TS n un ++ 
3 Isla) | A I ACEN © SAoa)-—I 
92 (4) 73 72 FE dr el 
8 8 16 16 
En J,(44) ve I, (44) a Jul) GE Jo(&) 
Ar 4d, dd I al—1 
(#4) J2(4 3) — 32 Jo(42) HT : 
EX) 
Sram + 2 Jo(4) 
4 
SERrE an) 
una | hen | hen | er 
Er = J,(@3) em Klar n BACH) 
d, 3d, CAR _ du @ @N—1 
= 377024) 37 722) zu 2 ANEBFN 3 
2 | 
d, In — d a 
ET, 21 a: et 
eh = — 2 | = ——Hel29) 


0 C [9] 
Be -Zhen | - her 


9, (4) = Pe 7,01 + [116 — 124, +8011,9,42 + [16m — dd +80 + 4a) ga /leN 


— gm 


07,0) — = 70 — 2,1931 d, +0 oe B od, +0]; 1 nan 
en ae ad | 


I. NA = 42 6 — 2 + [16 pg — 12d, + 80]- A l16m-4d, +80 +40] 


ET. pers 


N,f(@) - fi 2 | 5 Po — 2,1931 dı + A Zr B 9 — 119 4, °+ olim a 


Aa-m+& amt Is 


z Eingegangen am 7. Mai 1951. 


\k 
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Eine Fehlerabschätzung für die Verfahren von Adams 
und Störmer') 


Von Johannes Weissinger in Hamburg 


Nach der in [4] entwickelten Methode werden die von G. Schulz angegebenen Fehlerabschätzungen für 
das Adamssche Intewpolations- und das Störmersche Extrapolationsverfahren verschärft. 

Following the method of [4], the estimations of the error are improved that were given by G. Schulz 
for the method of interpolation (A da ms) and extrapolation (Störmer). 

Suivant la methode d&velopp£e en [4] les estimations des erreurs, donn&es par G. Schulz, pour la methode 
d’interpolation d’ Adams et la methode d’ extrapolation de Störmer sont perfectionn des. 


IIo merony, paspa60TaHHoMmy 8 [4], yıyamarorca naHuzıe T‘. ll yısıem omeHkH morpemmocreä 
AA HHTEPNOJAIAUMHOHHOTO (AnaMc) HU 9KCTPANONAUHOHHOTO (ÜlTöpMmep) cmoco0oB. 


1. Einleitung 


In einer früheren Arbeit [4] wurde gezeigt, wie man die v.Misessche Fehlerabschätzung [1] 
für das Extrapolationsverfahren' von Adams durch einen einfachen Kunstgriff wesentlich ver- 
schärfen kann. Dieser Kunstgriff besteht darin, daß man nicht wie bisher die exakten Fehler- 
gleichungen in Differenzenungleichungen verwandelt, indem man (im allgemeinen unter Verlust 
an Genauigkeit) den Betrag einer Summe durch die Summe der Beträge ersetzt, sondern daß 
man vorher die stets vorhandenen Relationen zwischen den in den Fehlergleichungen auftreten- 
den Integrationskoeffizienten ausnützt und durch Summation (oder andere lineare Kombination) 
der Fehlergleichungen neue Gleichungen herleitet, bei denen der Übergang zu den Beträgen 
keinen so großen Genauigkeitsverlust zur Folge hat. In dieser Formulierung ist der genannte 
Kunstgriff nicht auf den in [4] behandelten Fall beschränkt, sondern läßt sich allgemein auf die 
üblichen Differenzenverfahren anwenden. Das soll im folgenden noch für zwei weitere Ver- 
fahren durchgeführt werden, für die früher eine weniger scharfe Fehlerabschätzung von G. Schulz 
([2], [3]) gegeben wurde: Das Adamssche Interpolationsverfahren zur Lösung der Gleichung 


y'—=f(x,y) und das. Störmersche Verfahren zur Lösung von y"—=f(z,y). Dabei wird die Kennt- 


nis von [4], insbesondere der dort eingeführten Bezeichnungen, vorausgesetzt. 


2. Das Adamssche Interpolationsverfahren 


Der Gedankengang verläuft fast wörtlich genau so wie beim Extrapolationsverfahren, 
so daß wir uns kurz fassen können. Wir benutzen für Integrationskoeffizienten, Quadratur- 
fehler usw. dieselben Buchstaben wie früher, nur mit einem Querstrich versehen. 

Die Näherungslösung % berechnet sich aus 


rn hun hun hn—hd if a SAT a ER (2:1) 
o= = 
mit 
r | s 
% o\- ie 
%=(—1)! >63 ; a ıjewrn u Hal) du «IR ERBE 
o=Q —l 
so daß die Fehlergleichungen lauten 
| an=athls name) t+% 222.0. (23) 
und nach Summation über die Indizes s, i—1,...,i—1@W21+r—1, 1>r): 
r+li_ = U 
&n=änt h 2 bel fire —firi-e] +2: u 
mit r = 
- tat '+%, 0<o<r—1 
d.= 51 Zeh 
Hate EI HISOo<IHtr 


Daraus erhält man wie früher eine Fehlerabschätzung der Form 
a Dee 
u—yl<se + 1) NZ BR Er 
*) Über die Ergebnisse dieser Arbeit wurde auf der Jahrestagung 1950 der GaMMin Darmstadt berichtet.‘ 


Fu 
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wobei der Anfangsfehler &* durch die Forderung 


> |e; |, ed aeer ln ; reed 7 
definiert ist und z.B. nach G. Schulz [2] Horte werden kann. Der Süadrafürtekler q ist 
5 DT ee E . (2.8) 
und die Zahl A* berechnet sich aus 
r—1 


2 ih (In, bs t or +04. . (29), 


während 2, die positive Wurzel der Gleichung 
z1(zirl 1) 


Kh=-- = UNO) Na en ae ee 210) 
> |d.| „itr—o 
o=0 
ist. Es wird 
dx 1 
Ted tn N ten Aa er et a N 21 
dz |.-1 Dh ;% 
so daß näherungsweise S 
ATI ee en (212) 


gesetzt werden darf und der in [4] ausgesprochene Satz auch für das Interpolationsverfahren 
gültig bleibt. Dagegen haben die Funktionen x(z) genau so wie die entsprechenden Funktionen 


015 020 


"Bild 1. Die Funktionen x.(2) ( ) nach (2.10) mit /=[7, are mit den entsprechenden 
Funktionen ( u... ) aus [2] fürr = 1,2 ‚3, 4 


von G. Schulz [2] eine negative Zweite Ableitung, so daß in (2.12) das Zeichen nicht durch < 


ersetzt werden darf (vgl. Bild1). 


RN IN \ Y gr 


Ä Die Zahlen @, A* können der folgenden Tabelle entnommen werden, die zum Vergleich 
auch die entsprechenden Zahlen A von G. Schulz enthält: 


‚| (Ü+1)(4*—1) | Anach [2] | 


1 1.0000 1.0000 0.50000 
2 1.1667 1.1667 0.08333 
3 1.3333 1.4167 0.04167 
HH A 4 | 1.5444- 1.7861 0.02693 
5 0.01875 
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ie A ä it r wachsen als 
Ein Vergleich mit der Tabelle aus [4] zeigt, daß die A wesentlich schwächer mit 7 
die Zahlen A des Extrapolationsverfahrens. Daher ist der Gewinn an Ge Er 
neue Abschätzungsformel für das Interpolationsverfahren längst nicht so groß wie e 
Extrapolationsverfahren. Das zeigt auch die Durchführung der Abschätzung für das Beispie 


’ y—% de 
— ; y(0)=1. 

’ Tyta 
Wir wählen r—=4, 1=7 und können für &* die bereits in [2] berechnete Schranke für 
lals--- ‚ |&01 


er =:5673 10% 


nehmen. Der Tabelle entnimmt man 4*—1,193, also 27 =1-+0,006-1,193=1,00716, während | 
in [2] z=1,011 ist. Mit 2° —=2,041 sowie 9=1,061-10% bekommt man dann 


1,061 - 10 
Ey — 0,673 - 104 2,041 4 Goarg PM 


—=1,374 «10% + 1,543: 104 2,92-10, 


also keinen wesentlichen Gewinn gegenüber der in [2] errechneten Schranke 3,47-10. Bei 
größeren Schrittzahlen bzw. größeren Werten von Kh wird natürlich die neue Fehlerschranke 
erheblich mehr unter der alten liegen. — Vergleicht man die neuen Schranken für das Extra- 
und das Interpolationsverfahren, so ergibt sich ein Verhältnis von 45,1:2,92 #15:1, während 
in [2] für die alten Schranken ein Verhältnis von 574:1 errechnet wurde; d.h. das Extrapolations- 


verfahren ist doch nicht so ungenau, wie es nach der alten Abschätzung den Anschein haben * 


konnte. Er 


3. Das Verfahren von Störmer 


Die Näherungslösung % der Differentialgleichung y"=f(x,y) berechnet sich nach der 
Formel 


.. 


x ?: n r = 
Pr —=MW 2 ap fu—l I of ker. ee 
0o= 0= 
mit 
J 0 
x | 
ED do, ne ex a uu+1)... uto—I1)du . . . . (8.2), 


insbesondere «,—=1, a—=0. Die Fehlergleichung lautet dann 


Par = 2 late ti ker, V Eh en (3.3). 


Wir addieren nun nicht wie bisher eine feste Anzahl /-+1 dieser Gleichungen, sondern summieren 
über k=i, i—l,...,r: 


i N ir 

van Ve= BI be(fie— fie +20: Er ua 2 0 BER en (3.4), 
wobei auf Grund der leicht zu bestätigenden Relation 
| Digi 1 2 Ua A ER a (3.5) 

0=0 _ 
gilt S * 
GH +"'+0, 0<o<r—1 2 2 PR 
b=1J1 k r<oe<Iier (füri>2r) .. (3.6). 


N tt "+, diort1<osi 
Dabei ist >2r—1 vorausgesetzt. Allgemein kann man schreiben 


be 00; mit... 0 -Er ol Senne (3.7). 


b, hängt also auch von i.ab. Bezeichnet wieder K die Lipschitz-Konstante von S(&,y) und q 
eine obere Schranke für die einzelnen Quadraturfehler |g, |, so folgt: 


Iiral S lil+rl+ Iemıl + Et 3 IB PET RER TE ein. 2. . (88), | 
\ F o= “ £ ; h wer 
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Daher liefert ‚jede Lösung der sun 


Ai A. EX, +EAtERZ, Ib] Xi + @—r-+1)g, esor.r rl3 


mit der Anfangsbedingung 
pr, N (3.10) 
eine Fehlermajorante. Jede Lösung von (3.9) genügt auch den Gleichungen, die entstehen, wenn 


man je zwei aufeinanderfolgende Slelchungen (3.9) subtrahiert. Wählt man s=27—1, so lauten 
diese Gleichungen 


LGa=2%— X 7 + Kb) + (dl) Kit + (ld bl) Kr} +2] 
= 2X, —- X 1 + KR 1b — K)+ + |] es 8.11). 
v=2r, 2r41,. | 
Ist nun X,—=E; eine Lösung von (3.9), (3.10) und X;—=E,; eine Lösung von (3.11) mit der An- 
fangsbedingung 
E, 2 Es, se DR EM ee A NE VE X N A (3.12), 
so folgt offenbar 
Del Dei, BZ, Well, 2. een... (3.13) 
d.h. die E, bilden eine Fehlermajorante. 
Setzt man zur Abkürzung 


DB KLAR = Id) — |bi 9 E— ıR DEN Re Es OR (3.14), 
und bezeichnet mit z eine beliebige Wurzel der Gleichung 


R 
Pe Dr EEE RENbEFTeN N (3.15) 
o=0 
und mit (eine beliebige Konstante, so ist 
ei 
Be 02 rar ren BER erde .(3.16) 


eine Lösung von (3.11). Wegen PUZZERZB=ZERb= ER und P(@)= 


gibt es wenigstens eine Wurzel z>1. Wir Se mit z, hinfort ie kleinste Zeh von 1 
gelegene Wurzel. 

'Um die Anfangsbedingung (3.12) zu erfüllen, kann man so vorgehen. Bestimmt man 
ein e mit 
a N Sr Re x@:12) 


und setzt ke EN, 
E=e(i=0,1,...,2r—1), Er e[$+Km 12) MEERE EREN) 
so wird (3.12) befriedigt, wenn man 

C>e + ‚, 0(2(4,—1) > Er — Eys-ı=E +2 21%) +rgq=öd (3.19) 
also 


| 2 
(6 — Max le + en 5 en) To ee (3.20) 


wählt. Man kann offenbar auch weniger a 


= 
— ee re ER ER 3.21 
nehmen Ba a die A les 
l&l< ei IE —ä Hs EL (si 1125 RS (3.22) 
erhalten. Setzt man noch | 
£ = eh n IE x. RN, (3.23), 
> 5 
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so ist die Gleichung P(z)—=0 gleichbedeutend mit 


1-1 - re ear (3.24). 
Var" 


Wegen 4’(1)—=1 gilt also näherungsweise für kleine A 


el ti are ee (3.25) 
und dies in (3.22) eingesetzt liefert S | 
7 ö ” Zr | 
lel<e(t + AR) + (1+ aVR)i + ar Ki+arK%—1] . . . (3.36) | 
oder wegen h= a | 
K(2,—% io K(@,;—%) | 
ale IR (2;—%) Es ig au %) 2% a3 15 nn) 1) ee (3.27). 
Diese Ungleichung entspricht der Lösung 
K (0-2 0 YEtw-n) Kla—,) ae 
ya) + a ee Base in. (3.28) x 
der Differentialungleichung | 
\yrı<Klyl+ a a  O  ue (3.29) 
mit der Anfangsbedingung | 
|y(&) |< Yo » y(a)lsy - - - - - Er (3.30), 
k 0) Es,— E32, —ı ‚ . - P £ \ 
wobei e dem y, und re dem yo entspricht. Allerdings sind die Ungleichungen 


(3.26), (3.27) nicht streng gültig, da man wegen A"(1)<0 nicht ,<1-+-4 an Stelle von (3.25) 
schreiben kann. 


Die Lösung (3.28) des Problems (3.29), (3.30) ist recht grob. Die optimale Lösung lautet 
1 Nike 1} 0 \ —YRte-n) K («—%,) 
vn + 2%)? + e = 2 e® 


YK 2 ’K BR.» 


Auch hierzu kann man ein Analogon finden, wenn man beachtet, daß die Gleichung (3.24) | 
wegen der Zweideutigkeit der Wurzel zu gegebenem A >0 eine zweite Lösung z, besitzt, für die ’ 
. man näherungsweise 


schreiben kann. Der Ansatz 


= OGit0, in SE Er NR, (3.33) 
liefert dann das gewünschte Analogon. Das soll hier nicht weiter ausgeführt werden, wie über- 
haupt die Frage nach einer möglichst guten Anpassung an die Anfangsbedingungen nicht erörtert 
werden soll. \ 

no. 
Für r<5 stimmen die Koeffizienten ß, der chärakteristischen Gleichung (3.15) bzw. (3.24) 

mit. den a, überein, die ja wechselndes Vorzeichen haben. Die in (3.24) definierte Funktion A(z) 
geht also in die entsprechende Funktion von G. Schulz [3] über, wenn man in (3.24) die ß, 
durch |ß,| ersetzt, verläuft also oberhalb der Schulzschen Funktion, d.h. die Wurzel 2, ist 

stets kleiner als die entsprechende Wurzel in [3]. Das ist auch aus Bild 2 zu entnehmen. Man 
sieht dort auch, daß A(z) und damit z, im praktisch interessierenden A-Bereich nicht merklich 
von r abhängen (im Gegensatz zu den Schulzschen Kurven, die mit wachsendem r ungünstiger 

werden) und daß die Kurve nur wenig von einer Geraden mit der Steigung 1 abweicht, so daß 
(3.25) in guter Näherung erfüllt ist. “r 


Fa 
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Als „Quadraturfehler“ q kann man die in [3] (17) angegebene Zahl c, nehmen und als 
„Anfangsfehler“ e die Zahl X3,_ı, die sich durch rekursive Lösung der Gleichungen (3.9) ergibt, 
wobei wie üblich Schranken X,2 lei] für ©—=0,1,...,r bekannt sein müssen. Man kann sich 
e aber auch aus der expliziten Fehlerformel (30) in [3] verschaffen, nachdem man diese nicht ganz 
vollständige Formel] richtiggestellt hat. 


015 . 020 


) nach (3,24), verglichen mit den entsprechenden Funktionen 
(----) aus [3] fürr = 1,2, 3,4. i 


Bild 2. Die Funktionen 4 (2) ( 


Vergleicht man nämlich die Fehlerabschätzung (3.22) mit der ähnlich gebauten Abschät- 


zi entsprechender Term fehlt. 


zung (30) in [3], so fällt auf, daß in (30) ein dem Gliede en 


Das liegt daran, daß in [3] die Lösung (28) der Differenzengleichung (24) nicht vollständig den 
Anfangsbedingungen (27) angepaßt ist: Die Bedingung ö,,1—6,>0 ist nicht notwendig er- 
füllt. Um auch diese Bedingung zu befriedigen, muß man (29) durch 


6 2X 
— KH4 
O= Max I+ GE N) 
oder schwächer durch 
6 2X 
5 TIF4+r4, —1 


ersetzen und erhält dann statt (30) eine zu (3.22) ganz analoge Ungleichung, wie auch zu er- 
warten war, da der Gedankengang dieses Paragraphen weitgehend dem in [3] entwickelten nach- 


gebildet ist. 
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Zur Berechnung der Drucktendenz bei Wärmezufuhr 
innerhalb einer isothermen Atmosphäre von konstanter 
Grundgeschwindigkeit | 
Von H. Stümke in Denzlingen/Baden 


Die vom Verfasser vorgeschlagene Methode der ‚‚mittleren“ Zustandsänder ungen zur Berechnung der 
Wetterentwicklung führt beim vorliegenden Problem auf eine elliptische Differentialgleichung für die lokale 
Drucktendenz mit drei Parametern (geographische Breite, Temperatur, Grundgeschwindigkeit), die unter 
der Voraussetzung eines ebenen Erdbodens für den Halbraum mittels Green scher Funktion gelöst wird. In 
der Energiebilanz tritt gegenüber dem Fall der anfangs ruhenden Atmosphäre als neues Glied ein meleoro- 
logisch wichtiger Ausdruck für die Änderung der kinetischen Energie auf. 


The author has proposed a method of ‚‚mean‘‘ changes of state to compute the develpoment of the weather. 


Applied to Ihe problem mentioned in the title, the method leads to an elliptical differential equation for the 
local pressure tendency depending on three parameters (latitude, temperature, and Sundamental velocity). 
If ihe ground is supposed plane, the problem is solved for the half-space by aid of ihe Green’s function. Con- 
trary to the case of the initially resting almosphere a new term appears in ihe energy balance that repre- 
sents the change of the kinetic energy and is of meteorological importance. 

L’auleur propose une methode de changements d’öat ‚‚moyens“ ‚pour le calcul du de£veloppement du 
iemps. Appliquee au probl&me actuel, elle conduit & une Equation differentielle elliptique pour la tendance de 
pression locale avec trois parameires (latitude geographique, temp£rature, vilesse fondamentale). Suppose 
que le terrain est plain, ceile Equation est r&solue pour le demi-espace au moyen de la fonction de Green. 
Oppose au cas de l’atmosphere reposant d’abord, un terme nouveau repr£ösentant le changement de l’ Energie 

. kinetique ei &ant de grande importance met£orologique apparait dans le bilan d’£nergie. 


Ilpennaraembiä aBTOPOM MeTOAX „CPeNHHX‘ H3MeHEHHÜ COCTOAHHA AI BEINHC/TEHUA PasBHTuA 
NOTOAEI HPHBOAHT B PaCCMATPpHBaeMoH IIPoÖJIeMe K HLIHITHYECKOHy MubhepeHiuasHoMy ypaB- 
HeHHIO MIA MECTHOTO Pe3KHMa MaBıIeHuUA C TPeMA IlapaMeTpamu (reorpahuyeckau IIHPOTA, 
TeMlepaTypa, OCHOBHAA, CKOPOCTb), KOTOPOe pemaeTcH, IPeNNoNaTad 3eMHYIW IHOBEPXHOCTb 
INOCKON, AA MOAYUPOCTPAHCTBa TOCpeACcTBoM PyHEIHu I'puna. B IPOTHBONONOKHOCTb K 
TPAETOBAHHOMy llpasıTııem cııy4amw IePBOHAYAJIBHO CHOKOUHOH ATMOchepH, B BHEPTOÖANaHCE 
durypupyeT, B kayecTBe HOBOTO YIeHa, METEOPO.IOTHYECKH BA3KHOE BEIPAIKEHHE JULIA USMeHeHHA 
KHHeTHYeHOf 3HePpTCEH. 


1. Voraussetzungen 


In mittleren Breiten der nördlichen Erdhälfte stellt sich vielfach als Teilerscheinung der 
großen atmosphärischen Zirkulation ein Druckfeld ein, welches von Süden nach Norden abfällt 
und daher mit einer ausgeprägten West-Ost-Drift verbunden ist. Dieser Grundströmung sind 
im allgemeinen Störungen überlagert, deren Ursprung, von der orographischen Entstehungsart 
abgesehen, letzten Endes in horizontalen Temperaturunterschieden zu suchen ist. Es erscheint 
daher lohnend, in einem möglichst einfachen Schulbeispiel die Zustandsänderungen zu unter- 
suchen, die eine solche strömende Atmosphäre unter dem Einfluß einer Wärmezufuhr erfährt. 

"Wie im Falle der anfänglich ruhenden Atmosphäre!) sollen die beiden extremen Formen 
der Wärmezufuhr, nämlich einerseits die ‚stetig verteilte‘ Wärmezufuhr und andererseits die 
mit turbulenter Luftmassenumlagerung verbundene ‚„‚Konvektion“, gleichzeitig nebeneinander 
behandelt werden. Die stetig verteilte Wärmezufuhr pro Massen- und Zeiteinheit (= 5) tritt 
als Störfunktion in der Zustandsgleichung auf?). Hingegen macht sich die Konvektion als 
Massenergiebigkeit pro Massen- und Zeiteinheit (= e) in der Kontinuitätsgleichung bemerkbar®). 
Wir wollen uns hier auf den Spezialfall'beschränken, daß die Wärmezufuhr in einem endlichen 
Bereich stattfindet, und daß ferner die stetig verteilte Wärmezufuhr am Erdboden verschwindet. 

Um das Problem in mathematischer Hinsicht zu vereinfachen, ersetzen wir die Erdoberfläche 
durch ihre Tangentialebene an einen Punkt im interessierenden Bereich und legen im ungestörten 
Grundzustand eine isotherme, reibungslose Atmosphäre von einheitlicher Horizontalgeschwin- 

‚digkeit U zugrunde. Durch die zu U parallele x-Richtung und die dazu senkrechte y-Richtung 
in der Tangentialebene sowie die z-Richtung zum Zenit sei ein rechtshändiges kartesisches Koor- 
dinatensystem festgelegt, das sich gegenüber einem Inertialsystem mit der Winkelgeschwin- 


digkeit 4 = osin p um die z-Achse dreht (o=Winkelgeschwindigkeit der Erddrehung, p=geo- 


graphische Breite). 
Für den Ansatz der Differentialgleichungen soll die Methode der ‚‚mittleren“ Zustands- 
änderungen?) Verwendung finden. Sie besteht darin, daß man die Trägheitskräfte relativ zur 


!) H.Stümke: Zustandsänderungen durch großräumige Wärmezufuhr im Innern iner anfänglich 
ruhenden Atmosphäre. Erscheint in Gar Met. Rundschau. n : anfänglic 


’) H.Stümke: Theorie der Weiterentwicklung von‘ stabilen atmosphärischen Druckfeldern. Met. 


Rundschau (1950) 241. 


®) H.Stümke: Rationelle Behandlung der Konvektion bei d - 
entwicklung. Erscheint in der Met. Ferien, een ze ser 


az 


in erster Näherung (Index 1“): 


1 ME NEN 
Ri“ a 
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rotierenden Erde durch geeigrrete Führungskräfte ersetzt, die den Übergang von einer momen- 
tanen Gleichgewichtslage in die nächste bewerkstelligen. Schließlich wird die Voraussetzung 


kleiner Störungen erlauben, die Gleichungen bezüglich der Größen des Grundzustandes zu 
linearisieren. 


2. Differentialgleichungen, Formulierung der Randwertaufgabe 


Unter Beachtung der vorstehenden Annahmen verwenden wir folgende Bezeichnungen: 


BFRSTS U ee Druck, Dichte, Temperatur und Geschwindigkeit des Grund- 
zustandes. 


Pt»; P+o, T+49, U-+u= entsprechende Größen nach Eintritt der Störung, zu denen 
noch die Geschwindigkeiten v» und w in Y- bzw. z-Richtung 
treten. 


Für den ungestörten stationären Gleichgewichtszustand gelten die Gleichungen: 


Fre), Fe EL: 


gb’ P R’ RER: 
P,=—gP; P,=-1UP 


Dabei bedeutet P, den konstanten Bodendruck längs der x-Achse und b eine charakteristische 
Länge für die vertikale Schichtung. 


Betrachtet man nun das System kurze Zeit nach Beginn einer kleinen Störung (t>0), 
so erhält man nach Linearisieren die folgenden strengen Gleichungen: 


Bewegungsgleichungen: 


ABU ER =, —- PW US) rn (3), 
PBe-— u u H Un) nr or (4). 
Zustandsgleichung: 
pP b 
p+ Up +vP,+uwB— 5 la+U@+vP,+wPJ=(#—1)Pj +). 
Kontinuitätsgleichung: | 
6 TUT, + Pu. +(Po,+(Pw,=Pe :...... 0... (6). 


Die Störungen sollen nun so langsam erfolgen, daß die Trägheitskräfte relativ klein bleiben, 
also annähernd geostrophische Verhältnisse herrschen. Die Bewegungsgleichungen lauten dann 


ARu,— Pıa Ve Ki Se el A EN Wiesn Me WER ER L, % (2)1> 
A [eı U -+ Pu,] = — Piy ae a a ee Fee ai ee Fun ie S)1 
901 Ze Piz ee NE ER N ee NE NSÄR FE > Ave (4),- 


Um zu erreichen, daß diese Gleichungen näherungsweise in jedem Augenblick t>0 gelten, 
fordern wir, daß sie unter entsprechenden Voräussetzungen bezüglich Differenzierbarkeit auch 
nach zeitlicher Differentiation näherungsweise gültig sind, daß also die Beziehungen bestehen: 


APu:= ERS Re De ER ee u ı 
Ale: U + Pus]=— Pity a N an a eat, Te, Maid, /ı5 A RE (S)ır 
SR 901: —— ——Pp1 tz u Pa ea A nt denn ee en © (A): 


Setzt man diese Werte als Korrekturen in. die bislang vernachlässigten Relativbeschleunigungen 
der Gl.(2) und (3) ein?), so erhält man als zweite Näherungen (Index „2“): 


umso u ++ PO. wre e IS a RS EHER (2) 


4 [a + Pul=— pw —PUos 1 TE EEE (3)2- 
Aus den GI. (2)2> (3)2; (Mus (5), (6) findet man beim Grenzübergang {> 0 mit der Festsetzung, 


daß der Druck mit seinen räumlichen Ableitungen jedenfalls eine stetige Zeitfunktion ist, unter 


2 a Bde 
70 Stümke, Zur Berechnung der Drucktendenz bei Wärmezufuhr 24,3 a re: 
Fortlassung der Indizes die Gleichungen: 
.Pv= — Ua — E (= P«,) N a (7); 
een 2. (= — Po) ee ae ee (8), 
ga=—Me ann 2 gen a ne Se a a ei (9), 
PB : 
p+vP,+wP,—"5la+vP,twP]=@—1)Pj . (10) 
++ Puw,+ (Po), + (Pu), = Ps En ee u 5 (11). 


Diese Differentialgleichungen der ‚‚mittleren‘ Zustandsänderungen bilden den Ausgangspunkt 
der weiteren Rechnung. Man erkennt, daß sie im allgemeinen ein zusätzliches Geschwindigkeits- 
feld w, v, w definieren, welches sich nicht, wie man natürlicherweise erwartet, allmählich aufbaut, 
sondern unmittelbar nach dem Einsetzen der Störung dem geostrophischen Wind überlagert ist; 
ein kleiner Nachteil des Verfahrens, den man nach dem Streichen der wirklichen Beschleuni- 
gungen in Kauf nehmen muß. Die Größen «, v, w und o, sind durch den Wert von p, bereits 
eindeutig bestimmt. 'Eliminiert man sie aus diesen Gleichungen, so erhält man nach elementarer 
Zwischenrechnung für 9, die partielle Differentialgleichung zweiter Ordnung: 


. — 1) U? —1Ub 
A + | 2 fun + bBe 
a ie xgb AR 
Er ....(12). 
if | 
x 02 


- - Wir nehmen jetzt einen ebenen Erdboden an und schreiben demgemäß als Randbedingung das 
Verschwinden der Vertikalgeschwindigkeit w am Erdboden (z=0) vor. Gemäß Voraussetzung 
beschränken wir uns ferner auf den Spezialfall, daß am Erdboden keine stetige Wärmezufuhr 
stattfindet; die Größe 7 verschwindet daher am Boden, und wir erhalten aus Gl.(10) mit Be- 
nutzung der Gl.(1) und (7) die Randbedingung 


j x» —1) U? »—1) U 
ur mt. [1 —Piy he UM) Sure (13), 
in der bemerkenswerterweise eine Ableitung parallel zur Randfläche (=, ,) auftritt. Wegen der 
Unendlichkeit des Lösungsbereiches müssen wir außerdem für das Unendliche des Halbraums 
Regularitätsbedingungen vorgeben. In Anbetracht der Voraussetzung, daß die Störungen ganz 
im Endlichen liegen, fordern wir 


im ET SE NE Re (14a) 


Var+y+2? —o 


und zwar soll p,, damit die relative Kleinheit der Störungen sichergestellt ist, stärker als der 
ungestörte Druck verschwinden, also außerdem 


PETE RR 
lim In..e5\ y Hs TE (14b) 
Veryte- o 


gelten. Insgesamt ergibt sich somit die Aufgabe, die inhomogene Ditfgl. (12) unter der homogenen 
Randbedingung (13) sowie unter den Regularitätsbedingungen (14a) und (14b) zu lösen. 


3. Transformationen 


Abgesehen von den Größen g, , R und x, die man als universelle Konstanten des Problems 
ansehen kann, treten als stärker veränderliche Parameter die geographische Breite (in 


A=2osingy), die Temperatur T (in = sowie die Grundgeschwindigkeit U auf. Um 


diese Parameter möglichst weitgehend zu eliminieren und zugleich den Charakter der Diffgl. (12) 
festzustellen, führen wir zunächst folgende Transformationen ein: 


a: Teer al 


pı (2, y2)=W(X,Y,Z) | 


EN Wu 
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Die Diffgl. (12) nimmt dabei die Form an: 


DE Un day, 92.0 m BARS, 

ae Fre, (75,2): 2299) 

wobei I (X, Y, Z) das inhomogene Glied in den neuen Koordinaten ist. Aus dieser Form ersieht 
‚man leicht, daß die Diffgl. stets vom elliptischen Typ ist, da Be 1 ist j - 

Regularitätsbedingungen lauten: = Yi+p end 

lim [94 97% DO Near da 

lim a N en Re ET a ES j 

VERTIIZE>o Br 

a re ee ee Ne (14b). 


Vx+Y°+Z' +0 
Als Parameter tritt jetzt nur noch die Grundgeschwindigkeit U im Ausdruck u auf, der das 


/ #—1-fache der Machschen Zahl der Grundströmung bedeutet. Das Verhältnis der charakte- 
ristischen Längen 1: b ist etwa 200: 1t). 


Um die gemischte 2. Ableitung sowie die 1.Ableitung im homogenen Teil der Diffgl. (12)' 
zu beseitigen, verwenden wir weiter die Transformationen: 


Y=-Yii+w@+uZ, 


Bl | 
N A RE RE (16). 
Baorrer. org n 
Die Diffgl. für © lautet: | 
| | 
Der + Op +9 ZB-AKNDN.....2... (12)", 


wobei Q (X,9, Z) wieder das inhomogene Glied bedeutet. Die Rand- bzw. Regularitätsbedin- 
gungen nehmen die Form an: 


Kur An = mite A (13) 
z-0 22/1 + u2 
Zt (g4u9) 
lim * . eVıru Bl) En en De RR E (14a), 
VI®+J?+Z?>o \ 
ee I) 
lim Deiirn a N URRER EEE SA TSOHERTE (14b)". 


Vx°+Y+z' +» 
Der Parameter u tritt jetzt in der Diffgl. (12)’’ überhaupt nicht mehr auf und in der Randbedin- 


gung (13) nurin der Form Yl-+u?. Nun kann man unter den praktisch vorkommenden meteoro- 
logischen Verhältnissen die Machsche Zahl der Grundströmung und damit u als klein gegen 1 


u2 
voraussetzen, etwa Ma< 0,2, u< 0,13. Man kann daher schreiben YI+w=1-+ 2 und dem- 


2 
zufolge - als klein von zweiter Ordnung gegen 1 vernachlässigen. Das Problem geht damit in 


eine Form über, die für meteorologische Anwendungen in erster Näherung parameterfrei ist. 
Jedoch wollen wir im folgenden von dieser Vernachlässigung keinen Gebrauch machen, sondern 
mit den strengen Gleichungen weiterrechnen. Zu beachten ist, daß sich die Randbedingung (13)” 
gegenüber den vorangehenden Formen (13) bzw. (13) insofern vereinfacht hat, als jetzt keine 
Ableitung parallel zur Randfläche mehr auftritt. - 

Für die folgenden Betrachtungen ist es sehr instruktiv, die Lage der Ebenen 9 = const 
im physikalischen x, y, z-Raum aufzusuchen. Man findet 


Beispielsweise wird mit «—=1,4, A=10%/sec, b=8000 m, U=20 m/s, also 1= 1500 000 m, 


agb = V1,4-10-8000 = 335 m/sec,. Ma—0,06, 1—0,038 die Steigung 


bi+tm_ 8000 


et . T m 1500000.:0,0388 u? 


mit 
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man erhält also Ebenen, die nach negativen y-Werten, d.h. in Richtung zum höheren Boden- 
druck, mit der Steigung 0,14 flach ansteigen. Für U=0sind sie natürlich vertikal. Die Steigung 


erreicht theoretisch ein Minimum für u=1, jedoch hat diese Folgerung keine praktische Be- 


deutung, da so hohe Grundgeschwindigkeiten nicht vorkommen. Vergleichsweise lautet die 
Gleichung der Isobarenflächen: 


= 5 Yy + constt = — Yy Iu+ const. 


Ihre Steigung (= E „) beträgt bei den obigen Zahlenwerten nur 0,0002; die Isobarenflächen sind 
also unter normalen Verhältnissen etwa 1000 mal flacher als die Ebenen 9) = const. Für u—1 
würden sich die Steigungen der beiden Ebenscharen nur um den Faktor 2 unterscheiden und bei 
höheren Machschen Zahlen sogar asymptotisch zusammenrücken; jedoch haben diese Folge- 
rungen wieder nur noch theoretisches Interesse. 


4. Allgemeine Lösung der Randwertaufgabe 
a) Bestimmung einer Greenschen Funktion nach zwei verschiedenen Methoden 
Für die Lösung der Randwertaufgabe verwenden wir die Methode der Greenschen Funktion. 
Als solche soll jede Funktion @ der 6 Variablen X, 9, Z,€,n,& bezeichnet werden, welche die 
folgenden Eigenschaften besitzt: 

1. @ist mit Ausnahme einer kleinen Umgebung des Punktes X—=£, J=n, Z={ im oberen 
Halbraum Z>0, £>0 absolut beschränkt und gleichmäßig stetig, hat dort stetige Ableitungen 
bis zur 2. Ordnung und genügt in bezug auf die Variablen X, 9, Z der homogenen Diffgl. (12). 

2. @ wird im Punkte X=&,9=n, Z=Z wie 1/R(£) unendlich; 

RKI=-IR PM MHZ— 3. 

3. @ genügt der Randbedingung (13) sowie den Regularitätsbedingungen (14a)’’ und (14b)”. 

Unter Verwendung der Spiegelungsmethode, also im wesentlichen durch Probieren, 
findet man eine Greensche Funktion G; in der rechentechnisch brauchbaren Form?)®): 


a en Re SEHE 
AX-859—nZ De ee (17) 
1 9 N» ’ R(£) R(—t) J R(&) 1 . 
In G, tritt außer der Quellösung i 
GKX—EY—n Z, = DE (18), 


deren Singularität im Integrationsgebiet liegt, und ihrem Spiegelbild 6,(X—&, 9—n, Z, —{£) 
noch eine Quellstrecke auf, die auf der Geraden & = const, n = const liegt und von — w bis —& 


. * - 2 
reicht. Im physikalischen Raum besitzt diese Gerade die Steigung S —e gegenüber ihrer 


Projektion auf die Erdoberfläche, wie bereits in Abschnitt 3 gezeigt wurde. 

Auf mehr systematischem Wege kann man eine Greensche Funktion @, durch Super- 
position von Exponentiallösungen konstruieren. Man findet so, wenn man zur Quell- 
lösung 6, (X — 5, 9—n, Z, £) eine geeignete additive, für Z> 0 regulär-analytische Funktion 
K(X— 8,9 —n, Z, {) bestimmt, als zweite Form einer Greenschen Funktion: 


A -5,9-n2)-K-8,9—, ZI) + KX—8,9—n, Z, 0) 


KRX-&9—-,2,9= 77 W (a, ß, &) ei H On +rZga dB 
wobei Si | 
EangerS 
r=—Vi+e aa (9) 

1 
W(a,ß,d)= An(4+4,) If e-tar—ißs Vers, d) dr ds 
und 2 
7 (ns )= — im [420 +] 


*) H.Stümke: Rotationss 
Z. Geophysik (1940) 8. 127. 


ymmetrische Gleichgewichtsstörungen in einer isothermen- Atmosphäre, 


u 
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ist. Dain G,ein vierfaches Integral zwischen den Grenzen — oo, + oo vorkommt, ist diese Form 
für numerische Rechnungen unbrauchbar. 


b) Konstruktion der Lösung 


Nach diesen Vorbereitungen erhält man eine Lösung des Randwertproblems, wenn @ 
irgendeine Greensche Funktion ist, durch das Integral 


0899-7, ||| e&-5 9-2 906% 9ama ..... . 
>20 


wobei die Integration über den oberen Halbraum zu erstrecken ist. Vom inhomogenen Glied 
Q(X,9,Z) sei dabei vorausgesetzt, daß es nur endlich viele Sprungflächen besitzt, im übrigen 
aber stetig ist. ® wird dann mit seinen Ableitungen erster Ordnung überall stetig, also nach 
GI1.(7), (8), (10) auch die Geschwindigkeiten u, v, w. Die kinematische Grenzflächenbedingung 
ist damit überall erfüllt. 

Bezüglich eines Zahlenbeispiels sei auf den Bericht!) verwiesen, dessen Zahlenangaben und 
Diagramme man unschwer für das hier vorliegende Problem umrechnen kann. 


c) Eindeutigkeitsbeweis 


Um die Eindeutigkeit der Greenschen Funktion und damit der Lösung der Randwert- 
aufgabe zu zeigen, genügt der Nachweis, daß für Z>0 eine Funktion Z der 6 Variablen X, 9), Z, 
&,n, durch die folgenden Bedingungen eindeutig bestimmt ist: 

1. L genügt im oberen Halbraum bezüglich der Variablen X, 9), Z der homogenen Diffgl. 
(12). 


: 2. List im oberen Halbraum regulär-analytisch. i 
3. L genügt der Randbedingung 
lim [L- A+ Lo)]= —Uim [4:0 (X— 5 9—n, 2,467] mit >00... (al) 
z-0 = 


sowie den Regularitätsbedingungen (14a)’ und (14b)”. 


Vermöge dieses Satzes ist dann nämlich auch die Greensche Funktion @ eindeutig durch 
die Gleichung 


bestimmt. 

Die Existenz einer Funktion Z, die den vorstehenden Bedingungen genügt, geht bereits 
aus der Angabe ihres Wertverlaufs G6,—@, bzw. @,—6, in den Gl.(17) bzw. (19) hervor. Ange- 
nommen nun, die Eindeutigkeit von Z wäre nicht vorhanden, so müßte es auch eine nicht identisch 
verschwindende regulär-analytische Funktion M geben, die der homogenen Randbedingung, 


RER (23) 
z>I 


genügt und im übrigen dieselben Eigenschaften wie L besitzt. Wir bilden aus M und seiner 
Ableitung nach Z die Funktion 


ee A ne (24) 
N genügt der Randbedingung 
n lim N | 0 Dr ne ES ne Ren Zee EN are (25) 
2-0 


und besitzt im übrigen die gleichen Eigenschaften wie M bzw. L; insbesondere folgt aus den beiden 
Regularitätsbedingungen das Verschwinden von N im Unendlichen. Dies Problem hat aber 


Er: % en 5 
bekanntlich wegen des negativen Vorzeichens von 4 ® in Gl.(12)” nur die eine Lösung 


EEE (26) 

/ (vgl. Courant-Hilbert II). Folglich muß M allgemein für Z>0 der Dittgl. 
| = a er Basen 
genügen. Die allgemeine Lösung dieser Diffgl. ist 5 | ; : 
be ni Me KNEND:.:: EEE (28) 


= 


5 "wobei F im ganzen ‘oberen Halbraum regulär-analytisch ist und nach unserer Annahme nicht 
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identisch verschwindet. Für Zoo verschwindet M wie e 42, soll aber andererseits gemäß 


—Z 
der Regularitätsbedingung (14b)” stärker als e?Vit“* verschwinden. Es müßte also 
A> -———— oder N d.h. ->x 
2Y1-+u2 Ar 2 


sein, was den physikalischen Voraussetzungen widerspricht. Folglich muß F und damit M 
identisch verschwinden, womit die Eindeutigkeit bewiesen ist. 

Aus dem Eindeutigkeitsbeweis ergibt sich die Gleichheit der Funktionen @, und G,, nach 
G1.(17) bzw. (19), ein Satz, der auf direktem Wege nur schwer nachzuweisen sein dürfte. 

Alle vorstehenden Ausführungen zur Lösung des Randwertproblems lassen sich ohne 
Schwierigkeit auf das ebene Problem (Ebene senkrecht zur Grundströmung) übertragen. Die 
Lösung des Problems mit einer inhnomogenen Randbedingung (Wärmezufuhr am Boden) liegt 
ebenfalls vor, jedoch ist die Berechnung eines Zahlenbeispiels analog demjenigen in Bericht!) 
z.Z. noch nicht abgeschlossen, weshalb auf die nur formale Darstellung der Lösung hier verzichtet 
wurde. 


5. Vergleich der vorliegenden Rechnung mit einem Ansatz vonL. Prandtl 


Es sei an dieser Stelle eine kurze Bemerkung über das sonstige Auftreten einer Differential- 
gleichung vom Typ der Gl.(12) in der theoretischen Meteorologie gestattet. Die Diffgl.(12) 
stimmt im Falle U=0 und bei unterdrückter &- (oder y-)Komponente formal mit einer von 
L. Prandtl?) abgeleiteten Diffgl. für diejenigen endlichen Druckstörungen überein, die beim 
Übergang vom ungestörten zu einem gestörten Gleichgewichtszustand auftreten. Die von 
Prandtlangewandte Betrachtung läßt sich bei Unabhängigkeit aller Größen von der x-Richtung . 
(ebenes Problem) auch für den Fall U+0 durchführen. Dabei treten an Stelle der oben ange- 
schriebenen Tendenzen 7;,, 0, kleine, aber endliche Störungen des Drucks bzw. der Dichte und an 
Stelle der Geschwindigkeiten v, w kleine Verschiebungen. Mit diesen Modifikationen erhält man 
Beziehungen von derselben Form wie die Gl.(7), (9), (10), (11). Lediglich Gl.(8), deren beide 
Seiten im vorliegenden Falle identisch verschwinden würden, wäre beim Prandtlschen Problem 
durch die der Gl. (3), entsprechende Aussage zu ersetzen, daß auch im gestörten Zustand zwischen 
der Corioliskraft und dem Druckgradienten stationäres Gleichgewicht besteht. Da diese Gleichung 
aber beim ebenen Problem für die Ableitung der resultierenden Diffgl. nicht gebraucht wird, 
erhält man für endliche Druckstörungen ebenfalls eine Beziehung, die formal mit der obigen 
Gl.(12) für die.Drucktendenz übereinstimmt. Der entsprechende Ansatz nach Prandtl führt 
also auf eine Diffgl., die sich bei der hier verwendeten Methode als Spezialfall ergibt. 


u 


6. Energiesätze 


Wir wollen zunächst die Zustandsgleichung in Form einer Energiebilanz für die 
Masseneinheit schreiben. Unter Verwendung der Temperatur als Zustandsgröße an Stelle 
der Dichte lautet G1.(5) für 1>0 


+ U) SI + Up. HoPtwPl=i 2.2 ..2.. O% 


Für 0 wird daraus j 
1 . . - y 
9 — plr+ ei FePler VRRIE AED ENIUNE 
Diese Gleichung läßt sich in Verbindung mit den Gl.(1) und (7) in der Form schreiben: 


i+5=ohtm Utug a Re BEER (29). 


in der man sofort die Aussage, erkennt, daß bei jedem individuellen Teilchen die Wärmezufuhr 
und lokale Druckerhöhung pro Zeiteinheit dazu dienen, seine Enthalpie sowie seine kinetische 
und potentielle Energie zu erhöhen. Es sei hervorgehoben, daß hier gegenüber dem Fall der 
anfangs ruhenden Atmosphäre (U—=0)!) als neues Glied der meteorologisch wichtige Ausdruck 
u, U für die Änderung der kinetischen Energie der Masseneinheit auftritt. 

_ , Umdie Energiebilanz für ein ortsfest vorgegebenes Volumen V mit der Ober- 
fläche‘ F aufzustellen, müssen wir zum Ausdruck bringen, daß in der Zeiteinheit die Änderung 
der Gesamtenergie innerhalb des Volumens V gleich ist der Energiezufuhr von außen, vermindert 


"um die Energiemengen, welche durch den Massentransport und die dabei geleistete Arbeit gegen 


} 5) L.Prandt]l: Beiträge zur Mechanik der Atmosphäre (Memojire 6Al iati 
gique de ? U.G.G.I. Edimbourg — Septembre 1936). Dans 1089. RE Te 


re u RE 


Pe NV VE 00 
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n den äußeren .Drbeck verlorengehen. Die Energiezufuhr von außen ist!)®) 
T flfferrlenrr anti ro)der 
_ Bedeutet ferner S die Gesamtenergie pro Volumeinheit, 2 


ER Ba EN EEE (30) 
den Llisulssekter, so ist — an. | | SdV die Änderung der Gesamtenergie und 


er fs +P)v- df die Energie; die durch den ae or und die nach außen geleistete Arbeit 


En Hlorenbeht, ‚wobei die Normale von df nach außen gerichtet sein‘ "soll. Insgesamt muß also 
die nt bestehen: 


alllse- BELUZIE 


+3 P o).] av / (S4 P)v.di (1) 


er 

E+° oder in Differentialform: 

% : RS = -Pj- sn (‚Pr es > P U): — div US + Pte an 
E-; Zum‘ Beweis dieser Gleichung Eee wir, daß die nenn zur Zeit ?>0, wenn man 
Bi = En we Baer (U-+u)2 gleich streicht, j : 

Er % 1 | 
, S=eP+ HP ++ ierawin TER) 
ER und ihre zeitliche Anderung . ir / 


= 2094 p, et 5 +++) U HU Hr ” 3 (33) 
Beim Grenzübergang a5 und bei ea Linearisieren Wird’ - 5 Ze 
R $ a Auer Er 5 Br Beet rt Pum B ER ey 


nn IE S=matn; u ee ET EOS Re 


ion der Gl. 2) aber den alien oberen Haibrantn gewinnt ı man en ER 
eriite Energie gleich der Summe der Änderungen voninnerer, potentieller EA 
e a2 N ist, ‚Es sei noch darauf en daß die Es der 


ji 
5 
Az 
& 
gas 
= 
Sr 
nn 
Es%) 
oa 
—— 
ur 
=-2 
eo) 
en 
u 
'® 
un 
=: 
= 
+ 
& 
Br 
D|- 
> 
ee 
»D 
o 
B 
(un 
=% 
2== 
per) 
tn 
fo 
[g>} 
ler } 
2. 
3 


ng En 


m KIN mern a ea ri 
ng. ES andw wertaufgabe N 1 Drucktendenz gibt nur ee erste 
T I die r 


en en Idemen Zeit At und Addition 
en am Ende dee ARE N 


Ing op ein each. * Beispiel i 
x ie ne a 


; . . ‚ 2. . Math. Mech. 
76 Collatz, Einschließungssätze bei Iteration und Relaxation pa. 92 Nr. 2/3 Febr./März 1058 


Einschließungssätze bei Iteration und Relaxation 


Von L. Oollatz in Hannover 


Wenn bei Anwendung eines Iterationsverfahrens auf ein Gleichungssysiem für n Unbekannte sich bei 
einem Iterationsschritt alle Unbekannten in demselben Sinne ändern, ergibt sich die Frage, ob bei den weiteren 
Iterationsschritten die Unbekannten sich stets wieder in demselben Sinne ändern und ob man auf diese 
Weise von unten und von oben her die Lösung eingabeln kann. Das gilt keineswegs allgemein, aber in weiten, 
gerade für die Anwendung auf Randwertaufgaben (bei gewöhnlichen und partiellen Differentialgleichungen) 
wichtigen Klassen von Gleichungssystemen. Auch nichtlineare Systeme werden betrachtet. Besonders 
herausgegriffen werden die Iterationsverfahren in Gesamt- und in Einzelschritten und die Relaxation. 


The author considers the case that a system of equations with n unknowns is solved by iteration, and 
that, during one step, the correction terms of all unknowns have the same sign. It is asked, if, during further 
steps, the alteration of the unknowns proceeds steadily in the same sense, and if it is possible, by this way, 
to include the solution between two limits. That is not generally true but holds for wide classes of systems 
that are fortunately importänt for applications to boundary problems (connected with ordinary and partial 
differential equations). Non-linear systems are considered, too. T’he methods of total and successive steps 
and the relawation method are discussed more intensively. ” 


L’auleur examine le cas qu’un syst&me d’&qualions avec n inconnus est r&solu par itöralion et que, a un 
pas, lous les inconnus changent dans le m&me sens. Il se pose la question, si, aux pas suivanis, le change- 
ment des inconnus se fait constamment dans le m&me sens el si de cette maniere on peut resserrer la solution 
entre deux limites. Cela n’est pas la rögle generale, mais se fait dansdesclasses &lendues de syst&mes d’&qualions 
importantes jusiement pour l’application & des probl&mes de valeurs marginales (en connection avec des 
Equalions aux derivees ordinaires .et parlielles). Aussi les systömes non-lin£atres sont consider&s. Les 
methodes d’it£ration en pas totaux et successifs et la relaxation sont traitdes sp£cialement. 


Ecım, pH IpPHMeHeHHH HTEeP&IHOHHOTO CIIOCOD& K CHCTEMe YpaBHeHHÄ C n HeNsBecTHEIMH, 
Bce HeH3BeCTHEIE H3MEHAMTCA IIPH OAHOM Mare HTepaluH B O]HOM U TOM :Ke CMBICIe, TO 
BO3HHKAET BOIIPOC, H3MEHAIOTCH IH IIPH AAIBHeÄINIHX MATAX HeH3BeCTHBIe H BIIPeNb — HOCTOAHHO 
B OAHOM MH TOM :ke CMBICHe. — H_MOKHO JIH, TAKHM OÖ0PASOM, 3&KJIHOYHTb PelmeHue Meikıy 
BEPXHHM H HH:KHHM IIpenesaMH. JTO CIHPABeAIHBO OTHIOAb He BO BceX CIIyYaaXx, HO MOKkeT 
ÖbITb B IIHPOKOH Mepe HCIHONB30BAHO MIA KIIACCOB CHCTEM YPABHeHHÄ, BAKHEIX UI IIPHMe- 
HeHUA K KPAeBEIM 3ala4aM (CBASAHHBIM C IIPOCTEIMH H C YACTHBIMH MAbhepeHmHalnbHBIMH 
ypabHenusmn). PaccMaTpHBamTcH TAK>Ke HeJIHHeÜHBIe CHCTEMEI. OCo6oe BHHMAHHE YNeustercH 
HTePaAUHOHHEIM CIIOCODAM OÖINHX H IIOCJIENOBATEJIBHBIX IIATOB, A TAKIKe METOAY PeIaKCAIHH. 


1. Fragestellung 


Es seien für n Unbekannte &,,..., z, (lineare oder nichtlineare) Gleichungen 
ICE REN a et Fa An he ER es] 


vorgelegt; es existiere ein Lösungssystem &, &, .. ., &,, und es werde ein Iterationsverfahren. 


benutzt, nach welchem aus einem Näherungssystem x)”, ..., =) ein neues System z’*”,..., 


2%.+" perechnet werde (läßt man hierbei zu, daß bei diesem Übergang sich nur eine Unbekannte x, 


. geändert hat, so kann man hier auch die verschiedenen Arten von Relaxation einbeziehen, bei 


denen nicht notwendig die %,, ..., x, in festem Zyklus periodisch iteriert werden); beim prak- 


_ tischen Rechnen kommt es dann oft vor, daß sich bei einem Schritt alle Unbekannten in gleichem 


Sinn geändert haben, daß also z.B. alle Änderungen e% nicht negativ sind: 


ar) _ >0 (k=1, sn): Er (2) 


Es ergibt sich num die Frage, ob und wann man folgende Schlüsse ziehen darf: 
Aussage A: Aus ek >0 für alle k folgt & 20, d.h. die Änderungen erfolgen auch weiter- 


hin im gleichen (nichtfallenden) Sinn. 


Aussage A’: Aus ek < 0 für alle k folgt e®*” < 0. 

Aussage B: Aus ey, = 0 folgt &—aX "> 0, d.h. die Lösung liegt nicht unterhalb der letzten 
Näherung. : 

Aussage B': Aus ep <0 folgt K— a" < 0. : 

Aussage C (Einschließungsaussage): Hat man zwei Näherungssysteme x”, & mit den 
dazugehörigen Änderungen ex bzw. EX —=&%  — 5%), so folgt aus &) 20, 2 <0 fürk=1,...,n, 


w) 
k = 
daß die a und & ein Lösungssystem einschließen a <& Sa”. 
Daß diese Aussagen A, B, C keineswegs allgemein gelten, zeigen schon einfachste Beispiele; bei 
dem Gleichungssystem | 


et 
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mit 
1 0,3 0,2\ 157 1 
A — 0,5 1 — 0,1 9 de 1,4 3 T ze 1 
02 —01 1 1,1 1 
werde etwa ausgegangen vom Vektor 
1,3 
0 — | 0,9 
1 


und nach dem Iterationsverfahren in Gesamtschritten (welches in Nr. 5 nochmals erklärt ist)!) 
der iterierte Vektor 


1,05%, 
D==0,85 
0,93 


berechnet. Es sind bei 2 sämtliche Komponenten kleiner als bei 2%; trotzdem liegen die 
zweite und dritte Komponente der Lösung oberhalb der letzten Näherung. Es liegt hier ein 
Beispiel vor, bei dem das Iterationsverfahren gut konvergiert (die Hauptdiagonalelemente über- 
wiegen stark, das ‚‚Zeilensummenkriterium“ ist erfüllt), aber die Aussagen A’,B’,C gelten nicht. 
Trotzdem gibt es weite, gerade für die Anwendungen wichtige Klassen von Gleichungssystemen, 
bei denen die Aussagen A—C gelten. Im folgenden sollen einige genannt werden. 


Insbesondere ist natürlich die Eingabelungsaussage C von Wichtigkeit wegen der einfachen 
damit gegebenen Fehlerabschätzung. Man geht einmal von Näherungen aus, die sämtlich zu 
klein sind und nähert sich so von unten her dem Lösungssystem und ein zweites Mal ent- 
sprechend von oben her. 


2. Einfache hinreichende Bedingung für die Eingabelungsmöglichkeit 


Das Gleichungssystem (1) sei auf irgendeine Weise auf eine für die Iteration geeignete 
Gestalt 


= 9x (dı, er) Kell” EEE ER (3) 
gebracht; es lassen sich verschiedene Iterationsverfahren durchführen, z.B. irde man beim 


Iterationsverfahren in Gesamtschritten Folgen von Näherungen x%" nach | 
EEE EEE N k=1l,..,n:v=01..)....:() 

und bei der Iteration in Einzelschritten nach 
d—g (a: ae u at, N, BEN 2) EEE NT (5) 


berechnen. Bei einer der verschiedenen Arten von Relaxation berechnet man die Abweichun- 


gen, die „Defekte“: 


Na (6) 


und verbessert jeweils diejenige Unbekannte &;, bei welcher der Defekt den größten Betrag hat. 
Der Einfachheit halber werde in dieser Nummer das Iterationsverfahren i in Gesamtschritten 


“nach (4) zugrundegelegt. -Es mögen die gegebenen Funktionen g; in einem konvexen Bereich ® 


des 2; ... %,-Raumes, der ein Lösungssystem und die Näherungen 2%, 2.» Zenthäh, Bip- 
schitz- "Bedingungen genügen: es soll nichtnegative Konstanten K,, ; Br s0 ‚daß für beliebige 
Wertesysteme #,, %; aus. ® gilt: 


2 In) —n(2)| < Zul a IRRE (7), 
Aus (2) und (4) od | 
nn ea We (8) 

folgt dann Ä 

” e vr | < Rule ee en De. (0). 


9 R.v. ir: H. Pefiard ek- Geiringer. Z. angew. Math. Mech. 9.(1929), 8. 58—77. 
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Mit 
Er 
ew) — 3 2 N (Kr;) ee a ee ee Re ee (10), 
(v) 
In 
folgt also i 
ew+l) < Kem, ewıp) < Gy el Fr (11): 


Sind alle charakteristischen Zahlen xg von $ (d.h. alle Wurzeln von det(8—x &)= 0 
mit & als Einheitsmatrix) dem Betrage nach kleiner als 1, so folgt nach Sätzen der Matrizen- 


oo 
lehre?) die Konvergenz der 2% für »> 0, indem man ie | majorisiert. 
v 


() = f 
Sind überdies alle Ableitungen en >0in für j,k=1,..., n, so folgen die oben genann- 
j 
ten Aussagen A—C. Denn mit geeigneten Zwischenstellen & gilt mit (8) nach dem Taylor- 


schen Satz 


n 
09 
w+1) w) *k ,.(w) 2 
ek = ej du, (&# Be ee (12); 
j=1 


also ist dann mit e” >'0 auch e”*”> 0; entsprechend ist mit e” <0 auch eX*” <0. Man 
hat also dann die Möglichkeit des Eingabelns. 


3. Beispiel einer nichtlinearen Randwertaufgabe 
Es möge die Randwertaufgabe 


Y"+e=0, y(0) = yl)=0 EN EEE a (13) 
angenähert mit dem Differenzenverfahren behandelt werden?). 


"Um zunächst einen Überblick zu erhalten, werde mit der Maschenweite A— > gerechnet. 


Für den Funktionswert % 5) a erhält man die Gleichung 


a—2a--0 


0,1260 
h2 £ 


—e*—0 mit den Lösungen a = in 430 


Bei der Maschenweite A= . treten zwei Unbekannte a, b auf (Bild 1) mit den Gleichungen 
16 (b—2 e—=0; 
ae: ta; 
16 2a— 2b) +e—=0 
Wir verwenden die Iterationsvorschrift 
at) g— RN eu”) ni Es bw) 
32 ; 

(15). 


1 
bo) — ,— aM _ ep 
92 Q on 33 ed 


Bild 1 


; Das Schema auf Seite 79 zeigt die Durchführung der Relaxation; und zwar Zeile Nr.1—9 
für die Lösung mit den kleineren Werten und Zeile 10—15 für die mit den größeren Werten 
im Auszug. In Zeile Nr.1 wurde von a=0,1; b=0,13 ausgegangen und durch Anbringen von 
Änderungen versucht, die „Spanne“ |x|4-|ß| herabzudrücken. In Zeile Nr.3 sind beide Än- 
derungen a, ß positiv und in Zeile Nr.4 sind beide negativ. In diesem Falle kann man aussagen, 
daß die fragliche Lösung von diesen beiden Zeilen eingeschlossen wird. Eine gewisse Schwierig- 
keit besteht in der Abgrenzung des Bereiches ®, in dem die Lösung liegen muß. Hier kann man 
z.B. a eliminieren; dann genügt b der Gleichung f(b)=b— eo (1+ e »°)=0; man sieht so- 
fort f(0)<0, f(In2)> 0; es liegt also eine Nullstelle von b zwischen O-und In 2 und i 

i N ’ $ ’ man d S= 
kutiert leicht, daß dieselben Schranken auch für a bestehen, sofern man sich auf die Ehsung 


mit den kleineren Werten festlegt. Wir verwenden also 0 <a <In2, 0 <b<In2 als Bereich 8. 


2), Vgl. z.B. L.Collatz: Eigenwertaufgaben. Leipzig 1949, 8.311. 


®) Die Lösungen der Differentialgleich ind i i 
sichtlich nicht Aaron Gebrauch Beganbtwergehe a ee re  F a ara for 


Pe 
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Mit e®<s2, e<2 hat man von der dem Betrage nach größten charakteristischen Zahl 
der Matrix 


1-1 
16,72: 

R— i 
a 
16 


d.h. von ihrer Maximalwurzel nachzuweisen, daß sie kleiner als 1 ist. Nach dem Ein- 
schließungssatz?) genügt es, einen Vektor u mit positiven Komponenten anzugeben, so daß 
die Verhältnisse der Komponenten von Su zu denen von u kleiner als 1 bleiben. Ein 
solcher ist z.B. 


3 35 

4 64 
= a: Ku= 13)? 

16 


die Verhältnisse und 15 sind beide kleiner als 1. 


Da ferner die partiellen Ableitungen sämtlich positiv sind, gilt nach dem oben 


IR 
da or ap” 
Bewiesenen die Einschließungsaussage C, also mit den Zahlen der Zeilen Nr.3 und 4 

0,102 <a <= 0,107, 0,135 sb < 0,144. 


Natürlich kann man leicht die Genauigkeit steigern; ZeilenNr.5bis 7 zeigen die Ausdhrurs 
einer Korrekturrechnung mit den als klein angenommenen Korrekturen 6, &,_welche sich aus 
0,035 6+0,5 8 — ö= 0,00026 | 16 
Bose 0,0002. EEE, 12 
zu 
ö—= —0,000575, e—= — 0,000574 
ergeben. 3 
In Zeile 8 und 9 sind Korrekturen so angebracht, daß die Anderungen a, ß jeweils beide 
positiv, bzw. beide negativ sind. Damit hat man die Einschließung?). 


0,105446 <a < 0,105450 , 0,141443 <b < 0,141449.. 


Bei der anderen Lösung mit den größeren Werten bleiben die charakteristischen Zahlen 
der zugehörigen Matrix dem Betrage nach nicht unterhalb von 1, und man kann hier leicht 
Beispiele dafür geben, daß man nicht aus der Tatsache, daß alle Änderungen in einer Richtung 
erfolgen, schließen darf, die Lösung müsse auch in dieser Richtung liegen. Zeile Nr.11 ergibt 
beide Änderungen a, ß negativ, trotzdem ist der Wert b=3,5 zu klein, entsprechend sind in Zeile 
Nr.12 a und ß negativ und beide Werte a, b sind zu klein, Zeile Nr.13 gibt positive a, ß, trotzdem 


j Os = Sp anne 
ne a”) A] ga”) 5 »w) a®+D | pw+D ar+d _ My +b _;m | ja r + Bl 
ı 101 0,13 1,1052 [1,1388 | 0,09954 | 0,13559 | —0,00046 | 0,00559 | 0,006. 03 
2 |o,1 0,135 1,1062 [1,1445 | 0,10204 | 0,13577 0,00204 | 0,000 77 | 0,002 81 
3 [0,102 [0135 .|1,1074 11445 |0,10211 | 0,13777 0,000 11 | 0,002 77 | 0,002 88 
4 10,107  |0,144 1,1129 1,1549 | 0,106778 | 0,143.091 | —0,000 22 |—-0,000 91 | 0,00113 
5 | 0,106 10,142 1,1118 1,1526 | 0,105744 | 0,142019 | —0,00026 | 0,00002 | 0,000 28 
6 »+ö| „tel,+L11ö,,+1,15e| „+ 0,0356] „+ 0,036 
Eu +0,5€ Da 
7 | 0,105 425 0,141 426|1,111183|1,151 915 | 0,105437 | 0,141 422 0,000 012|—-0,000 004 | 0,000 016 
s | 0,105 445| 0,141 442|1,111 205| 1,151 934 | 0,105446 | 0,141 443 0,000 001| 0,000 001 | 0,000 002 
9 Io, 108 451/0,141 450| 1,111 212|1,151 943 | 0,105450 | 0,141449 | —0,000 001|—0,000 001 | 0,000 002 
10 4 20,09 54,60 2,628 4,71 —0,372 0,71 1,08 | 
Zu 3 ln 35 11,02 33.12 2,09 3,44 —0,31 0,06 0,37 
12 3 7,39 20,09 1,731 2,63 | 0,27 0,37 0,64 
x 2, 2 4 9,974 |54,60 2,3116 4,006 0,0116 0,006 0,018 
14 En „te |,+9976,,+5468€| ,, ER | a 
ee, j „IE 
15 | 2,305 | 3,984 10,024 |53,732' | 2,30525 | 3,98413 0,000 25 | 0,00013. | 0,000 38 


4) Vgl. z. B. Math. Z. 48 (1942), S. ‚221. ä Fa | 
5), Die Einschließung bezieht sich nur auf die Wa a, b, nicht auf y (4) und (5) i 


- 


% 


En 
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ist b=4 zu groß. Alle drei Zeilen 11 bis 13 geben Gegenbeispiele gegen Aussage B. In Zeile 13 
bis 15 ist eine Korrekturrechnung durchgeführt, die mit = 

0,3125 +0,5&e — 6 = — 0,0116 | ö= 0,00526 

1,1 e+ 5-2 =—0006 J e=— 0,0159 
die Spanne |a|]-+|ß] stark herabdrückt und zeigt, daß diese Korrekturrechnung auch zum 
Ziele führt, wenn das Iterationsverfahren in der gewöhnlichen Art versagen würde. 


4. Spezialfall des linearen Gleichungssystems 


Um Aussagen für lineare Gleichungssysteme zu erhalten, könnte man die Ergebnisse vom 
nichtlinearen Fall spezialisieren. Man bekommt jedoch etwas schärfere Resultate, wenn man die 
linearen Gleichungssysteme direkt in Angriff nimmt und die bei linearen Gleichungen weiter 
entwickelte Theorie benutzt. Es sei das Gleichungssystem vorgelegt 


Ar ==T 7a ee Le ee u (17). 
Mit A=98-+-€ und nichtsingulären Matrizen X, ® lautet die Iterationsvorschrift 
Bird. EHI EEE RT ee (18). 
Die Änderungsvektoren 
2 HN De a Se ee (19) 
genügen dann den Gleichungen 
RW) == % pr—) RR RY Hy) we Te Le a a (20) 
mit 
NY = — BE I ee Fr eh ee Te (21). 


Die Aussage Aist dann und nur dann bei beliebigem 5") richtig, wenn alle Elemente f;, der 
Matrix % nichtnegativ sind, wie aus d"+D —%5) unmittelbar folgt. (Wäre ein f„,<0, so wähle 
man d”>0, d” —0 für o+g, dann würd KH” = f,,% <0.) . 

Das Iterationsverfahren (18) heiße konvergent, wenn es bei beliebiger Wahl des Aus- 
gangsvektors x konvergiert. (Dafür ist bei nichtsingulären W, ® bekanntlich notwendig und 
hinreichend, daß alle charakteristischen Zahlen von % Beträge kleiner als 1 haben.) 

Dann gilt für ein konvergentes Iterationsverfahren die Aussage B, sofern f;.>0, % also 
nichtnegativ ist; es ist nämlich 


co 


Tr mE DIS EN ee (22) 
+ e=v+1 S 
mit 
6=3 %. 
e oe=0 
Da die charakteristischen Zahlen von % Beträge kleiner als 1 haben, konvergiert die 
Matrizenreihe, welche & definiert, und wegen f;, >20 hat auch & nur nichtnegative Elemente. 
Hat also d”+D nur nichtnegative (nichtpositive) Elemente, so auch — r”+D, womit auch hier 
die Möglichkeit des Eingabelns gegeben ist. R 


5. Eingabelung bei den Iterationsverfahren in Gesamt- und Einzelschritten 


Ein wichtiger Spezialfall ist der einer Matrix X mit positiven Hauptdiagonalelementen 
@;;>0 und nichtpositiven Elementen außerhalb der Hauptdiagonale a,.<O0 für j+k. Beim 
Iterationsverfahren in Gesamtschritten verwendet man als ® eine Diagonalmatrix mit den a,;: 


fa 0 Para 0) 0 da ... (un 
RE \; Ag +. 0 g— | 0 Bra 
0 Q, a 0 
Hier ist 5 ee 
1 
— 0 
Aıı 
g-1— 0 er 0 ? 
00 3 
Inn 


und 3=—B-1C offenbar nichtnegativ. 2 \ 
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Das gleiche gilt beim Iterationsverfahren in Einzelschritten. Hier verwendet man 


RN 0 0 0 
4 0 0 0 
= Qzı Ag Ag 0 ’ 
Anı Img (ng Ann 
[ 0 a Ms An 
= 0 0 Qg3 Ag 
\o 0 0 0 
Dann hat 3-! nur nichtnegative Elemente; denn die Vorzeichenverteilung ist bei ®: 
+00 || 
—+0...0 
(8) = > 
ED 


wobei an Stelle von negativen Elementen auch Nullen treten können; bei ®-! hat man das 
Schema 


205.0 0 
++0..0 
@ay=+++...0|, 
+++... + 


wobei wieder an Stelle der positiven Elemente außerhalb der Hauptdiagonale Nullen treten 
können. Die Richtigkeit der behaupteten Vorzeichen ergibt sich leicht spaltenweise durch voll- 
ständige Induktion. Somit hat auch jetzt %—=—®93-€ nur nichtnegative Elemente. 

Zusammenfassend gilt daher: Ist bei einer Matrix X mit a,;>0, a;.<0 für j+k das Itera- 
tionsverfahren in Gesamtschritten (bzw. ın Einzelschritten) konvergent, so gelten die Aussagen A bis CO; 
es besteht die oben genannte Möglichkeit des Evngabelns. 


6. Anwendung auf Randwertaufgaben bei gewöhnlichen und partiellen 
Differentialgleichungen 
Die Voraussetzungen dieses Satzes sind in vielen für die Anwendung wichtigen Fällen 
erfüllt. Wir geben einige Beispiele; 
1. Der Randwertaufgabe ; 


= (f@) +s@)y=r(@), ya=yw, YyO=W:-...-- ) 


mit stetigem g(x), r(x), stetig differenzierbarem f(x) und f(2)>0, g(#) Z0 in asz<b werden 
bei näherungsweiser Lösung mit dem gewöhnlichen Differenzenverfahren die Differenzen- 


gleichungen 
her Y,ut (+2+4;-3) 14-1379, Y=n;. : . (24), 
Yo = Ya» Y,=%» G=1, +’, n—1) 
zugeordnet. Dabei ist das Intervall <a,b> inn gleiche Teile der Länge h= 70-0) zerlegt, 


und der Wert einer Funktion f, g, ... an der Stelle »—=a+kh wie üblich durch Anhängen des 
Index % gekennzeichnet, wobei k nicht ganzzahlig zu sein braucht. Y; bezeichnet den durch - 
Lösung des linearen Gleichungssystems (24) erhaltenen Näherungswert für die exakte Lösung 
- y(z;). Die Koeffizientenmatrix des Gleichungssystems erfüllt die oben von 4 geforderten Vor- 
zeichenbedingungen, und es konvergieren sowohl Gesamtschritt- als auch Einzelschrittver- 


_ fahren®). Also bestehen hier die Aussagen A—C mit der Eingabelungsmöglichkeit. 


6) L.Collatz: Über die Konvergenzkriterien bei Iterationsverfahren für lineare Gleichungssysteme. 


Math. Z. 53 (1950), 8. 149—161. | 
- r E - 6 
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Entsprechend lassen sich die 2. und 3. Randwertaufgabe behandeln. Es seien etwa an 
Stelle der Randwerte %,, %, die Randbedingungen vorgegeben 


y(l)—Ay()=B, yO)+AYeÜ)=B:..- > 23) 
mit a<b, A,>0, A,>0. Dann verwenden wir (abweichend vom obigen) ein Gitter, bei dem die 


. . ® b—a ibt b 
Stellen a und 5 zwischen zwei Gitterpunkten liegen: Die Maschenweite Nee bleibt, aber 


die Gitterpunkte sind x =a+ (# — —) h. Den Randbedingungen entsprechen jetzt die Glei- 


chungen: , y Y Y 
en Sin, m 
oder & ; 5 i 4 nm 
A 1 {Ag FR 
(245) 5) Br (+3) nl 3): @). 


h sei bereits so klein gewählt, daß die runden Klammern positiv ausfallen. Dann sind alle Vor- | 
aussetzungen (auch das „schwache Zeilensummenkriterium“, welches die Konvergenz des 

Iterationsverfahrens sichert) erfüllt, und es gelten wieder A bis C. 

2. Bei der elliptischen Differentialgleichung 
u u ou u_- 

= 4 -— u — ——— GUEER  Inee Fle 27 
L[u] GE ee, gu=r (27) 
seien die gegebenen Funktionen a, c, d, e, g, r in einem Bereich der x-y-Ebene stetig, a, c positiv 
und g nichtnegativ. In einem rechteckigen Gitter mit den Maschenweiten A, !, den Gitterpunkten 


ah | Re 
Ye Yo+ kl 0 1,2, ...) | 
und Bezeichnung der Funktionswerte an einer Stelle z,, y, durch Anhängen der Indizes 5, k 
entspricht der Differentialgleichung (27) die gewöhnliche Differenzengleichung für die Näherungs- 
werte U;r | 
20; , 26 1 h 1 Be 
(ae age) u ln Es) Una len En 
ET | I 14 I 
| 5%) (8-1 (+58) OH 40 


Es seien A und /so klein gewählt, daß alle runden Klam- 
mern nichtnegativ sind. 

‚ Sind längs einer Randkurve /’ die Werte von % vor- 
geschrieben, so treten mit den Bezeichnungen des Bildes 2 
noch Gleichungen der Form hinzu 


($+A)u(A)—8u(O)—=hu(B) . . . . (29), | 
wobei u(B) gegeben ist. Auch jetzt sind die Vorzeichen- 
forderungen erfüllt; es konvergieren Gesamtschritt- und 
Einzelschrittverfahren, und es bestehen daher die Aussagen 
Bild 2 A—C mit der Eingabelungsmöglichkeit. Auch hier las- 
sen sich die’ Betrachtungen leicht auf allgemeinere Dif- 


- ferentialgleichungen (z.B. in mehr Dimensionen) und auf andere . Randbedingungen über- 
tragen. + 


! 


(28). 


?. Beispiel für Durchführung der Fehlerabschätzung bei der Relaxation 
Es sei für einen quadratischen Bereich O<z<s1, 0<y<si 
4du=0 .. im Innern 


und auf dem Rande u=y für <=0, «=1 und y=1, = =0für y=0 vorgegeben; uR,y) 


läßt sich in bekannter Weise deuten als stationäre Temperaturverteilung (an drei Randgeraden 
' vorgegebene Temperatur mit linearem Anstieg und längs der vierten Randgeraden Wärme- 


isolation, Bild 3). Bei Benutzung des Differenzenverfahrens mit der Maschenweite h =+ hat ° 


\ 


ärz 1952 "Col r atz, Einschließungssätze bei Iteration und Relaxation 788 < 
man 10 unbekannte Funktionswerte. Durch Rechnung mit der gröberen Maschenweite } = - E- 


und graphische Interpolation wurden die Ausgangswerte in Schema I jeweils links oben er- 
- halten, z.B. an der Stelle A der Wert 0,63; aus Symmetriegründen. ist nur die eine Hälfte des 
Quadrats gezeichnet. Die Mittelbildung aus den vier 
Nachbarwerten ergibt die Zahlen rechts oben, also 0,6275 
an der Stelle A; die Änderungen sind in Einheiten der 
- 4.Dezimale rechts unten hingeschrieben, bei A also —25. 
- ——_ Nach Art der Relaxation wurden die Korrekturen von 
Schema II an. den Funktionswerten angebracht, welche die 
Änderungen des Schemas III hervorrufen; in Schema III 
sind die von Schema I übernommenen Änderungen durch 
Unterstreichen gekennzeichnet. Bringt man die Korrek- 
turen von Schema II an den Funktionswerten von 
Schema I an, so erhält man die neuen Funktionswerte 
links oben in Schema IV, wo also der maximale Än- 
derungsbetrag von 200-10-* bereits auf 50-10-2 heruntergedrückt ist. 
Durch Fortführung der Relaxation wurden die Änderungen auf 10-10- (SchemaV) und 
dann auf 2-10 (SchemaVI) herabgedrückt. Wollte man sie noch weiter verkleinern, so müßte 
2 man eine weitere Dezimale hinzunehmen. Nun interessiert, wieviele Dezimalen gesichert sind. 
2 In Schema VI sind alle Änderungen nichtnegativ, also liegt die gesuchte Lösung der Differenzen- 
4 


- gleichungen nirgendwo unterhalb der Funktionswerte in SchemaVI. Nun wenden wir zweimal 
Blockrelaxation’) an, Schema VII und VIII, und zwar erhöhen wir alle Zahlen der Spalte «0,4 


II IV 


Br: 1m 1,00 


- ER 
en = Ta > 


E:; 0,80. | 051 8100 0,81 | 8150 1410 8100 |8115 8160| 8175 
& = 25 |2 BET: +15. 
+10 +15 
30 | 6300 [6285 6440| 6450 
= =.50. 150-410: "15 |+0 
+25 +25 . | Re: x 
+25 -.|4600 14575 4900| 110 
50 —1001 150 —75 5 |4+10 
NE Pe ra 
+25 or 3100 3150 3700139725 
4200 +50-125+900 | I-+50 +25 
122009- °.\250..;* Bee 
EiET +50 2300 [2350 320013225 


A De 
[ı+r210-4 | © ER 
2 
| 8148| 8240 
ar ae 
| 6340| 6544 
oo & 
4664, 5044 
0 BL 
3272| 3924 


2500) zu 3 je 
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m 8 Einheiten der 4.Dezimale, und dann alle Zahlen beider Spalten 20,2 und x—=0,4 noch- 
Hals um 12 Einheiten, und wir erhalten in SchemaIX durch Addition der Änderungen von 
SchemaVIII neue Änderungen, die sämtlich nichtpositiv sind. Wenn man also an den Funk- 
tionswerten von Schema VI alle Werte der linken Spalte 20,2 um 12-10 und alle Werte 
der rechten Spalte &—=0,4 um 20-10 erhöht, hat man in den Zahlen von SchemaX sicher 
obere Schranken für die Lösung der Differenzengleichungen, und diese ist somit auf ‚einfache 
Weise in Schranken eingeschlossen, also z.B. an der Stelle A: 0,6328 zu <0,6340. (Die Block- 
relaxation ermöglicht so sehr rasch eine Einschließung; man könnte natürlich auf sorgfältigere 
Weise die großen Änderungen —3, —7 in SchemaIX herabdrücken und engere Schranken 


gewinnen.) 
Eingegangen am 31. Mai 1951. 


KLEINE MITTEILUNGEN 


Eine neue Zwangsführung zum Nyströmschen 
Stieltjesplanimeter.!) 


Das Stieltjesplanimeter dient bekanntlich ?,3)4)5) 
zur Auswertung von Integralen der Form 


ta ta 
1= 5 sw gar= S’syanth, 


welche als Stieltjesintegrale bezeichnet werden. Neben 
den Spezialinstrumenten für bestimmte Arten dieser 
Integrale, z. B. Momentenplanimetern, harmonischen 
Analysatoren usw., gestattet das Stieltjesplanimeter 
eine Auswertung bei beliebigen Funktionen f(t) und 
h(t) und zwar in folgender Weise: Der Fahrstift A 
(Bild 2) wird von einer ersten Bedienungsperson auf 
der Kurve y=f(t) bzw. auf der Abszissenachse bewegt; 
eine zweite Bedienungsperson verfolgt mit einer Fahr- 


lupe die auf dem beweglichen Tisch 7’ befestigte Zeich- . 


nung der Kurve z=hft). E.Laurila) ersetzt 
die Tätigkeit der zweiten Person durch eine mecha- 
nische Zwangsführung. Er bildet die Kurve x =(1l) 
durch ein Stahlband nach, das auf dem beweglichen 
Tisch 7 befestigt wird, und setzt an Stelle der Fahr- 
lupe eine Rollenführung, in welcher das Stahlband 
gleitet. Damit die Vorrichtung ohne Nachhilfe arbei- 
tet, darf die Steigung der nachzufahrenden Kurve 
nicht größer als 70° werden. Ferner muß der Krüm- 
mungsradius der Kurve h(t) stets größer als der Rollen- 
radius bleiben, 

Die neue Zwangsführung (Bild 1, 2 u. 3) vermeidet 
die hierin liegenden Beschränkungen. Sie verwendet 
Blechschablonen & in Form der Kurven 


REIT 


[7 


Bild1. Prinzipskizze der neuen Zwangsführung mit Blechschablone 
(a = Bewegungsrichtung von E, b = Bewegungsrichtung von T) 


2) Aus dem Institut für Praktische Mathematik d 
ne N, Prof. Dr. A. Walther, wg 
"FR... illers; Mathematisch 
Berlin Naht. r sche Instrumente, (2. Aufl.) 
) W.Meyer zur Capellen; Mathematische Inst: 
3. Aufl. Leipzig (Akadem, Verlagsgesellschatt) 1949, darge 
BANN T, Nyström; Ein Instrument zur Auswertung von 
nen. Boc. Sci. Fennica. Commentat, phys.-math, 


°) E. Laurila; Über das Nyströmsche Stieltjesplanimeter, 


_ Boc. Sci. Fennica, Commentat, phys.-math, X, 7. (1939) 


z—h(t). Am Rand gleitet ein Metallstift 81, der durch 
ein Gewicht G mit dem Seilzug $, gegen den Rand 
von S gezogen wird. Bei großen a hebt man 
von Hand mit dem Seilzug S, den von @ herrührenden 


Zug auf. Diese Betätigung von 8, muß gleichzeitig mit 
dem Nachfahren der Kurve y=f(t) erfolgen, so daß 
der Stift Stständig den Rand der Schablone berührt. 
Dies kann man bei einiger Übung ohne Schwierig- 
keiten erreichen. 


Au 


Bild 2. Gesamtansicht des Stieltjesplanimeters mit der neuen 
Zwangsführung 


Damit ist es möglich, die Zwangsführung auf alle 
praktisch- vorkommenden Kurven z2=h(t) anzuwen- 
den, auch solche mit Ecken und mit nahe benach- 


barten Kurventeilen, wie sie durch Herausschneiden - 


von Polen zustandekommen. Auch lassen sich die 
Schablonen 8 leicht auswechseln und verschiebbar 
anordnen wie etwa im nachfolgenden Beispiel. Ihre 
Herstellung ist oft einfacher als die Nachbildung der 
Kurve mit einem Stahlband. , 


Anwendungsbeispiel: Auf Anregung von 
H. Witti ch (Karlsruhe) wurde mit Hille des 
Se REBetinieteee die Konvergenz der Iterationgs- 

olge 


EN 


In+ı (9) = ar 


Sla+gn (a)]ete——? da 


So 


mit 
(9) =hy1-(1- Msn und 90(9) = 0 


für verschiedene Werte von k untersucht. Der Inte- 
grand hat für «= einen Pol 1.Ordnung. Deshalb 
wurde das Integral zwischen den Grenzen 9—e und 
+8 durch Reihenentwicklung des Integranden nähe- 


Tungsweise berechnet. Dabei wurde & auf Grund einer 


Fehlerabschätzung beim ersten Iterationsschritt zu 
e=1? festgelegt. 


= 


7 
En 


z 


\ 
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Um auf die Gestalt 
wir um: 


27 
fra 
0 2n 


= [m [l—(1— 12) sin? 0]4(21n 
0 


S f(a) dh(&) zu kommen, formen 


—(1—%2)sin? O]ctg* de 


i a 
sin —.— | R 


Die nachzufahrenden Kurven sind also 


y=f(a)=iIn[1—(1— %2)sin? 0] 
und 


x =h(a)=21n sin? 


In Polnähe' tritt von a=9—1° bis a=p-+1° an 
Stelle von z=h(a) eine Parallele zur &-Achse. 

Die für die Auswertung verwendete Zwangsführung 
(Bild 3) setzt sich folgendermaßen zusammen: Die 
Schablone S aus 2 mm starkem Duraluminiumblech 
wird mit Hilfe einer Führungsleiste L in der Nut der 
auf den beweßlichen Tisch T geschraubten Platte P 


verschieblich gelagert. Die Arretierung erfolgt mit 


Hilfe der beiden Schrauben M. Durch einen besonders 


angefertigten Einsatz E ist der Stift St fest mit dem 
beweglichen Schlitten B verbunden. Da man erreichen 
will, daß sich die mittlere Achse von St exakt auf der 


meer =  reani 
F 


Bild 3. Einzelteile der neuen Zwangsführung 


Kurve r=h(a) bewegt, legt man bei der Härstellung 
der Blechschablone die Parallelkurve im Abstand des 
halben Durchmessers von St zugrunde. Um Abwei- 
chungen an den Ecken der Kurve zu vermeiden, muß 
der Durchmesser von St möglichst klein gehalten wer- 


den. Er beträgt im vorliegenden Beispiel 0,6 mm. 


Auf diese Weise war es möglich, die Konvergenz der 
vorgegebenen Iterationsfolge für k=0,9 bis 0,5 nach- 
zuweisen. Die mit der beschriebenen Zwangsführung 
erzielten Ergebnisse weichen bei Kontrollumfahrungen 


im allgemeinen nur 1 bis 2 Noniuseinheiten in den am 
Polarplanimeter abgelesenen Werten voneinander ab. 
Diese Genauigkeit kann eine zweite Bedienungsperson 
beim Nachfahren mit der Lupe nur in seltenen Fällen, 
unter erheblich größerem Zeitaufwand und mit starker 
Ermüdung erzielen. 


Darmstadt. W. Knappe. 


Iterationsverfahren zur Berechnung komplexer 
Nullstellen von Gleichungen. 


Herr Collatz hat in der Z. angew. Math. Mech. 
30 (1950) 8.97—101einIterationsverfahren zur numeri- 


‚schen Lösung algebraischer Gleichungen angegeben, 


das in vielen Fällen wesentliche Vorzüge gegenüber 
sonst üblichen Methoden besitzt. Der Zweck der fol- 

enden Mitteilung ist einerseits eine Vereinfachung des 
Ee ibusbeweises von Collatz, verbunden mit 


einer Erweiterung auf beliebige (nicht notwendig alge- 
braische) Gleichungen, andererseits eine -Verallgemei- 
nerung auf Gleichungssysteme für mehr als eine Un- 
bekannte. 

I. 


Wir beginnen mit dem Fall einer (nicht notwen- 
dig algebraischen) Gleichung 
EEE Er 
Jede solche Gleichung läßt sich mit 
o@)=2m Fi); Fe)=zm — ol) . .(1.1a) 
auch schreiben in der Gestalt 


mir =ol)- tt 2..... (1.2) 
oder 
1 
z=Ö() mit 5le)= [ple)-ar®tr, (1.3). 
Das Iterationsverfahren 
TIL RD (1.4) 


zur Berechnung einer Wurzel 2z=£ der Gl. (1.1) kon- 
vergiert unter der Voraussetzung eines hinreichend 
guten Anfangswertes 2, immer dann, wenn 


IFaO|<1 (1.5). 


Nun folgt aber aus (1.3) mittels logarithmischer Dif- 
ferentiation (Nenner + 0 vorausgesetzt) . 


ee) EN HR. 


p %| 


Für die Lösung z={Z ist nach Gl. (1.1a) o(d) = 5” 
und nach Gl. (1.3) D(d)=L. Für 2=L£ folgt aus 
Gl. (1.6) daher sofort 

In [eo 


le mr rt 


+] a 


Bei gegebenem m (und gegebenem » bzw. F) ist daher 
nicht nur Gl.(1.5) erfüllt, sondern sogar ®’(d)—=0, 
wenn nur 

2:49) 


Er (1.8) 


en 
und wenn zugleich r-#—m (und +0.) 
r=—m=— = würde bedeuten m + {X=1 _ o’(&) 


=(0, d.h. wegen Gl. (l1.1a), daß zugleich mit 

F(£)=0 auch F’(d)=0 wäre. Schließt man diesen 

Fall einer ‚„‚mehrfachen“ Wurzel (und den Fall {=0) 

aus, soist der in Gl.(1.8) angegebene Wert von r für 

das Iterationsverfahren Gl.(1.4) immer brauchbar. 

Man hat daher den folgenden 

Satz 1: Zu jeder Gl. (1.1) und jedem Exponenten m 

gibt es unter der Voraussetzung, daß 2=& 
eine einfache von 0 verschiedene Lösung be- 
deutet, immer solche Werte von r, für die 
das Iterationsverfahren Gl. (1.4) unter der 
Voraussetzung eines hinreichend guten An- 
fangswertes konvergiert!). 
Der günstigste Wert von r bei gegebenem m 
folgt dabei aus Gl. (1.8). 

Zum Beweis des Satzes 1 war es nicht notwendig, 
F=0 als algebraische Gleichung vorauszusetzen. Für 
die Brauchbarkeit des Verfahrens ist nun allerdings 
weniger die Aussage des Satzes 1 entscheidend, die 
hinreichend gute Anfangswerte voraussetzt; wichtiger 
ist in diesem Zusammenhang die für algebraische Glei- 
chungen bereits von Oollatz an Beispielen ge- 
zeigte verhältnismäßig geringe Empfindlichkeit der 
Rechnung gegenüber schlechten (von allen Lösungen & 


2) Zurgeeigneten Auswahlunter den Werten der (im allgemeinen 
mehrdeutigen) Funktion ® vgl. Collatz. 

Die von Collatz für algebraische Gleichungen behandelte 
Frage, inwieweit man. mit reellem r auskommt, ist hier beiseite 
gelassen, n 
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weit entfernten) Anfangswerten: während z.B. das 
Newtonsche Verfahren bei solchen ‚schlechten‘ An- 
fangswerten praktisch unbrauchbar ist, kann man bei 
algebraischen Gleichungen aus einem von 
allen Wurzeln £ weit entfernten Anfangswert häufig 
mit wenigen Schritten des Verfahrens (1.4), die 
sämtlichen Wurzeln approximieren. Dies kann bei 
transzendenten Gleichungen natürlich nicht im glei- 
chen Maß gelten, und zwar nicht nur deswegen, weil 
unendlich viele Wurzeln möglich sind; denn bei einer 
transzendenten Funktion F bzw. p kann im 
allgemeinen die Multiplikation mit einer Potenz ?’, 
wie sie in Gl.(1.2) geschehen ist, überhaupt kein so 
kräftiges Mittel sein, um über weite Bereiche der 
z-Ebene hinweg soweit ausgleichend zu wirken, daß 
®’ 20 würde, 

Bei Beschränkung auf geeignete Bereiche der 
z-Ebene wirken sich trotzdem die Vorteile des Ver- 
fahrens auch für transzendente Gleichungen aus. Dies 
gilt besonders dann, wenn man (vgl. Collatz) die 
gegebene Gleichung statt in der Gestalt der Gl.(1.2) 
noch allgemeiner in der Gestalt 


(Fo) =yle) rn Base. (1.9) 


schreibt und Gl.(1.9) als Ausgangspunkt für ein Ite- 
rationsverfahren verwendet (a=0 gibt mit veränder- 
ten Bezeichnungen dasselbe wie Gl.(1.2)). Für die 
Praxis kommt es natürlich entscheidend darauf an, 
die gegebene Gleichung mit einer möglichst einfachen 
Potenz von 2 zu multiplizieren und damit zu erreichen, 
daß (1.5) bzw. die entsprechende B:dingung im Fall 
ei für möglichst viele Lösungen gleichzeitig er- 
üllt ist. 


17; 
Wir betrachten als Beispiel die Gleichung 
F=2—-A422—3e=0..... (2.1). 


GesuchtseienLösungenmit negativem Real- 
teil. 
Wir multiplizieren mit 2 und erhalten 


Hide . 2.2... (2.2) 


(zB —2i = —4+3ze ... .(23). 


Dies hat die Gestalt der G1.(1.9). Das zugehörige 
Iterationsverfahren lautet nun 


2y+1 = ) 2i-+ de Izye” = vy(2y) (2.4). 


Wir gehen aus von ,=—1 (für diesen Wert ist 
y’=0) und erhalten bei Beschränkung auf solche z, 
deren Realteil nicht positiv ist, die Möglichkeiten 


oder 
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Sielassensich mit 
pla;y) = amy" — Ka; y) ; ... (8.18) 
y(z;y) = arys — la; y) 
auch schreiben in der Gestalt 
amtp. yitv=olay) an y= ee (3.2) 
arte. ysto = yla;y)- we yo=Ple;y) 
oder auch (x y=+ 0 vorausgesetzt) 
(m + u)-loge + (n + v)log y=log D(x;y) | (3.3). 
(r + o)-loge + (+ o)logy=log Plz;y) } 


Das zugehörige Iterationsverfahren (rechts z= x, und 
y= ys, links = 0y+1und y= yy+1 eingesetzt) erlaubt 


bei gegebenen x, und y» die Auflösung nach zy+1 und . 


Yv+1, wenn 
D- m+tyu, n+v 


er INSERIEREN: 


Man hat dann letzten Endes ein Iterationsverfahren 


sy =flarsy); Ww+ı=glaßy) - » (3.5). 


Dieses Verfahren konvergiert aber . bei hinreiehend 
guten Anfangswerten für eine gesuchte Lösung (&;n) 
sicher, wenn nur die partiellen een vonfundg 
nach z und y im betrachteten Bereich klein genug 
sind. Dafür ist es nach (3.3) unter der Voraus- 
setzung (3.4) am günstigsten, wenn für z=$; 
y=n zugleich 

9 =9d, = Pr. =YP,=0 ..e0.0. (3.6) 


(weil dann auch die partiellen Ableitungen von fund g 
nach x und yan der Stelle (£;n) verschwinden). Dies 
bedeutet aber nach G1.(3.2) bei gegebenen m; n; r; 8 
nichts anderes, als daß (mittels logarithmischer Dif- 
ferentiation wie bei Gl. (1.6)) 


9 [4 Py 1 
etz 917 


Yz e Yy c 
I isn tee 
= - wi 0 


fürz=&; y=n(&*n+ 0 vorausgesetzt). 
Der Sonderfall, daß mit den De Wer- 
ten von u4;v;0;0 die Determinante D=0 wi 
ergibt 


BIP HR) en 


füre=f; y-m 


| - | = | ee 122. @) 
ET, = erh! 5% ; F 
%g 1,621 e 1,621-e ? 0,637 e 5 De ni 
5 ade h F +0,00) er = (z +0,008) a En G 0,065) E 
2 Kg: (£ 40,008) rel Bi 0,074) Ei 


2, im Fall A, z, in den Fällen B und © geben bereits 
a el Een je a Lösung von Gl.(2.1); 
i a wird der teil von 2, positiv, 
in der Nähe von, bleibt. a > 


II 


‚Gegeben seien nun für die (reellen oder auch be- 
liebig komplexen) Unbekannten x und y die beiden 
Gleichungen 


Faey)=0; Cay)=0%)....(. 


?) Die Verallgemeinerung aufn Gleichungen mit kann 
li An der folgenden ohne weiteren Und bedentetim Grund 
ewwas mehr 
Ableitungen 1. Ordnung best = Eat F und @ sollen stetige 


Durch Differentiation 31 
diese Stelle z.B. von (3.1a) folgt aber für 


Ya = mam-1yn _ P,— ERS), 


_:9\_ ,2F%r, 
a 


Da Entsprechendes für F,; G, und @, ei 
G1.(3.8) nichts anderes als £ v gilt, besagt 
RE We.) 
?-Y Fay= 5 GaFy 


9) Wegnr@;m=0, 


oder 3 , 

ne: F2:Fy= 62:4, für =; y=n (3.9), 
d.h. an der betrachteten Stelle eine gemeinsame Tan- 
 gente derKurven F(2;y)=0 und G@(2;y)=0 oder 
_ mindestens für eine der beiden Flächen 2=F(x;y) 
und 2=@(x;y)die xy-Ebene als Tangentialebene (so 

daß (&; n) eine ‚„‚mehrfache Lösung“ wäre). 
Das Ergebnis läßt sich zusammenfassen im folgenden 
Satz 2: Zu jedem Gleichungssystem (3.1) und jedem 


E- Se ‘ Quadrupel von Exponenten m; n; r; 5 gibt 
ec es immer zugehörige Exponenten u; v; 0; o, 


so daß das’ aus (3.3) gewonnene Iterations- 


(x) 


I 3 


Pe 
J 


verfahren (3.5) unter der Voraussetzung 
hinreichend guter Anfangswerte „konvergiert, 
wenn eine im Sinn der Überlegungen zu 
-@1.(3.9) ‚‚einfache‘‘ Lösung mit &n=#0 ge- 
ht ist. BR 5 BT, B 


EV BA 
eit bei der Anwendung des Verfah- 
von chungen wird besonders ein- 
wählen kann, daß. 
en Gleichung 
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Bild 1. Beispiele für Lösungen des Gleichungssystems (4.1). Sowohlin der komplexen @-Ebene wie in der komplexen 
. y-Ebene bedeutet jeweils der Punkt O den Ursprung und der Punkt P den Punkt Y2 auf der reellen Achse 


lichkeit, aus ‚‚schlechten‘“ 1.Näherungen verhältnis- 


‚gelangen, die in der Nähe einer Lösung liegen. 


‚Der E 
Körper in einer Flüssigkeit. 
befaßte sich 8. Neumark!) mit dem Beschleuni- 
gungswiderstand geradlinig bewegter Körper in einer 
id 
ge + 
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Ein Vergleich der einzelnen Glieder in (4.1) läßt 
Beträge von x und y zwischen 1 und 2 vermuten. Wir 
wählen daher als rechnerisch einfache „1.Näherung“ 
die Werte %ı =Yı=)2. Einige mit Hilfe des zu 
(4.2) > gehörigen Iterationsverfahrens ' aus dieser 
‚‚1.Näherung“ gewonnenen Lösungen sind in Bild 1 
eingezeichnet (die Näherungen durch kleine Kreise, 
die richtigen Werte durch kleine Sterne angedeutet). 
Dabei ist natürlich z.B. für C bzw. D die Rechnung 
nicht mehr neu durchzuführen, wenn sie für A bzw. B 
schon vorliegt, da mit (z;y) auch (—2; — y) immer 
eine Lösung von Gl.(4.1) darstellt). 


) 


N 


") 


_ Wird größere Genauigkeit der gesuchten Wertevr- 
langt, so kann es zweckmäßig sein, daß man die Rech- BR: 
nung gegebenenfalls auf das Newtonsche Verfahren 
umstellt, wenn man einmal in der Nähe einer Nul- 
stelle ist. Der entscheidende Vorteil des hier behan- ° 
delten Verfahrens liegt, wie bereits gesagt, inderMög- 


mäßig schnell zu solchen Werten der Unbekannten zu 


München, Et Wenz 1: 


*) Bund F als Beispielefür Lösungen von G1.(4.1),fürdiearez 
und are y stark voneinander abweichen. RR, 


L 3 a > 


_ Beschleunigungswiderstand bewegter 


- 1. In einer interessanten Arbeit aus dem Jahre 1936 


alen Flüssigkeit. Wird der Körper in der x-Rich- 
beschleunigt mit der Geschwindigkeit u=u(l) 
ıvor rühende ideale Flüssigkeit bewegt, 
ekehrt der ruhende Körper in der x-Richtung 
echend beschleunigt angeströmt, so gilt, wenn 
(Druck-)Widerstand im ersten Falle mit W’im 
alle mit-W bezeichnet wird, nach S. Neu- 


a ee 
unigungswiderstand geradlinig beweg, 


on Flüssigkeit; Z.angew. Math. Mech. 6 


a 4 .- F 


88 Kleine Mitteilungen 


wo M die Masse der verdrängten Flüssigkeit ist; von 

äußeren Kräften ist hierbei abgesehen. Wir zeigen im 

folgenden, : 

a) wie man das Neumarksche Ergebnis ohne Rech- 
nung erhalten kann, Zr 

b) daß es auch für eine Klasse nichtgradliniger Bewe- 
gungen gilt, und schließlich 

c) daß dasselbe Ergebnis auch 
dann noch zu Recht besteht, 
wenn man die Zähigkeit be- 
rücksichtigt. Am Schluß folgt 
noch eine Bemerkung zur Föt- 
tingerkuppelung. 

2. x, y, z sei das raumfeste 
System, x’, y’, 2’ das dazu par- 
allele körperfeste System. Der Kör- 
per bewege sich in der x-Rich- 
tung mit der Geschwindigkeit & 
durch die im Unendlichen ru- 
hend gedachte ideale Flüssig- 
keit; sein Widerstand (im &, %, 
2-System!) sei W’. Gehen wir 
jetzt auf das System x’, y’, 2’, so ; 
wirkt dort in der «’-Richtung die scheinbare Kraft 

du 
keit). Der Körper wird jetzt beschleunigt angeströmt, 
und nach wie vor ist der Widerstand W’. ‚‚Löscht“ 
man aber nun das scheinbare Kraftfeld, indem man 
sich die beschleunigte Anströmung allein durch 
Druckkräfte hervorgerufen denkt, so kommt zu dem 


pro Volumeneinheit (o—=Dichte der Flüssig- 


bislang herrschenden Druckgefälle — 2 noch hinzu 


| du ) 2. - 
E35 0° ‚d.h. der Druck erhöht sich von p 


auf p+o er Der neue Widerstand W ist also um 


so viel größer geworden, als der ‚‚Auftrieb‘‘ des Kör- 
pers im vorher bestehenden scheinbaren Kraftfeld 
ausmacht, mit anderen Worten es gilt Gl.(1). 

Diese Betrachtung ergibt von selbst mehr. Wird 
der Körper beliebig bewegt, doch so, daß das körper- 
feste System x’, y’, z’ stets parallel dem System z, y, 2 
bleibt, ist seine Führungsbeschleunigung by, und sein 
Widerstand W’, ist weiter W der Widerstand des 
ruhend gedachten Körpers, wenn man die .betr. An- 
strömung allein durch Druckkräfte hervorruft, so gilt 


EEE PETE N 


wo der Auftrieb des Körpers vom System «’, y’, 2’ 
aus in dem scheinbaren Kraftfeld {— by} (pro Massen- 
einheit) ist. 

Und offenbar gilt die Argumentation auch dann 
noch, wenn die Zähigkeit berücksichtigt wird. 


3. Aus dem Vorstehenden sieht man auch deutlich, 
warum ein ähnlicher Satz nicht mehr gilt, wenn man 
Drehbewegungen des Körpers zuläßt. Denn dann 
treten vom System #’, y’, #’ aus betrachtet Coriolis- 
kräfte auf, die sich nicht aus einem Potential ableiten 
lassen, so daß man die Anströmung jetzt nicht mehr 
durch Druckkräfte allein realisieren kann. In der Tat 
ist ja bei einer Drehung mit dem Drehvektor d die 
Relativströmung von Wirbeln der konstanten Dichte 
— 2b durchsetzt. So betrachtet erscheint es uns auch 
nur bedingt richtig, wenn man, wie das mitunter ge- 
schieht, die Wirkungsweise der Föttinger-Kupplung 
dahingehend erklärt, daß man sagt, die Corioliskräfte 
würden gegen die Wände des oberen Blätterkranzes 
(diesen als ruhend gedacht) drücken; denn eine Druck- 
kraft können diese Kräfte nicht ergeben. Das Bild 
zeigt stark schematisiert eine Kammer des oberen 
Blätterkranzes, welcher an der Arbeitsachse hängt, 
von oben betrachtet. v, ist die (relative) Geschwin- 
digkeit, & deutet die Corioliskraft und 8 die Zen- 
trifugalkraft an; letztere wird hier im Gegensatz zu 
oben nicht ausgelöscht. Wendet man nun den Impuls- 
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satz in der y-Richtung auf den in dem Bild gestrichelt 
abgegrenzten Bereich an, so liefert hierfür weder der 
Impulstransport noch die Zentrifugalkraft (in erster 
Näherung) einen Beitrag, jedoch die Corioliskraft 
einen Beitrag —0. Dieser Beitrag muß also allein 
durch die Druckkräfte auf den gestrichelten Flüssig- 
keitskörper kompensiert werden, d.h. also umgekehrt 


Bildi 


die Flüssigkeit drückt alles in allem auf die Wände 
nach oben in der y’-Richtung. 


Berlin. Ernst Mohr. 


Ermittlung der Wirk- und Blindleistung bei 
mittelbarem Antrieb des einfachen Schwingers. 


Die Kenntnis der Wirk- und Blindleistung des An- 
triebes ist von Interesse z.B. bei der Bestimmung des 
Wirkungsgrades des Antriebes von schwingungstech- 
nischen Arbeitsmaschinen und der Bemessung der 
Größe des Schwungrades zur Begrenzung des Un- 
gleichförmigkeitsgrades auf ein bestimmtes Maß. Un- 
ter Wirkleistung l„(l) verstehen wir das Produkt der 
mit der Geschwindigkeit v(tl)=V cos (Ri+e) des 
Kraftangriffspunktes in Phase liegenden Komponente 
p’(t) der Antriebskraft p(t) = P cos (Ri), unter Blind- 
leistung l;ı(t) das Produkt der mit v(t) in Gegenphase 
befindlichen Komponente 9’’(t) der Antricheirei 


= (+ 008290; In=— 7 sin2Q. 


Bei mittllbarem Antrieb ist die Auslenkung «(t) 
bzw. die Auslenkungsgeschwindigkeit des Kraftan- 
griffspunktes gegeben, die Größe und Phase der An- 
triebskraft selbst dagegen nicht. Zunächst behandeln. 
wir das einfache Beispiel des durch eine zweite Feder 
mittelbar angetriebenen Schwingers (vgl. Bild 1). Der 
Schwingungsvorgang findet seine mathematische For- 


‚mulierung in der D-Gleichung 


eg rbgtag=aWw—g)..... (l), 
die auch wie folgt geschrieben werden kann 
oa +bga + to)g=uu... . (le), 


Für die Bestimmung der Auslenkung, Auslenkungs- 
geschwindigkeit und -Beschleunigung der Masse a ist 
also ein ra entsprechend Bild 2 gleichwertig 
mit dem in Bild 1. Auch für die Berechnung der Wirk- 
leistung sind in diesem Falle Bild 1 und 2 einander 
Beh weei da die Wirkleistung hier die im Dämp- 

er bin Wärme umgesetzte Antriebsleistung bedeutet, _ 
die lediglich von der Auslenkungsgeschwindigkeit der 
Masse a abhängt. Bei der Ermittlung der Blindleistung 
führt jedoch im vorliegenden Falle die Benutzung des 
Ersatzsystemes gemäß Bild 2 zum falschen Ergebnist). 
Man sollte daher stets die Bestimmung der Leistungs- 
größen auf Grund des „Original“-Ersatzsystems 
(Bild 1) und nicht unter Bezugnahme auf ein „abge- 
leitetes“ Ersatzsystem (Bild 2) durchführen. 


2. Aufl.. Berlin, Göttingen, Heidelberg 1951, 8. 213, 


_— 


!) Vgl. K. Klotter: Technische Schwingungslehre, 1. Bd., 
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In komplexer Schreibweise läßt sich die komplexe 
Amplitude ® der Antriebskraft und die komplexe Am- 
plitude ® der Geschwindigkeit ihres Angriffspunktes 
wie folgt ausdrücken?) 


ee) 


mit 
Ba ze 
TI H+4DN 70 
Vin. 


Darin bedeuten 


c=c,+ c, Federkonstanten 

y, die hier gültige reduzierte kinetische Einflußzabl 
2 Kreisfrequenz der Erregung 

& Kreisfrequenz der Eigenschwingung des Systems 


D= 5 2 dimensionsloses Dämpfungsmaß, 
ac 
Die mit ® in Phase liegende Komponente der An- 
triebskraft beträgt also 
A ; 
en! (2a), 


e dm HDi" 


Ber; 


Cult) 


eg, 
pl): 
6,.u(t) 
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Zur Ermittlung der den Gl. (3a) und (3b) entspre- 
chenden Ausdrücke bei -mittelbarem Antrieb des 
Schwingers über einen zusätzlichen Dämpfer (vgl. 
„Original“-Ersatzsystem Bild 4) darf wiederum nicht 
das ‚‚abgeleitete‘‘ Ersatzsystem nach Bild 5 der Rech- 
nung zugrunde gelegt werden. 

Als komplexe Amplitude der Antriebskraft erhalten 
wir 

b 
P=ih,2(U-QA)= u _ vs) 
mit 
2Dnf2Dn+i(l—n?)] 
y,= (LP F4Dim -’ 
für die Geschwindigkeit ihres Angriffspunktes 
HAIR 


Darin bedeuten 
b=b, + b, Dämpfungskonstanten 
4, die hier gültige reduzierte kinetische Finflußzahl. 


Entsprechend Gl. (2a) und (2b) ergeben sich die in 
Phase bzw. in Gegenphase mit ® befindlichen Kom- 


Bild 3 Bild 4 


Bild 5 Bild 6 , B Bild ? 


die in Gegenphase befindliche 


wre u(ı aa, mn) (2b). 
er Led 
Daraus folgt nach einfachen Umformungen als Am- 
plitude der Wirkleistung?) 
03 .c3 2Dn? 
Lv 


Say Amer 


und als Amplitude der Blindleistung?) 


Be, nd 3b). 
ae e (G—-m?+4D?r? (3b) 


Für den Fall, daß der mittelbare Antrieb nicht über 


eine Zusatzfeder erfolgt (vgl. Bild 3) ist in Gl. (3a) und 
(3b) c, durch c zu ersetzen. 


2). Vgl, ebenda, 8.183 ff. 


ponenten der komplexen Kraftamplitude 


ge b 4.D2 m 
Wil (1- 7 - rad) (4a) 
WON ram (4b) 


Wir erhalten Somit für die Amplitude der Wirk- 
bzw. Blindleistung®) nach einfacher Zwischenrechnung _ 


v2 1 /c® (b,\ SD 2 2 (53) 
W»=7Yy,\) ap +raDm b, 


u21/&(b, 4De®(l m) 
Ze a\b (1—2)?+ 4D?n? 


an 8 (Dh), 


3) Vgl. G1.(63.8) auf 8.213 des angezogenen Buches von 
K. Klotter. ; 
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Erfolgt der mittelbare Antrieb ohne Zusatzdämpfer 
(vgl. Bild 6), so ist in Gl. (5a) und (5b) b, durch b zu 
ersetzen. 

Als letzten Fall untersuchen wir den des Antriebes 
durch Massenkräfte (vgl. Bild 7). Die Resultierende 
der komplexen Amplituden der Beschleunigung der 
rotierenden Massenpunkte a, beträgt 


—- (OU :...:... (6). 
Die komplexe Amplitude der an a, in Schwingungs- 
richtung angreifenden Antriebskraft zur Überwindung 


des durch die Beschleunigung (6) hervorgerufenen 
Trägheitswiderstandes ergibt sich also zu 


Ba, HN) =a,NRU ( 2) . 9) 


mit 
a=4+0, 
und 
__TA—-d)+2Drfi 
Ya (1 — ?)2 + 4 .D2n? 

Die mit der komplexen Amplitude iQ U der Relativ- 
geschwindigkeit der umlaufenden Massen in Phase 
bzw. in Gegenphase befindliche Komponente von 
(7) ist damit 

Ve N he 2Dn?i 
B = 2? u 1-7 —P)? +4Dn% a: 4. Dir? u (8a) 


bzw. 


siert 
= 1,21 f nz r (1 Se Eu u „| (8b) 


Nach einigen einfachen Umformungen erhalten wir 
dann für die Amplitude der Wirkleistung®) 


508 & (a,\” 2Dn® 
=: %) am rad 08) 
und für die Amplitude der Blindleistung®) 


_Ul,)ı/e|a , Pl—?) 
iu-3:(@) VE a" Famradm| PD 
Hannover. E. Pestel.. 


*) Vgl. G1.(63.9) auf 8.213 des angezogenen Buches von 
K. Klotter, 


DerKrümmungseffekt beidickwandigen Hoch- 
druckrohren beliebiger asymmetrischer 
Querschnittsformen!). 


Zur Messung hoher Flüssigkeitsdrücke werden in 
der Technik neben elektrischen Verfahren vor allem 
rohrförmige elastische Meßorgane benutzt. Bei einem 
in neuerer Zeit bekannt gewordenen elastischen Hoch- 
druckmeßorgan nutzt man die Durchkrümmung eines 
exzentrisch gebohrten Rohres unter der Wirkung eines 
Innendruckes aus. Dieser Effekt kommt dadurch zu- 
stande, daß der Schwerpunkt O (Bild 1) der Ringfläche 
und der herr M der auf die beiden Enden 
wirkenden Axialkraft nicht zusammenfallen, so daß 
ein Biegemoment entsteht, das eine Krümmung des 
ursprünglich geraden Rohres verursacht. Diese Wir- 
kung wird jedoch durch den Einfluß der Querkontrak- 
tion geschwächt, da auf der schwächeren Wandseite, 
wie in Bild 1 durch die Länge der Pfeile angedeutet 
ist, die Seit, größer sind als auf der 
stärkeren Wandseite. Damit werden aber wegen der 
Querkontraktion die Axialspannungen auf der schwä- 
cheren Wandseite verstärkt, auf der anderen Seite ge- 
schwächt. Die genaue Rechnung?) hat gezeigt, daß 
der Einfluß der Querkontraktion sich bei einem ex- 
zentrisch gebohrten Rohr in einem Faktor (1—2/m) 


*) Vorgetragen auf dem Mechanik-Kolloquium der Nieder- 
sächsischen Hochschulen in Göttingen 7. Juli 1951. 

») W. Wuest; Theorie des Hochdruckmeßrohres mit ausmit- 
tiger Bohrung, Ing.-Arch. 19 (1951), 8. 12—21. 
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auswirkt. Versuche haben im Einklang damit gezeigt, 
daß der ursprüngliche Krümmungseffekt, der ja ele- 
mentar zu berechnen ist, durch die Querkontraktion 
auf die Hälfte herabgesetzt wird. Die Vermutung liegt 


2 


Bild 1. Krümmungseffekt beim exzentrisch gebohrten Rohr 


nahe, daß der Einfluß der Querkontraktion sich ganz 
allgemein bei beliebigen asymmetrischen Querschnitts- 
formen des Hochdruckrohres in einem Faktor (1—2/m) 
äußert,”mit dem der primäre Krümmungseffekt zu 
multiplizieren ist. Diese Vermutung wird durch die 
folgende Untersuchung bestätigt. 

Wir nehmen also einen beliebigen Rohrquerschnitt 
an, von dem wir zunächst nur voraussetzen, daß er 
eine einzige Bohrung enthält, die von einem Flüssig- 
keitsdruck p beaufschlagt wird. Die Querschnitts- 
fläche ist also zweifach zusammenhängend. Die Länge 
des Rohres sei groß gegenüber den Querschnitts- 
abmessungen, so daß man mit genügender Näherung 
einen ebenen Formänderungszustand voraussetzen 
kann. Die beiden Rohrenden üben offenbar eine Axial- 
kraft pQ; aus, wenn man mit Q; denlichten Innenquer- 
schnitt bezeichnet. Weitere äußere Kräfte oder Mo- 
mente seien dagegen nicht vorhanden. Das Koordi- 
natensystem kann in der Querschnittsebene natürlich 
beliebig gelegt werden. Um die Rechnung zu verein- 
fachen, legen wir es so, daß der Koordinatenursprung 
mit dem Schwerpunkt O der Ringfläche Q zusammen- 
fällt und die y-Achse durch den Schwerpunkt M der 
Innenfläche Q; geht. Die Axialspannung o, ist dann 
gegeben durch: 


%=(k+iy+ ja) B+l(@y+o,). .Ä{). 


Hierbei ist % die noch unbekannte axiale Dehnung 
und i,j sind die Krümmungen in y und z2-Richtung 
(d.h. die reziproken Krümmungsradien in der zy- und 
xz-Ebene). 

Bei Einführung der Airyschen Spannungsfunktion F 
ist die Spannungssumme in bekannter Weise durch 
0y+0z=4AF gegeben. Die drei Konstanten sind aus 
folgenden Gleichgewichtsbedingungen zu bestimmen: 


S 0udQ = (ip 
Say dr (2) 
JS 022:dQ =0 f 


oder nach Einsetzen von Gl. (1): 
QEE+SAFQ=-Gp — 
2 1 
WiE+WE+ SAFyR=Qap\ (9 
| Int E+jI,E+ 2 SAR2dQ =0 
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Hierbei bezeichnet Q die Ringfläche und dQ ein 
Flächenelement der Ringfläche, J, und J, sind die 
Trägheitmomente, Jy, das Zentrifugalmoment der 
Ringfläche. Die Integrale sind über die Ringfläche zu 
erstrecken. Die Flächenintegrale können aber sofort 
in Integrale über die beiden geschlossenen Randkurven 
umgeformt werden: 


Fr or 
[ara = (+ za) ara 
Q Q 


27 2 
[arsao = (+ 35) vava 
Q 


er oF or 
== rede a 
8 


8 


or oF 
[Ara (+ 32) 7% dz 
Q 


Q 
or oF 
plane 


Die Vorzeichen sind dabei so gewählt, daß die beiden 
Randkurven im Zeigersinn von Bild 2 zu umlaufen 
sind. Die Spannungsfunktion F und ihre ersten par- 


Bild 2. Rohrquerschnittbeliebiger Form. Der Koordinatenursprung 
ist in den Schwerpunkt O der Ringfläche @ gelegt 


tiellen Ableitungen OF/dy und OF/öz sind aber auf den 
Randkurven durch die Randbedingungen gegeben. Es 
gelten nämlich bei ebenem Spannungs- oder Form- 
änderungszustand folgende Oberflächenbedingungen: 


oy cos (n,Y) + Tzy cos (n,2) = Py 

Pr LT) 
TyzC08 (N, Y) + 02 COS (N,2) = Pz 

"Setzt man (vgl. Bild 3): 


02 
_ c08 (n,y) = 35; 
cos (n, 2) = — 23 5 
so gehen die Bedingungen GI. (7) mit 
0°F o2F 
ua «za 
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und 
BES RL EN 
1 ze 0902 
über in: 
er Dry ar 
02° 08  Oyozds 0%: 7 
SFr 0roy Dr 


020y0s ' Oy? ds say. ei 


Bei reiner Normalbelastung der Oberfläche durch 
einen gleichmäßigen Flüssigkeitsdruck ist aber: 


02 
I I Fe FR 
d 
m=—pcos(n.= pP 


so daß also auf der inneren Randkurve gilt: 


or 
| ars =—py+Cı.. (8) 
8 


or 
a | re t+a= net 8) 
8 


or or 
(tz) 
8 


= (y+P)+09y+ 03240, 


(10). 


An der äußeren Randkurve ist p=0. Jedoch tritt 
dort ein weiteres Konstantentripel C7, C%, 05 auf. 
Ebenso treten bei mehr als einer Bohrung an jeder 
neuen Randkurve weitere Konstantentripel auf, von 
denen aber eines willkürlich gewählt werden kann. 
Wir setzen also im vorliegenden Fall nur einer Boh- 


Z n 


y 


Bild 3. Normalenrichtung an der Berandung des Ringquerschnittes 


rung (4 =C}=(C7=0 und brauchen die Integration 
also nur noch über die innere Randkurve zu erstrecken. 


Das Integral Gl. (4) nimmt damit folgende Forman: 
[Araydz = $(—py+ CD) de— f(—pz+ 0,)dy. 
Q 8 8 


Die Glieder mit C, und C, verschwinden bei der In- 
A : dz 3 
tegration und mit da = (2) dy erhält man: 


ke —2p fzdy=2pQ:. 
U) 


Damit liefert die erste Gl. (3) für die axiale Deh- 


‚nung k die einfache Beziehung: 


N 
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Ganz entsprechend können auch die Integrale 
Gl. (5) und (6) ausgewertet werden. 


[ Ar ydydz = fy(—py+ Cdde 
() 8 
+f2W+9)-0y-02—-0 dz 
8 
+ $u(pz— 0,)dy 
8 
= f2pyzdy=2padi. 
8 


Hierbei ist a wiederum der Abstand des Schwer- 
punktes M der Innenfläche Q; vom Schwerpunkt der 
Ringfläche. 


In entsprechender Weise findet man 


[ Arzdya = 2pfyrda=0. 
Q 8 


da die y-Achse eine Schwerpunktsachse des Innen- 
querschnittsist. Man erhält damit für dieKrümmungs- 
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komponenten: 
Ne 2\ aQids pP 
lm) —I, EB 


a 1 aQiJyz pP 
Ir, m)JyJz — JS? E° 


Bei einem Ringquerschnitt mit einfacher u 
trie (z.B. exzentrisch gebohrtes Rohr) ist das Devia- 
tionsmoment Jyz=0 und die Krümmungsebene fällt 
mit der Symmetrieebene zusammen. 2 

Die vorstehenden rlegungen können ohne wei- 
teres auch auf mehrfach durchlochte Querschnitte an- 
gewandt werden, da die Konstanten C,, C,, C, im Er- 
gebnis nicht vorkommen. Auch wenn die verschie- 
denen Bohrungen mit verschiedenen Flüssigkeits- 
drücken beaufschlagt werden, äußert sich die Quer- 
kontraktion stets in einem Faktor (1—2/m) für die 
Krümmung. Diese Gesetzmäßigkeit gilt natürlich 
nicht, wenn das Rohr nicht mehr durch reinen Flüssig- 
keitsdruck beaufschlagt wird. Doch können auch in 
diesem allgemeineren Fall die Längsdehnung und die 
Krümmungen durch einfache Quadraturen berechnet 
werden, ohne daß die Spannungsverteilung im ein- 
zelnen bekannt ist. 


Göttingen. Walter Wuest. 
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Arthur Linder (Prof. a. d. Univ. Genf und der ETH. 
Zürich), Statistische Methoden für Natur- 


: wissenschaftler, Mediziner und Ingenieure 


(Lehrbücher und Monographien aus dem Gebiete der 
exakten Wissenschaften. Mathematische Reihe Bd. 3), 
2. erweiterte Auflage. 238 S. mit 42 Abb.: Basel 1951. 
Verlag Birkhäuser. Preis brosch. 26 Fr., geb. 30 Fr. 


Zweck dieses Buches ist, eine Einführung in die 
neueren statistischen Prüfverfahren zu geben, die in 
den letzten dreißig Jahren in den angelsächsischen 
Ländern entwickelt sind. Da das Buch sowohl für 
den Praktiker wie für den theoretisch Interessierten 
bestimmt ist, wird zunächst in den ersten drei Ab- 
schnitten (Statistische Maßzahlen, statistische Prüf- 
verfahren, Streuungszerlegung) die Art der Anwen- 
dung der Verfahren in der Hauptsache an Beispielen 
gezeigt, während erst das letzte (Theorie der Stich- 
proben) die mathematischen Ableitungen durch- 
gehend unter Benutzung der n-dimensionalen Geo- 
metrie gibt. Obwohl sich Verfasser auf die grund- 
legenden Verfahren beschränkt, scheint mir die Dar- 
stellung in diesem letzten Abschnitt etwas zu kon- 
zentriert; man möchte wünschen, daß etwas mehr 
Zwischenrechnungen eingefügt würden. Am Schluß 
hat Verfasser eine Reihe von Tafeln insbesondere 
solche der verschiedenen Standardverteilungen an- 
gefügt, die für die praktische Anwendung erforderlich 
sind. 

Die vorliegende zweite Auflage ist gegenüber der 
ersten etwas erweitert. Neuaufgenommen sind Ab- 
schnitte über die Regressionsmethoden, die Trenn- 
verfahren, den verallgemeinerten Abstand und ins- 
besondere über die Streuungszerlegung. — Das Buch 
kann sowohl dem aktibchen Statistiker wie dem 
Mathematiker empfohlen werden, 


Dresden. Willersz 


R. Sauer (o. Prof. a. d. TH. München), Ecoule- 
mentsdesfluidescompressibles. XVI 


+ 307 8. mit 150 Abb. Paris et Liöge 1951. Librairie 


polytechnique Ch. B6ranger. Preis geb. 3,900 fr. 
Der Autor hatte 1943 im Verlag Springer, Berlin, 


‘ die „Theoretische Einführung in die Gasdynamik“ 


herausgegeben (1951 in zweiter Auflage erschienen). 
Diese beschränkte sich jedoch nur auf die stationären 


Probleme. Als Ergänzung hierzu sollte gegen Ende 
des Krieges ein Buch erscheinen, das die nicht- 
stationären Probleme zusammenfaßte. Es war bereits 
gesetzt, doch ist der Ausdruck nicht mehr erfolgt. 

Die Mußestunden, die die erste Nachkriegszeit dem 
Autor zwangsweise brachte, hat er dazu benutzt, die 
beiden Abhandlungen zu einer einzigen zusammen- 
zufassen. Diese hat daher den großen Vorteil, daß die 
nichtstationären Probleme von vornherein an ge- 
eigneten Stellen mit eingearbeitet sind. Dieses 
Manuskript ist dann ins Französische übersetzt 
worden. 


Gegenüber der „Einführung in die theoretische. 


Gasdynamik“ (1951) enthält das Buch folgende neuen 
oder wesentlich erweiterten Kapitel: 

Grundlagen der instationären Rohrströmung;? Po- 
tential- und Stromfunktion der isentropischen Strö- 
mung; linearisierte Druckwellen; stationäre} Über- 
schallströmung mit Rotationssymmetrie oder? nicht 
isentropisch;; isentropische Druckwellen großer Inten- 
sität; starke Rohrwellen bei variablem Querschnitt 
oder variabler Entropie; Strömungsdiagramm für 
K=+1lund K=3; Gebrauchsanweisung für die 
Charakteristiken-Theorie; Grundgesetze derf Stoß- 
welle; ebene Strömung mit homogenen Druckwellen; 
ebene und sphärische en 

Das Buch bietet in einer klaren Form eine aus- 
gezeichnete Einführung in ein modernes Gebiet, das 


sonst nur in einer Reihe oft unzugänglicher Artikel in. 


wissenschaftlichen Zeitschriften oder — noch weni 
erreichbar — in Berichten wissenschaftlicher Institute 
seinen Niederschlag gefunden hat. 


Weil a. Rh. Schardin. 


R. Oleg (Ph.D., MathematicalEn- 
gineering Analysis. XIV-+426 S. mit 
220 Abb. New York 1950. The Macmillan Company. 
Preis 6,— $. - 

In der Praxis liegen die Schwierigkeiten für die 
Anwendung der Mathematik auf physikalische und 
technische Probleme vielfach im Herausarbeiten des 
mathematischen Kernes und der Aufstellung von 
Gleichungen, die der mathematischen Behandlung 
zugänglich sind.. Dazu sind nicht nur gründliche 
Kenntnisse der theoretischen Physik und.der in Frage 
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kommenden technischen Wissenschaften, sondern 
auch eine“ langjährige Erfahrung nötig, wie sie der 
Verf. hat. Der Zweck dieses Buches ist, diese Erfah- 
rung den Studierenden zu übermitteln. In fünf Ka- 
piteln ‚Mechanik starrer Körper, Elektrizität und 
Magnetismus, Wärmelehre, Elastizität und Strömungs- 
lehre“ stellt Verf. zunächst die grundlegenden Gesetze 
der einzelnen Gebiete dar und zeigt anschließend an 
einer großen Zahl von Problemen, wie die mathe- 
matischen Ansätze zu machen sind, wie man dabei zu 
Gleichungen zwischen dimensionslosen Größen kommt, 
welche Vernachlässigungen in den einzelnen Fällen 
erlaubt sind, damit die mathematische Behandlung 
des Ansatzes mit den üblichen Methoden möglich 
wird usw. Auf die mathematische Behandlung selbst 
geht er nicht ein. Am Ende eines jeden Abschnittes 
sind eine große Anzahl von Problemen zusammen- 
gestellt für die am Schluß des Bandes die mathe- 
matischen Formulierungen bzw. Hinweise auf ihre 
Behandlung gegeben sind. So hat Verf., fußend auf 
seine längjährigen Erfahrungen, ein für den jungen 
Techniker sehr nützliches Buch geschaffen. 


Dresden. ‘:Willers. 


M. V. Wilkes, D. J. Wheeler and $. Gill, The 
preparation ofprogramsforanelec- 
tronie digital computer, with spezial re- 
ference to the EDSAC and the use of a library of sub- 
routines. X + 1708. m. 2 Abb. Cambridge (Mass) 
1951. Addison-Wesley Press. Preis geb. 5,— $. 

Bei Ausführung längerer Rechnungen treten immer 
wieder die gleichen Operationsfolgen auf. Diese kann 
man für die Ausführung mit einem Rechenautomaten 
ein- für allemal programmieren, um sie als Teilpro- 
gramme dann den einzelnen Gesamtprogrammen ein- 
zufügen. Das bietet den Vorteil, daß dadurch die Ar- 
beit des Programmierens und die Fehlermöglichkeit für 
diesen Teil der Rechnung beträchtlich vermindert wer- 
den. Das vorliegende Buch enthält nun ins einzelne 
gehende Beschreibungen der Sammlung von Unterpro- 
grammen, wie sie in Cambridge für den EDSAC benutzt 
werden, und schildert die Art, wie die Programmierung 
mit ihrer Hilfe vorgenommen wird, Gerade in Cam- 
bridge hat man auf diesem Gebiet wohl die reichsten 
Erfahrungen und verfügt über die umfangreichste 
Sammlung solcher Unterprogramme. Diskutiert wird 
hier der beste Weg, solche Unterprogramme für nume- 
rische Quadratur, genäherte Berechnung von Integral- 
kurven und andere Prozesse herzustellen. Mehr theo- 
retische Fragen, wie die Konvergenz iterativer Pro- 
zesse, die Anhäufung der Abrundungsfehler usw. wer- 
den nicht behandelt. Die Darstellung bezieht sich zu- 
nächst auf den EDSAC, doch sind Ideen und Technik 
auch für andere Maschinen, die nach den gleichen 
Grundsätzen gebaut sind, anwendbar, so daß jeder, der 
mit großen Rechenautomaten zu arbeiten hat, das Buch 
mit Vorteil benutzen wird. 

Dresden. Willers. 

Dr. Ing. habil. Kurt Nesselmann, „Die Grund. 
lagenderangewandtenThermodyna- 
mik“, XI + 320 S. mit 311 Abb. und 5 Diagrammen, 
Berlin-Göttingen-Heidelberg 1950. Springer-Verlag. 
Preis geb. 18,— DM. 

Der Verfasser unternimmt wohl erstmalig den Ver- 
such, die gesamte technische Thermodynamik ge- 
schlossen in einem einzigen Buche darzustellen. 
Während in der Regel die Autoren sich auf die all- 


a Thermodynamik, d.h. die Lehre von den, 


asen und Dämpfen, beschränken, allenfalls noch 
die Wärmeübertragung streifen, bezieht der Verfasser 
in der richtigen Erkenntnis der immer mehr in den 
Vorder dtretenden chemischen Erfahrungstechnik 
auch die sogenannte chemische Thermodynamik ein. 
So behandelt das Buch außer der allgemeinen Ther- 
modynamik und der Wärmeübertragung u.a. die 


Gleichgewichte der Mehrstoffgemische der festen, 
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flüssigen und gasförmigen Phase, die destillativen 
Trennverfahren, die Wärmetönung, das chemische 
Gleichgewicht einschl. dritten Hauptsatz und den 
Stoffaustausch durch Diffusion. 

Beirder Fülle des Stoffes konnte es nicht ausbleiben, 
daß das Werk in manchen Teilen infolge der Gedrängt- 
heit der Darstellung skizzenhaft erscheint; verschie- 
dentlich wäre ein tieferes Eindringen in die Materie 
wünschenswert gewesen. Manche wichtigen Gebiete 
wurden allzu knapp behandelt (z. B. die Reaktions- 
kinetik auf einer halben Seite, die Gasstrahlung auf 
einer Seite), andere fehlen ganz (z. B. die destillative 
Trennung durch Teilniederschlag, die Solventtren- 
nung, die unvollständige Verbrennung). Das Buch 
enthält bedauerlicher Weise eine große Zahl von 
Druckfehlern, die das an sich nicht einfache Studium 
empfindlich stören, und ist auch nicht ganz frei von 
Irrtümern. Trotz dieser bei einer Erstbearbeitung 
verständlichen Schwächen muß anerkannt werden, 
daß es dem Verfasser durch straffe, übersichtliche 
Gliederung, weise Beschränkung auf die rein thermo- 
dynamischen Entwicklungen, verbunden mit einer 
klaren Darstellung des Stoffes, und die Ergänzung 
desselben durch eingestreute geschickt gewählte 
Rechenbeispiele im ganzen gelungen ist, das gesteckte 
Ziel einer einheitlichen Behandlung der gesamten 
technischen Thermodynamik zu erreichen und der 
Fachwelt .damit ein, dankenswertes Hilfsmittel zu 
schenken. 


Dresden. Faltin. 


Th. L. Wade, The Algebra of Vectors 
and Matrices. XII-+ 189 S. m. 10 Abb. Cam- 
bridge (Mass.) 1951. Addison Wesley Press. Preis geb. 

Das Buch bietet eine elementare, relativ einfache, 
gut lesbare Einführung in die für alle Anwendungs- 
gebiete wichtige Algebra der Vektoren und Matrizen, 
wobei die Grundbegriffe der modernen Algebra, wie 
Gruppe, Feld, Ring, Basis, Dimension, Isomorphis- 
mus usw. eingeführt werden, so daß es sowohl als Vor- 
bereitung für die höhere Algebra wie die Vektor- und 
Matrizenrechnung dienen kann. Vorausgesetzt wird 
eine gewisse Vertrautheit mit der analytischen Geo- 
metrie und den einfachsten Sätzen über Determi- 
nanten. Eine Übersicht über den behandelten Stoff 
mögen die Kapitelüberschriften geben: Grundbe- 
griffe; Vektoren im zwei- und im dreidimensionalen 
Raum; Vektormethoden in der analytischen Geo- 
metrie; lineare Abhängigkeit von Vektoren; Vektoren 
in n-Dimensionen; elementare’ Eigenschaften von 
Matrizen; spezielle Matrizen; Gruppen, Matrizen und 
Transformationen; die charakteristische Gleichung 
einer Matrix; Matrizen. und algebraische Formen; 
einige Anwendungen der Matrizenalgebra. Überall 
eingefügte Übungsaufgaben dienen der Einübung und 
der Vertiefung des dargebotenen Stoffes. 

Das ausgezeichnete Büchlein gibt so den Stoff einer 
Einführungsvorlesung, wie sie für die ersten Semester 
gehalten werden sollte. Es sei vor allem den Anfangs- 
semestern empfohlen. 


Dresden. Willers. 

Dr. F. Neiss, (Prof. a. d. Humboldt-Univ. Berlin), 
Analytische Geometrie. VIII + 167 8. mit 
64 Abb. Berlin-Göttingen-Heidelberg 1950. Springer- 
Verlag. Preis brosch. 9,60 DM. 

Das vorliegende Buch arbeitet weitgehend mit den 
Hilfsmitteln der linearen Algebra. Da es jedoch als 
Lehrbuch für Studierende gedacht ist, ist der eigent- 
lichen Darstellung eine elementare Einführung voran- 
gestellt, in der noch keine höheren Hilfsmittel ver- 
wendet werden, die aber die Tragweite der elemen- 
taren Methoden aufzeigt und dem künftigen Lehrer 
wertvolle Hinweise für die Gestaltung seines Unter- 
richts.geben kann. Für den wesentlichen Teil seines 
Buches befindet sich Verf. in der günstigen Lage, 
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ständig auf sein bewährtes Buch über ‚‚Determinanten 
und Matrizen‘ hinweisen zu können. Es gelingt so, 
den gesamten Stoff auf reichlich 100 Seiten zu bringen 
und den Leser davor zu bewahren, daß die Aufmerk- 
samkeit durch Erörterungen über den Kalkül ab- 
gelenkt wird. Die nötigen Rechnungen können mit 
großer Eleganz durchgeführt werden, und der geome- 
trische Inhalt steht stets im Vordergrund. 

Nach einer Überleitung (Geometrie der Geraden 
und Ebene) folgen zwei Kapitel über kongruente und 
ähnliche Abbildungen sowie über die projektive Geo- 
metrie der linearen Gebilde, das Zentralstück des 
Buches hinsichtlich Umfang und Bedeutung. An- 
schließend werden die Kurven und Flächen zweiter 
Ordnung behandelt. Jedem Kapitel sind eine Anzahl 
Übungs- und Ergänzungsaufgaben beigegeben, z. B. 
Aufgaben zu den kongruenten Abbildungen, die auf 
gruppentheoretische Gesichtspunkte hinweisen. 

So ist ein Buch entstanden, das der Leser nicht nur 
mit großem Gewinn, sondern auch mit großem Genuß 
durcharbeiten wird, Die Ausstattung des Buches 
hinsichtlich Druck, Papier und Deutlichkeit der 
Abbildungen ist sehr gut. 


Dresden. G. Opitz. 


Heinz Rutishauer, Ambros Speiser, Eduard Stiefel, 
Programmgesteuerte digitale Rechengeräte 
(elektronische Rechenmaschinen). (Mitteilungen 
aus dem Institut für angewandte Mathematik an der 
eidgenössischen Technischen Hochschule in Zürich, 
Nr.2.) 102 S. mit 28 Abb. Basel 1951. Verlag Birk- 
häuser. Preis 8,50 Schweizer Franken. 


Es handelt sich bei dieser Veröffentlichung um 
einen Separatdruck eines Berichtes der Verfasser in 
der Z. angew. Math. Phys., Bd.I, S. 277—297, 340 
bis 362; Bd.II, S.1—25, 63—92. 

Erstmalig wird hier von berufener Seite in deut- 
scher Sprache ein umfassender Überblick über die 
Entwicklung und den Aufbau programmgesteuerter 
(meist elektronischer) Ziffernrechenautomaten ge- 
geben. 

Nach einer Einführung in die allgemeinen Grund- 
prinzipien werden Organisation und Arbeitsweise, die 
arithmetischen Prinzipien, Fertigung von Rechen- 
plänen und schließlich auch die physikalischen Grund- 
lagen bestebender Maschinen behandelt. Überall ist 
dabei die reiche persönliche Erfahrung der Verfasser 
über den Gegenstand spürbar. In übersichtlicher und 
straffer Darstellung wird der Leser bis unmittelbar an 
die zugehörige Spezialliteratur herangeführt, deren 
Aufsuchen ein umfangreiches Literaturverzeichnis 
erleichtert. 

Eine Tafel mit 18 bekannt gewordenen Rechen- 
'automaten und Mitteilungen über behandelte Auf- 
gaben illustrieren Bedeutung und den rasch zuneh- 
menden Umfang dieser automatisierten Rechen- 
technik. Esist hervorzuheben, daß in diesem Bericht 
neben den amerikanischen und englischen Pionier- 
arbeiten auch die davon unabhängigen und erfolg- 
reichen Konstruktionen des deutschen Ingenieurs 
Konrad Zuse gewürdigt werden. 


Dresden. N. Joachim Lehmann. 


Lynde Phelps Wheeler. Josiah Willard 

ibbs. The history ofa great mind. 
XI + 2648. mit 9 Taf. New Haven 1951. Yale 
University Press. Preis 4,— $. 

Nachdem Gibbs’ Schüler Lynde Wheeler 
(Assoc. Prof. of Physics an der Yale University) vor 
vier Jahren an der Schrift über ‚‚The early work of 
Willard Gibbs in applied mechanics‘“ mitgewirkt hat, 
legt er jetzt zur 250-Jahrfeier von Yale eine gerundete 
ns von J.W. Gibbs (1839—1903) vor. 
Mit Fleiß und Umsicht hat er Gib bs’ wissenschaft- 
lichen Nachlaß und Familienpapiere benutzt und ein 
auch das Menschliche treffendes Bild von Gibbe’ 
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Leben und Werk gezeichnet. Wir sehen nunmehr, daß 
Gibbs nicht der ‚‚Mathematiker in einem Turm“ 
gewesen, sondern in regem geistigen Austausch mit 
den Fachgelehrten der Neuen und der Alten Welt ge- 
standen hat, und daß Gi b b s’ Bedeutung als Thermo- 
dynamiker und als Vater der Phasenlehre schon zu 
Lebzeiten von Männern wie Lord Kelvin, 
Wilhelm Ostwaldund JamesClerkMax- 
wellerkannt worden ist. 

Hingewiesen sei noch‘ auf die z.T. methodisch 
interessanten Anhänge des Buches: 1. Genealogische 
Tafel; 2. Gibbs’ erste bislang unveröffentlicht geblie- 
bene Arbeit vom Jahre 1866 über ‚‚The proper magni- 
tude of the units of length, and of other quantities 
used in mechanics‘; 3. Katalog der Wissenschaftlichen 
Korrespondenz von Gibbs; 4. Gibbs’ Postlisten 
seiner von ihm an Personen und Zeitschriften- 
Redaktionen verschickten Bücher und Abhandlungen; 
5. Referat für Koenigsbergers ‚‚Repertorium‘‘ über 
Gibbs’ Abhandlung ‚‚On the fundamental formulae 
of dynamics‘ (1879); 6. Bibliographie der Arbeiten 
über Gibbs, die — zwar nicht vollständig!) — 
das wachsende Interesse an dem großen Amerikaner 
zeigt. — Den reichen Inhalt schließt ein Register auf, 
in das leider nicht die von Gibbs hier und da falsch 
geschriebenen Personennamen seiner Korrespondenten 
und der Postlisten aufgenommen sind. 


Radebeul/Dresden. Rudolph Zaunick. 


1) Wir vermissen hier ungern u.a. A. Euckens aufschluß- 
reichen ‚‚Versuch, das Lebenswerk W. Gibbs, der Gegenwart nahe 
zu bringen“ (in: Naturwiss. Bd. 26 (1938), 8. 230—233 und W. 
Nernsts ‚‚Erinnerung an den 100, Geburtstag von Gibbs‘ 
(ebenda Bd, 27 (1939), S..393 f.) 


Dr.-Ing. R. C. Oldenbourg und Dr.-Ing. H. Sartorius, 
Dynamik selbsttätiger Regelungen. 
Band I. Allgem. und mathem. Grundlagen. Stetige 
und unstetige Regelungen. Nichtlinearitäten. 2. Aufl. 


“2588. mit 112 Abb. und 1 Tafel. München 1951. 
Verlag R. Oldenbourg. Preis geb. 26,— DM. 


Die erste Auflage dieses Buches erschien 1944 und 
war schnell vergriffen; 1948 wurde in den USA ein 
Nachdruck hergestellt. Beides beweist, wie not- 
wendig eine Darstellung der theoretischen Regelungs- 
technik benötigt wird, die mathematische Hilfsmittel 
verstärkt auszunutzen versteht. } 

Gleich zu Beginn wird mit Hilfe der Laplacetrans- 
formation der Zusammenhang zwischen den üblichen 
Beschreibungsformen von Regelgliedern mit Differen- 
tialgleichungen und Integraldarstellungen, sowie dem 
Frequenzgang hergestellt und in ingenieurmäßiger« 
Form das Rüstzeug zur Durchrechnung automatischer 
Regelaufgaben bereitgestellt. Im Hinblick auf die 
physikalische, anschauliche Bedeutung wird dabei be- 
sonders die Übergangsfunktion der Regelsysteme 
herangezogen. Für die Stabilitätsuntersuchungen aus 
dem Frequenzgang wird das Nyquist-Kriterium zur 


Verf gestellt. } 
Einen Überblick über den reichhaltigen Inhalt ver- 


mitteln einige Stichworte: | 

Mathematische Grundlagen, Stabilitätskriterien, 
Anwendung von rgangsfunktion und Frequenz- 

ang, Güte der Regelung, Stetige Regler (Elemente, 

tabilisierung, Optimaler Regelverlauf), Nichtlineari- 
täten (Ansprechempfindung, Ba und Lose), Un- 
stetige Regelungen, Anwendungsbeispeile. 

„Mit diesem reichhaltigen Inhalt und der besonderen 
Zielsetzung beim Einsatz der mathematischen Hilfs- 
mittel hat das Buch trotz aller Fortschritte der Regel- 
technik bis heute seinen Wert behalten. 

Zum Schluß soll auf einige Fehler hingewiesen 
werden: Die RE 38,3 für den ge- 
schlossenen . reis hat nur die Lösung z= 0 
(mit 17,7—8 als homogene Volterra-Gleichung 2. Art), 
e8 muß von 39,4 ausgegangen werden. Das bene 
Nyquist-Kriterium ist nur unter zusätzlichen Voraus- 
setzungen allgemein richtig. N.d. Lehmann, 


Die besprochenen und angezeigten. Bücher sind durch den Buchhandel zu beziehen, 
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EINGEGANGENE BÜCHER 


Bei der Schriftleitung sind folgende Bücher eingegangen (ausführliche Besprechung bleibt vorbehalten). 


Dr.-Ing. habil. Dr. Horst Müller, Mechanik für 
Ingenieure. (Westermanns Fachbücher der Inge- 
nierkunde), 254 S. mit 268 Abb. Braunschweig/ 
Berlin/Hamburg/Kiel 1951. Georg Westermann Ver- 
lag. Preis geh. 16,60 DM., geb. 18,60 DM. 


Hermann Hömberger, Einführung in die Elek- 

trotechnik. (Westermanns, Fachbücher der Elek- 

‚ trotechnik). 302 8. mit 277 Abb. Braunschweig 

De uchitg/Kiel 1951. Georg Westermann Ver- 
ag. 


: Walter Schnee, Über magische Quadräte und 
lineare Gitterpunktprobleme. (Berichte über die 
Verhandlungen der Sächsischen Akademie der Wissen- 


schaft zu Leipzig. Mathematisch-naturwissenschaft- 


liche Klasse, Bd. 98, Heft 1) 48 S. Berlin 1951. Aka- 
demie-Verlag. Preis brosch. 4,65 DM. 


; _ H. Goldstein, Ph. D., Classical Mechanics. 
0. XII + 398 8. mit 71 Abb. Cambridge (Mass.) 1951. 
N Addison Wesley Press. Preis geb. 7,— $. 


D. J. Struik, (Prof. a. Massachusetts Inst. of Tech- 
nology). Lectures on Classical Differential 
Geometry. VIII + 221 S. mit 120 Abb. Cambridge 
(Mass.) 1950. Addison Wesley Press. Preis geb. 


’ 


+ L. Landau und E. Lifshitz, The Classical Theory 
: of Fields (aus dem Russischen übersetzt von M. Ha- 
mermesh). IX + 354 S. mit 12 Abb. Cambridge 
'(Mass.) 1951. Addison Wesley Press. Preis geb. 

7,50 $. i a RT Te 

e o; r 

5 Dr. K. Federhofer (o. Prof. a. d. TH. Graz), Prü- 
Br fungs- und Übungsaufgaben aus der 
Mechanik des Punktes und des starren 
Körpers. I. Teil: KinematikundKine- 


+ 


” 


Te Viktor Blaess } Ra 


_ Kurz vor Vollendung seines 75. Lebensjahres, die 


1951 der emeritierte Professor der Mechanik an der 
Technischen Hochschule Darmstadt, Dr.-Ing. Viktor 
 Blaess. Fast durch 4 Jahrzehnte, von 1911 bis 1919 


licher Tiefe 


_ _ eram 12. September gefeiert hätte, starb am 16. August 


tikdes Punktes. 113 Aufgaben mit Lösungen. 
IV + 1038. mit 105 Abb. III. Teil: Kinematik 
umd Kinetikstarrer Systeme. 149 Auf- 
gaben mit Lösungen. V-+1398. mit 191 Abb, 
Wien 1951. Springer-Verlag. Preis je 9,60 DM. 


Dr. L.F. Ryan, Experiments on Aero- 
dynamic Cooling und M. Degen, Unter- 
suchungen an einem Gegenlaufpro- 
pellerimWindkanal. (Mitteilungen aus dem 
Institut für Aerodynamik an der ETH. Zürich, Nr. 18) 
65 S. mit 47 Abb. Zürich 1951. Verlag Leemann. 


R. Richter (em. Prof.a.d. TH. Karlsruhe), Elek- 
trischeMaschinen,BandI:Allgemeine 
Berechnungselemente. Die Gleich- 
strommaschine. 2.verbesserte Auflage. XVI+ 
630 S. mit 453 Abb. Basel 1951. Verlag Birkhäuser. 
Preis brosch. 45,25 Fr., geb. 49,40 Fr. 


Dr. H. von Sanden (o. Prof. a. d. TH. Hannover), - 
PraktischeMathematik (Teubners mathe- 
matische Leitfäden Bd.44). 2. erweiterte Auflage. 
120 S. mit 25 Abb. Leipzig 1951. B.G. Teubner, 
Verlagsgesellschaft. Preis kart. 3,80 DM. 


F. Lösch-F. Schoblik, Die Fakultät. (Ga’m’ma - 
funktion)undverwandteFunktionen 
mit besonderer Berücksichtigung 
ihrer Anwendungen. (Bearbeitet von Dr.EF. 
Lösch, Prof. a. d. Techn. Hochschule Stuttgart.) _ 


“ VI+ 2058. mit 22Abb. Leipzig 1951. Teubners- 


Verlagsgesellschaft. Preis geb. 16,80 DM. 


Dr.L.Bieberbach (ehem. Prof. a. d. Univ. Berlin). 
Einführungindie Funktionentheo- 
rie. Zweite neubearbeitete Aufl. 220 S. mit 43 Abb. 
Bielefeld 1952. Verlag für Wissenschaft und Fach- 
buch. Preis 12,60 DM. 
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h zurückzuführen. Genannt sei hierbei sein in der ZaMM 


1926 erschienener Aufsatz über den ‚‚Massenausgleich 
rasch umlaufender Körper“. Bekannt ist auch der von 
ihm entwickelte mechanische Resonanz-Schwingungs- 
messer. 


Von seinen bis zur letzten Anschaulichkeit vordrin- 


‚genden Bemühungen legen die vielen von ihm konstru- 


ierten Modelle, z. B. des ebenen Spannungstensors und 
solche zur Darstellung von Schwingungserscheinungen 
beredtes Zeugnis ab. Sie vertieften den Vorlesungsstoff _ 
ungemein und gestalteten seine Darlegungen für seine 


vielen Schüler zu einem unvergeßlichen Erlebnis. 


Weiterhin fand Blaess im Jahre 1934 ein beson- hy 
_ deres Verfahren zur angenäherten Lösung von Diffe- 


rentialgleichungen, das als das. Blaesssche Verfahren 


Runge-Kutta-N yström das der Praxis wohl am 


besten zugängliche ist. AR 

> Von seinen Arbeiten auf anderen Gebieten sei hier 
nur noch sein 1911 im Verlag Oldenbourg erschienenes 

Buch über die Berechnung verzweigter Rohrleitungen 

für Luft und Gase erwähnt, das besonders damals von _ 

17 großer Wichtigkeit war. 

in Wer das Glück hatte, Prof. Blaess trotz seiner 
etwas schwer zugänglichen Wesensart persönlich näher 

kennen zulernen, dem offenbarte sich ein selten lauterer 
‚ Charakter, ein freundlich gütiger und liebenswürdiger 
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Mensch, dessen wunderbarem Einfluß man sich nicht 
entziehen konnte. Mehr als sonst ist bei Blaess der 
Grundsatz ‚‚mehr sein als scheinen‘ verwirklicht ge- 
"wesen. 


Darmstadt. K. Karas. 


IV. Kongreß der Italienischen Mathematiker-Ver- 
einigung in Taormina vom 25. bis 380. Oktober 1951. 


. Unter dem Protektorate des Präsidenten der Ita- 
lienischen Republik hatte sich ein Ehrenkomitee von 
18 Persönlichkeiten aus den Spitzen der Verwaltung, 
der Industrie und der Wissenschaft gebildet, das die 
materiellen Vorbedingungen für das Zustandekommen 
dieses} Kongresses geschaffen hat. Der alle Erwar- 
tungen überschreitende Verlauf wird allen, Teil- 
nehmern unvergeßlich bleiben. Es waren über 300 
italienische und mehr als 30 auswärtige Mathematiker 
erschienen, die sich zu einer Einheit zusammenfanden, 
die in jedem Sinne als eine Keimzelle störungsfreier 
internationaler Zusammenarbeit angesehen werden 
kann. Außer einer Anzahl von allgemeinen Vorträgen 
in allgemeinen Sitzungen wurden in 5 Sektionen: 
Analysis, Geometrie, Mechanik u. math. Physik, Wirt- 
schaftsmathematik, Philosophie und Didaktik über 
100 Vorträge geboten. Durch die auswärtigen Teil- 
nehmer waren die folgenden Staaten vertreten: 
Belgien (2), Frankreich (6), Holland (1), Österreich (8), 
USA (3), Schweden (1), Spanien (2), Deutschland (8), 
Polen (2). Der Wunsch, daß sich auch die übrigen 
Staaten der Welt auf dem Boden der mathematischen 
Wissenschaften zusammenfinden mögen, blieb leider 
diesmal noch unerfüllt, doch konnte hier an einem 
Beispiel ersehen werden, daß eine solche Vereinigung 
durchaus im Bereich der Möglichkeiten liegt und auch 
über die fachlichen Interessen hinaus den Boden zu 
einer allgemeinen Verständigung der Menschheit 
bilden könnte. Was das Land Italien hier geboten hat, 
lag nicht nur auf dem Boden der Mathematik, sondern 
auch auf dem der Kultur und der Geschichte des 
Landes, die den einzigartigen Rahmen für eine Ver- 
anstaltung abgaben, der sonst kaum irgendwo in der 
Welt etwas Ähnliches an die Seite gestellt werden 
kann. Th. Pöschl, Karlsruhe. 

Im Anschluß an den Kongreß der Unione Mate- 
matica Italiana (25.—30. 10. 51 in Taormina) folgten 
Prof. Dr. L.Collatz, Hannover, und Prof. Dr. H. 
Görtler, Freiburg i. Br., am 2. und 3. 11.51 einer 
Einladung zu Vorträgen am Mathematischen Institut 
der Universität Triest. 


Tagung der Deutschen Mathematiker - Vereinigung 
in Berlin vom 16. bis 21. September 1951. 

Die Tagung war mit rund 200 Teilnehmern wohl 
die am stärksten besuchte, welche die DMV je ver- 
anstaltet,hat. Eine Reihe ausländischer Gäste konnte 
begrüßt werden, darunter auch Kollegen, die früher 
in Berlin tätig waren und nun wieder an der Tagung 
der DMV teilnahmen, 

Rund 50 Vorträge legten Zeugnis ab von dem 
lebendigen Leben und den Fortschritten in der mathe- 
matischen Forschung. Einstündige Hauptvorträge 
hielten die Herren 
Kähler-Leipzig: Die Einheitlichkeit der mathe- 

matischen Wissenschaften. 

R.v. Mises-Cambridge, USA: Die Theorie der 
„Statistischen Funktionen. 
H. Hopf-Zürich: Vom Bolzano-Kroneckerschen 

Satz zur Homotopietheorie der Sphären. 
Eichler-Münster: Arithmetische Theorie der 

quadratischen Formen und Thetafunktionen. 
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Außerdem wurden Vorträge aus folgenden Gebieten 
gehalten: Komplexe Funktionentheorie, Geometrie 
einschließlich Diff.-Geom. und algebr. Geom., Inte- 
gralgleichungen, partielle und gewöhnliche Diff.- 
Gleichungen und Rand-- und Eigenwertaufgaben, 
Algebra, Gruppentheorie, Zahlentheorie, Summierungs- 
verfahren, Mengenlehre und andere Einzelfragen. - 

Kamke. 


Ein internationales Symposion über nicht- 
lineare Schwingungen fand vom 18. bis 22. Sept. 
1951, veranstaltet von der IUTAM*) und unterstützt 
von der Unesco, auf der Ile de Porquerolles (im Mittel- 
meer bei Toulon) statt. Aus 12 Ländern nahmen 35 
einzeln Eingeladene teil, darunter aus Deutschland 
(je mit einem Vortrag) die Herren Collatz, Grammel 
und Mettler. Präsident war Professor P6rös- 
Paris, die Organisation lag in den Händen eines 
Komitees mit Professor Koiter - Delft als 
Sekretär. Diesem Organisationskomittee gehörte 
deutscherseits Herr Grammel an. Die 18 Vorträge 
umfaßten alle Teile des in den letzten zwei Jahr- 
zehnten stark weiterentwickelten Gebietes der nicht- 
linearen Schwingungen von einem bis zu unendlich 
vielen Freiheitsgraden, und zwar von der abstrakt 
mathematischen Theorie bis zu ihren Anwendungen. 
Der persönliche Kontakt zwischen den Teilnehmern 
war sehr eng und herzlich und wurde durch den 
genius loci dieser schönen, abseits von Lärm und Ver- 
kehr liegenden Insel gefördert. Ein Kongreßbericht 
mit ausführlichen Vortragsauszügen wird in Frank- 


reich erscheinen. 
# 


2) International Union of the theoretical and applied 
Mechanics, 


Stuttgart: Der o. Prof. an der Techn. Hochschule 
Aachen Dr. Günther Schulz hat einen Ruf auf den 
durch die Emeritierung von Prof. Dr. Friedrich 
Pfeiffer frei gewordenen Lehrstuhl für Mathematik 
an der Techn. Hochschule Stuttgart angenommen. 


München: Am 4.Dez. 1951 verstarb im 
90. Lebensjahre Geh.-Rat Prof. Dr. Sebastian 
Finsterwalder. Er war seit 1891 o. Prof.a.d. 
Techn. Hochschule München (em. 1932) und seit 1903 
Mitglied der Bayrischen Akademie der Wissenschaften. 


Der apl. Prof. WilhelmMaak wurdeals o. Prof. 
der Mathematik an die Univ. München berufen. 


Dresden: An der Techn. Hochschule Dresden 
habilitierte sich Dr. N.J. Lehmann für das Fach 
der angew. Mathematik. 


Wien: Inder Zeit vom 9. bis 14. Sept. 1952 findet 
in Salzburg der 3.Österreichische Mathematiker- 


kongreß statt, der wieder als internationales Mathe- 
matikertreffen gestaltet werden soll. Für seine 
wissenschaftl. Arbeit werden folgende Sektionen ge- 
bildet: Analysis, Geometrie und Topologie, Algebra 
und Zahlentheorie, angew. Mathematik, Geschichte 
und Philosophie, 


Bloomington: Von Januar 1952 ab wird eine 
neue Zeitschrift ‚Journal of Rational Mechanics and 
Analysis“ von den Herren T. Y. Thomas und 
C, Truesdell vom Graduate Institute for Applied- 
Mathematics der Indiana Universität in Blooming- 
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